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Abstract 

We  formulate  and  test  a  mathematical  program  to  select  Air  and  Sea  Ports  of  De¬ 
barkation  and  intermediate  logistical  distribution  centers,  through  which  we  route 
military  supplies  over  a  directed  transportation  network  to  meet  aggregated  weekly 
demands  by  military  units  conducting  a  steady  state  contingency  operation.  The 
multi-objective  model  seeks  to  minimize  a  weighted  combination  of  the  total  risk  en¬ 
countered  by  transported  supplies,  the  total  distance  traveled  by  supplies,  and  the 
maximum  per  capita  workload  supported  by  transportation  assets  at  a  given  echelon 
(i.e.,  port-to-distribution  center  versus  distribution  center  to  demands).  Within  our 
formulation,  we  account  for  capacities  on  arc  flows  and  node  throughputs,  with  the 
latter  enabling  the  representation  of  material  handling  equipment  limitations  at  the 
ports  and  distribution  centers.  For  our  model,  we  develop  and  test  an  Excel-based 
decision  support  tool  that  invokes  the  commercial  solver  CPLEX  (Version  12.5)  to 
solve  the  underlying  mixed-integer  linear  program,  and  we  demonstrate  its  efficacy  on 
a  representative  instance.  For  this  instance,  we  identify  extreme  points  and  selected 
interior  solutions  on  the  Pareto  efficient  frontier  and  examine  the  model’s  sensitivity 
to  selected  parameters.  We  conclude  by  discussing  how  the  model  can  account  for 
intra-theater  airlift  and  outline  modifications  that  can  account  for  expected  pilferage 
losses  within  the  distribution  system. 

Key  words:  Distribution  Network,  Facility  Location,  Supply  Chain  Design 
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THE  MILITARY  THEATER  DISTRIBUTION 


NETWORK  DESIGN  PROBLEM 

I.  Introduction 

“Amateurs  talk  about  tactics,  but  professionals  study  logistics.” 

Gen  Robert  Barrow,  USMC  [28] 

This  chapter  explains  the  background  for  this  research.  First,  in  Section  1.1  we 
describe  the  significance  of  logistics  in  historical  campaigns  and  provide  an  overview 
of  the  Army’s  distribution  system.  Next,  in  Section  1.2  we  formally  define  the  problem 
and  present  the  problem  statement.  In  Section  1.3  we  outline  the  research  objectives. 
Finally,  in  Section  1.4  we  provide  an  overview  of  the  remaining  chapters  in  this  thesis. 

1.1  Overview 

The  purpose  of  this  research  is  to  facilitate  the  movement  of  Army  cargo/materiel 
from  large  depots  through  a  distribution  system  to  the  user  level  where  it  is  needed 
for  warfighting.  The  success  or  failure  of  military  campaigns  can  provide  case  studies 
in  logistics.  Examples  are  Hannibal’s  crossing  of  the  Alps,  Napoleon’s  march  towards 
Moscow,  Sherman’s  March  to  the  Sea  during  the  American  Civil  War,  and  the  Al¬ 
lied  Armies’  movement  to  Berlin  after  the  successful  landings  on  Normandy  beaches. 
More  modern  examples  include  the  humanitarian  response  to  the  2010  Haiti  earth¬ 
quake  and  the  International  Security  Force  Assistance  (ISAF)  movement  into  the 
most  remote  villages  of  the  Hindu  Kush.  The  success  of  the  United  States  military  in 
the  field  is  purposefully  reliant  on  the  joint  interdependence  each  service  has  on  the 
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other  services’  capabilities.  Army  Doctrine  Publication  3-0  (ADP  3-0)  Unified  Land 
Operations  explains  the  Army’s  involvement  in  the  joint  environment.  This  reliance 
encompasses  functions  such  as  intelligence,  operations,  and  logistics,  and  it  allows  the 
Army  to  focus  efforts  on  its  core  tasks  rather  than  expend  effort  on  tasks  the  other 
services  are  well  suited  to  perform.  For  example,  the  United  States  Air  Force  (USAF) 
provides  “worldwide  cargo  and  personnel  airlift,  air  refueling,  and  aeromedical  evacu¬ 
ation,”  whereas  the  United  States  Navy  (USN)  conducts  “over-the-shore”  operations 
to  up-load,  transport,  and  down-load  cargo”  [15].  As  the  USAF  and  USN  conduct 
these  core  tasks  of  their  respective  service,  the  Army  is  better  able  to  focus  on  its 
core  tasks.  However,  although  joint  interdependence  exists  by  design,  Title  10,  U.S. 
Code  requires  that  each  service  retain  responsibility  for  the  sustainment  of  forces  it 
allocates  to  a  joint  force  [15]. 

The  Department  of  Defense’s  (DOD)  Global  Distribution  process  is  a  joint  endeavor 
in  which  each  service  contributes  to  distribute  cargo/materiel  to  meet  the  demands 
of  the  warfighter.  The  United  States  Transportation  Command  (USTRANSCOM) 
is  responsible  for  planning,  resourcing  and  executing  the  Global  distribution  process 
which  is  the  mechanism  USTRASNSCOM  uses  to  move  cargo/materiel  from  points  of 
origin  to  an  aerial  or  sea  port  of  debarkation  (A/SPOD)  [7].  The  last  segment  of  the 
transport  of  cargo/materiel  is  called  intra-theater  movement  and  is  the  responsibility 
of  each  service.  This  thesis  considers  the  intra-theater  distribution  of  cargo/material 
by  Army  systems  to  support  Army  forces  after  the  supplies  arrive  at  an  A/SPOD. 
Cargo/materiel,  referred  to  as  classes  of  military  supplies  in  Army  Doctrine  Publi¬ 
cation  4-0  (ADP  4-0)  and  shown  in  Figure  1,  will  heretofore  be  referred  to  as  supplies. 
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Class 

Item 

1 

Subsistence  (food)  and  gratuitous  (free)  health  and 
comfort  items 

II 

Clothing,  individual  equipment,  tentage,  organizational 
tool  sets  and  kits,  hand  tools,  unclassified  maps, 
administrative  and  housekeeping  supplies,  and 
^tjuigmen^ 

III 

Petroleum,  oil  and  lubricants  (package  and  bulk): 
petroleum,  fuels,  lubricants,  hydraulic  and  insulating  oils, 
preservatives,  liquids  and  gases,  bulk  chemical  products, 
coolants,  deicer,  antifreeze  compounds,  components, 
additives  of  petroleum  and  chemical  products,  and  coal 

IV 

Construction  materials,  including  installed  equipment  and 
all  fortification  and  barrier  materials 

V 

Ammunition  of  all  types:  bombs,  explosives,  mines, 
fuzes,  detonators,  pyrotechnics,  missiles,  rockets, 
propellants,  and  associated  items 

VI 

Personal  demand  items  (such  as  health  and  hygiene 
products,  soaps  and  toothpaste,  writing  material,  snack 
food,  beverages,  cigarettes,  batteries,  and  cameras — 
nonmilitarv  sales  items] 

VII 

Major  end  items  such  as  launchers,  tanks,  mobile 
machine  shops,  and  vehicles 

VIII 

Medical  materiel  including  repair  parts  peculiar  to  medical 
equipment 

IX 

Repair  parts  and  components  to  include  kits,  assemblies, 
and  subassemblies  (repairable  or  nonrepairable)  required 
for  maintenance  support  of  all  equipment 

X 

Material  to  support  nonmilitary  programs  such  as 
agriculture  and  economic  development  (not  included  in 
Classes  Uhrou2h[X} 

Misc 

Water,  salvage,  and  captured  material 

Figure  1.  US  Military  Classes  of  Supply 


The  Army’s  concept  of  sustainment  involves  moving  and  supporting  forces.  As  ex¬ 
plained  in  ADP  4-0,  sustainment  is  built  around  three  elements;  logistics,  personnel 
services,  and  health  services,  as  shown  in  Figure  2: 


Figure  2.  Army  Sustainment  Operations 


The  distribution  of  supplies  occurs  within  the  logistics  field.  Figure  3  from  ADP 
4-0  displays  the  tasks  Army  logisticians  plan  and  execute.  This  thesis  addresses  the 
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distribution  of  supplies  from  an  A/SPOD  through  logistic  hubs  to  meet  the  demands 
of  Army  maneuver  units  in  the  held. 


Figure  3.  ArmyLogistics  Operations 


In  order  to  fight  effectively  in  the  held,  the  Army  must  conduct  unified  land  oper¬ 
ations,  a  term  which  describes  how  the  Army  fights  simultaneously  in  the  offense, 
defense,  and  stability  operations.  To  conduct  unified  land  operations  the  Army  must 
sustain  itself  to  enable  operational  reach  (operating  over  long  distances  from  an  A/S¬ 
POD),  freedom  of  action  (the  depth  and  breadth  of  the  operational  area),  and  pro¬ 
longed  endurance  (extended  durations)  [16]. 

The  remainder  of  this  section  describes  the  movement  of  military  supplies  through  hve 
echelons  from  the  A/SPOD  to  the  maneuver  companies  in  the  held  in  order  to  sup¬ 
port  unified  land  operations.  However,  this  research  will  only  consider  the  movement 
of  supplies  from  the  A/SPOD  to  brigade  combat  teams  (BCTs)  at  forward  operating 
bases  (FOBs).  The  hve  echelons  from  highest  to  lowest  are  the  theater  sustainment 
command  (TSC),  the  expeditionary  sustainment  command  (ESC),  the  sustainment 
brigade  (SB),  the  brigade  support  battalion  (BSB),  and  the  forward  support  company 
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(FSC). 


The  highest  echelon  is  the  Army’s  TSC,  which  is  a  fixed  headquarters  that  provides 
operational-level  support  to  Army  forces  or  a  joint  task  force  (JTF).  ADP  4-0  explains 
the  TSCs  organizational  structure  and  a  diagram  is  shown  in  Figure  4.  The  TSC  can 
support  one  or  more  SBs  as  well  as  one  or  more  ESCs.  The  TSC  is  responsible  for 
distributing  supplies  from  A/SPODs  to  the  next  echelon,  the  ESC. 


The  second  echelon,  the  ESC,  is  a  headquarters  which  may  either  augment  a  TSC 
in  the  event  of  a  very  large  and  enduring  operation  or,  more  typically,  replace  it  for 
a  smaller-scale  or  short-term  operation.  An  ESC  is  established  when  multiple  SBs 
are  required  to  support  the  GCC,  or  the  theater  of  operations  is  far  away  from  the 
TSC  such  that  an  additional  command  node,  the  ESC,  is  necessary  to  facilitate  the 
distribution  of  supplies.  The  ESC  then  becomes  a  headquarters  to  SBs. 
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Figure  4.  TSC  Organizational  Chart 


The  SB  is  the  next  organization  that  is  responsible  for  receiving  and  moving  mili- 


5 


tary  supplies  through  the  distribution  system.  The  SB  supports  division-sized  units 
(  20,000  soldiers).  Like  the  TSC  and  ESC,  the  SB  is  modular,  and  its  design  varies 
depending  on  the  needs  of  the  GCC.  Generally  speaking,  the  SB  will  have  a  spe¬ 
cial  troops  battalion  (STB)  and  one  or  more  combat  sustainment  support  battalions 
(CSSBs),  as  shown  in  Figure  4.  The  physical  transport  of  supplies  from  the  TSC/ESC 
to  the  SB  is  accomplished  by  CSSBs.  CSSBs,  displayed  in  Figure  4  are  tasked  orga¬ 
nized  with  various  types  of  transportation  untis  based  on  the  operational  requirements 
to  move  supplies  within  these  echelons. 

As  explained  in  ADP  4-0,  the  next  organization  responsible  for  movement  of  military 
supplies  in  the  distribution  system  is  the  BSB,  and  a  representative  diagram  is  shown 
in  Figure  5.  The  BSB  is  assigned  to  the  maneuver  BCT  both  in  a  deployed  and  non- 
deployed  environment.  The  relationship  between  the  SB  and  the  BSB  is  maintained 
through  staff  coordination  for  the  requisition  and  movement  of  supplies  between  the 
two  organizations.  The  SB’s  support  operations  officer  (SPO)  coordinates  with  the 
BSB  support  operations  officer  (SPO)  and  the  maneuver  brigade  S3  to  move  military 
supplies  to  the  FSCs,  the  final  echelon  in  the  distribution  of  military  supplies.  The 
other  units  depicted  in  Figure  5  provide  the  full  complement  of  support  necessary 
for  a  BSB  to  perform  the  logistics  operations  depicted  in  Figure  3.  From  the  BSBs, 
supplies  are  transported  using  internal  unit  haul  assets  to  the  last  echelon,  the  FSCs, 
each  of  which  supports  a  battalion  -sized  unit  (~600  soldiers). 
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Figure  5.  BSB  Organizational  Chart 


Figure  6  from  ADP  4-0  shows  the  FSC  organization  which  moves  supplies  to  the 
maneuver  companies.  The  FSCs  are  assigned  to  the  BSB  and  are  under  operational 
control  (OPCON)  of  the  supported  maneuver  battalion.  The  FSCs  move  military 
supplies  directly  to  the  maneuver  companies.  The  other  units  depicted  in  Figure  6 
provide  the  full  complement  of  support  necessary  for  a  FSC  to  perform  the  logistics 
operations  depicted  in  Figure  3. 


Figure  6.  FSC  Organizational  Chart 
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This  research  considers  only  the  flow  of  military  supplies  from  the  TSC/ESC  echelon 
through  the  SB  echelon  and  to  the  BSB  locations.  In  terms  of  physical  locations, 
this  movement  comprises  the  flow  of  supplies  from  A/SPODs  through  intermediate 
distribution  centers  to  Forward  Operating  Bases  (FOBs).  Underlying  this  research 
scope  is  the  assumption  that  only  a  TSC  or  an  ESC  is  utilized  but  not  both;  however, 
the  model  we  propose  may  be  readily  adapted  if  the  assumption  does  not  hold. 

1.2  Problem  Definition  and  Problem  Statement 

Army  planners  use  experience,  training  center  exercises,  and  models  to  plan  real 
world  deployments.  The  Army’s  Training  and  Doctrine  Command  Analysis  Center 
at  Ft.  Lee,  Virginia  (TRAC-LEE)  is  an  analysis  agency  that  provides  planners  with 
support  in  planning  sustainment  operations.  This  research  was  initiated  to  comple¬ 
ment  TRAC-LEE’s  support  of  LInited  States  Army  Africa  (USARAF)  to  enhance  its 
logistic  planning  capabilities,  but  it  could  also  be  used  by  other  major  commands 
(MAJCOM).  TRAC-LEE  employs  a  tool  called  Logistics  Battle  Command  (LBC) 
that  “has  planning  and  decision  support  features  to  enable  a  simulated  logistics  bat¬ 
tle  captain  to  (1)  monitor  the  logistics  common  operating  picture,  (2)  forecast  demand 
for  most  classes  of  supply,  and  (3)  initiate  and  adjust  missions  to  distribute  supplies 
and  perform  sustainment  functions”  [26].  LBC  is  a  dynamic  stochastic  model  that 
uses  demands  and  nodes  to  assess  the  utility  of  A/SPOD  and  intra-theater  logistic 
hub  locations.  TRAC-LEE  would  be  better  able  to  assist  USARAF  if  it  was  possible 
to  quickly  model  a  logistics  infrastructure  system  that  supports  operations  within  its 
AOR.  This  research  will  explore  ways  to  enhance  LBC  by  developing  good  starting 
solutions  based  on  analytic  rigor. 


In  support  of  steady  state  supply  operations  for  a  military  contingency  operation  and 


given  an  existing  road  network,  a  set  of  demands  for  military  supplies,  and  a  finite 
number  of  logistical  hubs  at  the  various  echelons,  we  seek  to  locate  logistical  hubs 
and  route  the  flow  of  military  supplies  to  provide  effective  logistical  support.  Within 
this  context,  we  consider  alternative  means  to  measure  effective  logistical  support,  to 
include  minimizing  the  total  risk  encountered  by  transporting  supplies,  minimizing 
the  total  distance  traveled  by  supplies,  balancing  the  work  load  between  echelons  of 
command  or  the  span  of  control  at  each  mission  command  (MC)  node.  In  order  to 
adhere  to  doctrinal  implementation  of  logistical  support  operations,  we  restrict  our 
solutions  based  on  the  assumptions  that  material  supplies  must  transit  the  established 
MC  echelons  (skipping  levels  is  not  allowed),  the  span  of  control  for  a  headquarters 
(HQ)  is  2-5  subordinate  units,  and  the  model  is  built  to  support  steady-state  opera¬ 
tions. 

Given  the  preceding  doctrinal  explanation  for  the  movement  of  supplies  through  the 
Army  distribution  system,  the  following  is  an  operational  explanation  of  the  move¬ 
ment  of  supplies  that  will  be  used  in  this  research.  A  TSC  or  ESC,  but  not  both,  will 
be  used.  The  TSC/ESC  will  be  the  coordinating  HQ  for  the  movement  of  supplies 
from  CONUS  into  the  theater.  A  SB  will  be  co-located  at  an  A/SPOD,  and  supplies 
will  be  moved  to  regional  warehouse  distribution  centers  by  CSSB  assets.  Hereafter, 
we  use  the  term  regional  warehouse  distribution  center  (DC)  in  lieu  of  the  colloquial 
‘logistical  hubs’  to  represent  locations  established  to  support  the  movement,  storage, 
repackaging,  and  transshipment  of  military  supplies.  The  DCs  are  an  extension  of 
the  SBs  and  are  controlled  by  the  SB  at  the  A/SPOD.  In  some  instances  the  CSSB 
may  move  supplies  directly  from  the  SB  at  an  A/SPOD  to  meet  demands  at  a  FOB 
where  a  BSB  is  located.  However,  in  our  model  this  is  not  allowed  as  using  DCs 
allows  us  to  seek  more  efficient  methods  moving  supplies  from  the  A/SPOD  to  the 
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FOB.  The  DC  is  controlled  by  a  movement  control  team  (MCT)  from  the  SB,  and 
tailored  supply  packages  to  meet  BCT  demands  at  a  FOB  are  generated  by  material 
handling  teams  (MHE)  provided  by  the  SB.  Once  the  tailored  supply  packages  are 
created,  they  are  moved  to  a  FOB  by  CSSB  assets.  Finally,  this  research  will  use 
DCs,  but  the  warehousing  aspect  of  the  distribution  system  will  not  be  considered. 

With  the  above  discussion,  the  following  problem  statement  is  proposed.  Given  an 
existing  capacitated  road  network,  seaports,  and  airfields;  a  geographically  dispersed 
set  of  demands  for  supplies  (i.e.,  BCT-sized  FOB  locations);  a  bounded  number  of 
TSC/ESC-operated  APODs  and  SPODs;  and  a  bounded  number  of  SB-operated 
intermediate  distribution  centers;  we  seek  to  determine  good  solutions  for  locating 
A/SPODs  and  distribution  centers,  as  well  as  the  routing  of  supplies  over  the  road 
network,  in  order  to  meet  demand.  Because  ‘good’  is  an  ill-defined  term,  we  consider 
three  alternative  metrics  to  minimize,  either  individually  or  in  combination:  the  ac¬ 
cumulated  risk  to  supplies  being  transported,  the  total  distance  traveled  by  supplies, 
and  the  maximum  per  capita  workload  supported  by  transportation  assets  at  a  given 
echelon  (i.e.,  A/SPOD-to-distribution  center  versus  distribution  center  to  demands). 


1.3  Research  Objectives 

There  are  three  proposed  research  objectives: 

1.  To  develop  a  math  programming  model  that  identifies  capacitated  A/SPOD  lo¬ 
cations  and  intermediate  DCs  in  order  to  maximize  multi-commodity  flows  over 
existing  in-theater  infrastructure  networks  in  support  of  contingency  operations. 

2.  To  solve  the  model  by  enabling  a  simple-to-use-and-understand  user-input  for  crit¬ 
ical  information  about  A/SPODs  that  may  be  used. 
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3.  To  demonstrate  the  model  by  using  a  scenario-based  case  study. 

1.4  Overview  of  Remaining  Chapters 

Chapter  2  reviews  work  that  has  already  been  published  on  the  topic  of  this  thesis 
-  facility  location  multicommodity  transshipment  problems.  A  detailed  literature 
search  of  this  topic  will  help  narrow  the  scope  of  the  research  and  provide  clarity  of 
the  research  objectives  and  question.  In  this  way  it  will  be  shown  that  the  work  in  this 
thesis  is  original.  Chapter  3  presents  our  modeling  assumptions,  introduces  notation, 
and  formally  presents  our  mathematical  program  with  several  variants.  It  further 
addresses  the  limitations  of  our  modeling  approach.  Chapter  4  provides  a  description 
of  the  model  validation  using  a  scenario-based  case  study.  We  conclude  Chapter  4 
with  a  presentation  and  comparison  of  the  solutions  yielded  by  our  modeling  variants. 
Chapter  5  summarizes  the  major  contributions  of  this  research,  addresses  some  of  its 
limitations,  and  proposes  directions  in  which  to  continue  and  extend  the  research  for 
even  greater  utility. 
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II.  Literature  Review 


“Gentlemen,  the  officer  who  doesn’t  know  his  communications 
and  supply  as  well  as  his  tactics  is  totally  useless.” 

Gen.  George  S.  Patton,  USA  [25] 


2.1  Overview 

This  chapter  provides  a  review  of  literature  related  to  selected  optimization  modeling 
research.  The  first  part  of  this  chapter  provides  and  discusses  problem  formulations 
for  several  selected  facility  location  models  including  the  set  covering  problem,  the 
maximal  covering  location  problem,  the  p-center  problem,  and  the  p-median  problem. 
The  second  part  of  this  chapter  provides  and  discusses  problem  formulations  for  sev¬ 
eral  common  network  flow  models  including  the  shortest  path  problem,  the  maximum 
flow  problem,  and  the  minimum  cost  flow  problem.  Three  additional  network  flow 
models  are  discussed  as  well:  the  transportation  problem,  the  assignment  problem, 
and  the  transshipment  problem.  The  last  section  is  a  review  of  the  facility  location 
and  network  flow  problem.  This  review  is  not  exhaustive;  rather,  it  provides  the  back¬ 
ground  to  highlight  and  understand  the  work  that  has  previously  been  accomplished 
in  these  areas. 


2.2  Facility  Location  Models 

In  general,  facility  location  is  a  process  whereby  a  commodity  is  placed  in  a  specific 
location  because  the  commodity  provides  an  advantage  to  the  decision  maker.  The 
location  decision  is  based  on  conditions  in  the  environment  such  that  the  the  most 
benefit  is  derived.  For  example,  early  settlers  in  the  eastern  United  States  established 
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their  communities  in  locations  that  provided  the  most  benefit  in  terms  of  physical  pro¬ 
tection,  water  sources,  agricultural  opportunities,  and  transportation  systems.  Today, 
retail  companies  build  distribution  centers  and  stores  in  locations  that  meet  demands 
and  maximize  their  profits.  When  location  problems  become  complex,  as  in  the  case 
of  designing  a  supply  chain  for  a  contemporary  retail  company,  mathematical  pro¬ 
gramming  models  can  suggest  solutions. 

One  of  the  first  to  consider  facility  location  models  was  Hakimi  [22],  He  analyzed 
a  network  from  the  perspective  of  the  most  geographically  centered  node.  From  the 
central  node  he  calculated  the  distance  to  each  node  in  the  network  and  summed  the 
distances.  He  then  designated  an  adjacent  node  as  the  “central”  node  and  repeated 
the  process  until  all  nodes  in  the  network  had  been  selected  as  the  “central”  node.  He 
compared  all  the  iterations  and  noticed  that  the  shortest  summation  was  not  alwys 
the  node  in  the  geographic  center.  His  work  essentially  became  what  is  today  defined 
as  the  p-median  problem. 

Hakimi’s  successors  (Geoffrion  &  Graves  [20]  and  Telia  &  Lee  [37])  noted,  that  by 
moving  away  from  the  center  of  the  network  to  other  locations  within  the  network,  it 
was  possible  to  maximize  and  minimize  different  measures  of  demand.  By  consider¬ 
ing  alternative  assumption-based  objectives  and  constraints,  several  facility  location 
models  were  developed.  In  this  literature  review,  we  will  address  the  set  covering  loca¬ 
tion  problem  (SCLP),  the  maximal  covering  location  problem  (MCLP),  the  p-median 
problem  (PMP),  and  the  p-center  problem  (PCP).  While  each  model  can  be  used  to 
identify  an  optimal  solution  to  a  problem,  each  differs  with  regard  to  its  underlying 
assumptions  and  constraints. 
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In  a  facility  location  model,  the  feasible  region  for  facility  location  can  be  modeled 
as  either  discrete  or  continuous.  A  map  of  an  area  that  a  distribution  network  cov¬ 
ers  will  depict  rivers,  lakes,  roads,  road  junctions,  railroad  tracks,  built  up  areas, 
and  other  map  features.  In  a  ground  transportation  distribution  network,  the  road 
junctions  are  nodes.  In  a  discrete  facility  location  problem  there  are  a  finite  number 
of  locations  to  locate  a  facility.  In  other  words,  a  facility  can  be  located  only  at  a 
node.  In  a  continuous  facility  location  problem,  a  facility  can  be  located  in  one  of  two 
ways.  First,  it  can  be  restricted  to  any  node  or  anywhere  along  an  arc  (a  retailer  for 
a  chain  of  stores).  Secondly,  if  the  network  structure  is  considered  to  be  irrelevant, 
then  the  facility  can  be  located  anywhere  in  M2  (a  cell  phone  tower).  This  research 
will  consider  only  discrete  models. 

The  purpose  of  a  facility  in  a  supply  chain  network  is  to  ship,  store,  or  receive  goods. 
The  capabilities  of  a  given  facility  affect  problem  formulation  and  ultimately  the  op¬ 
timal  solution.  A  source  node  sends  goods  through  the  network  to  a  sink  node  which 
requires  goods.  All  nodes  between  the  source  node  and  the  sink  node  are  potential 
facility  locations  with  the  capability  to  temporarily  store  goods.  In  an  uncapacitated 
facility  location  problem  (UFLP),  the  upper  temporary  storage  limit  is  determined 
“without  any  budgetary,  technological,  or  physical  restrictions.”  Alternatively,  the 
capacitated  facility  location  problem  (CFLP)  does  incorporate  an  upper  temporary 
storage  limit.  In  both  the  UFLP  and  the  CFLP  each  node’s  temporary  storage  limit 
may  not  necessarily  be  the  same  [41],  This  research  will  incorporate  aspects  of  CFLP, 
in  modeling. 

Sometimes  in  a  network  there  is  a  reason  to  give  preference  to  the  use  of  a  facility 
over  others,  or  preference  to  the  service  of  a  given  demand  over  other  demands.  Per- 
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haps  using  a  specific  facility  will  result  in  a  tax  credit  or  a  public  relations  advantage. 
Maybe  the  cost  of  using  that  facility  is  cheaper  than  the  other  facilities  and  it  us 
best  to  maximize  the  use  of  the  lowest  cost  facility  first  before  maximizing  the  use  of 
any  other  facility.  Whatever  the  reason  for  wanting  to  use  one  facility  over  another, 
weighting  a  facility  will  accomplish  give  one  facility  preferene  over  another  facility. 
By  weighting  a  facility,  the  model  will  seek  to  utilize  the  most  heavily  weighted  fa¬ 
cility  first.  By  weighting  a  demand,  the  model  will  seek  to  utilize  the  most  heavily 
weighted  demand  first.  This  research  will  consider  only  weighted  models. 


2.2.1  Set  Covering  Location  Problem. 

The  SCLP  is  focused  on  minimizing  the  number  of  facilities  based  on  specific  “cover” 
criteria  to  meet  demands.  The  term  cover  refers  to  a  minumum  or  maximum  allowable 
distance  or  time  that  must  be  met.  In  an  SCLP,  each  demand  must  be  covered  by  a 
facility.  The  objective  in  an  SCLP  is  to  cover  all  of  the  demands  using  the  minimum 
number  of  facilities  [23]. 

Model  Development 
The  following  notation  is  used. 

Set  Notation 

•  /  =  the  set  of  demand  nodes  indexed  by  i. 

•  J  =  the  set  of  facility  locations  indexed  by  j. 

•  Ni  =  {j\dij  <  Dc}  =  the  set  of  nodes  j  within  distance  or  time  D  of  node  i\ 
these  are  the  nodes  eligible  to  house  facilities  which  “cover”  node  i. 

Parameters 

•  dij  =  the  distance  between  demand  node  i  and  location  j. 
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•  Dc  =  a  distance  or  time  standard;  a  facility  sited  at  some  node  j  within  the 
standard  of  a  demand  node  i  is  eligible  to  serve  the  demand  node. 

Decision  Variables 

{1  if  we  locate  at  site  j 
0  if  not 

With  the  given  set  notation,  decision  variables,  and  parameters  the  SCLP  forumlation 
is  shown  below. 


Minimize  Xj 

(1) 

J'6  J 

such  that  Xj  >  1, 

Vi  G  /, 

(2) 

jeNi 

Xj  e  {0,  l}, 

V j  e  J. 

(3) 

The  objective  function  (1)  minimizes  the  total  number  of  facilities  needed.  Constraint 
(2)  ensures  that  every  demand  in  the  system  is  covered  by  at  least  one  facility.  Con¬ 
straint  (3)  restricts  the  decision  variable  to  binary  values. 

One  of  the  first  set  covering  problems  involved  the  location  of  emergency  service  facil¬ 
ities.  Toregas  et  al.  [38]  asked  the  question,  “where  should  city  planners  have  located 
emergency  service  facilities  so  that  all  households  (the  demand)  had  equal  access  to 
coverage?”  The  critical  consideration  was  the  maximum  time  or  distance  a  house¬ 
hold  was  from  an  emergency  care  facility.  For  example,  a  household  could  be  located 
no  more  than  20  minutes  from  an  emergency  service  facility,  or  a  household  could 
be  located  no  more  than  15  miles  from  an  emergency  service  facility.  There  was  an 
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assumption  that  every  household  must  be  covered,  and  so  a  minimal  but  unbounded 
number  of  facilities  would  be  built  to  meet  the  demand. 

The  weighted  set  location  covering  problem  (WSCLP)  incorporates  a  weight  hj  that 
allows  the  decision  maker  to  give  preference  or  penalty  to  potential  facility  locations. 
The  notation  for  the  WSCLP  and  the  constraints  are  the  same  as  the  SCLP,  but  the 
objective  function  would  change  as  shown  below: 

Minimize  hjXj  (4) 

J'6  J 

If  we  capacitate  an  SCLP,  it  becomes  a  type  of  assignment  problem,  which  we  will 
discuss  in  Section  2.2.3. 

2.2.2  Maximal  Covering  Location  Problem  (MCLP). 

The  MCLP  is  related  to  the  SCLP.  The  differences  are  that  the  MCLP  requires  that 
all  demand  nodes  must  be  covered  and  the  number  of  facilities  to  be  located  is  fixed. 

Model  Development 

The  following  additional  notations  beyond  those  designated  in  the  SCLP  formulation 
are  described  below. 

Decision  Variables 

{1  if  demand  node  i  is  covered 
0  if  not 

Parameters 

•  hi  =  the  demand  at  node  i. 

•  p  =  the  number  of  facilities  to  locate. 
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With  the  given  set  notation,  decision  variables,  and  parameters  the  MCLP  formula¬ 
tion  is  shown  below. 


Maximize  Zi 

(5) 

iei 

such  that  V''  Xj  —  Zi  >  0, 

Vi  g  /, 

(6) 

j£Ni 

= v, 

(7) 

j&J 

Xj  e  {0,  l}, 

V j  e  J, 

(8) 

Zi  G  {0, 1}, 

Vi  g  /. 

(9) 

The  objective  function  (5)  maximizes  the  total  demand  covered.  With  an  MCLP 
there  is  no  guarantee  that  every  demand  node  is  covered.  Constraint  (6)  ensures 
that  a  demand  node  is  not  counted  as  being  covered  unless  it  is  sufficiently  close  to 
a  facility.  Constraint  (7)  fixes  the  number  of  facilities  to  be  located.  Constraints  (8) 
and  (9)  are  binary  logical  constraints  on  facility  siting  decisions  and  demand  node 
coverage,  respectively  [23]. 

Church  &  Velle  [8]  used  the  MCLP  algorithm  to  help  a  school  district  with  their  bus¬ 
ing  strategy  that  had  a  goal  to  have  a  certain  percentage  of  its  school  age  population 
within  a  certain  walking  distance  of  school.  With  more  students  within  walking  dis¬ 
tance,  the  school  district  could  save  money  on  busing  capital  and  operationals  costs. 
However,  satisfying  the  goal  required  more  schools  be  built  than  the  budget  would 
allow.  Essentially,  the  MCLP  imposed  a  budget  constraint  by  limiting  the  number  of 
schools  that  could  be  built. 


18 


The  weighted  maximal  covering  location  problem  (WMCLP)  uses  a  weight  hj  that 
allows  the  decision  maker  to  give  preference  or  penalty  to  desired  locations.  Assume 
that  a  location  equates  to  a  specific  city  and  larger  cities  are  considered  more  impor¬ 
tant.  In  the  MCLP  all  cities  are  considered  equally  important.  The  objective  function 
covers  as  many  cities  as  possible  without  regard  to  city  size.  In  the  WMCLP  larger 
cities  would  be  covered  before  smaller  cities.  The  notation  for  the  WMCLP  and  the 
constraints  are  the  same  as  the  MCLP,  but  the  objective  function  would  change  as 
shown  below: 


Maximize  htZ{  (10) 

iei 

If  we  capacitate  an  MCLP,  it  becomes  a  type  of  assignment  problem  which  we  will 
discuss  in  Section  2.23. 

2.2.3  p-Center  Problem. 

The  PCP  requires  a  fixed  number  of  facilites  and  minimizes  the  maximum  distance 
that  a  demand  is  from  its  closest  assigned  facility. 

Model  Development 

The  following  additional  notations  beyond  those  designated  in  the  SCLP  formulation 
are  described  below. 

Decision  Variables 

•  W  —  the  maximum  distance  between  a  demand  node  and  the  facility  to  which 
it  is  assigned. 

fl  if  demand  node  i  is  assigned  to  a  facility  at  node  j 
0  if  not 
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With  the  given  set  notation,  decision  variables,  and  parameters  the  MCLP  formula¬ 
tion  is  shown  below. 


Minimize  W 

(ii) 

such  that  Xj  =  p, 

(12) 

3&J 

Vij  = 

Vi  G  /, 

(13) 

ieJ 

o' 

VI 

B 

1 

Vi  G  I,j  G  J, 

(14) 

\v  y  }  dijUij  v  o, 

Vi  G  /, 

(15) 

ieJ 

xj  G  {0,1}, 

Vj  G  J, 

(16) 

yn  e  {0,1}, 

Vi  G  I,j  G  J. 

(17) 

The  objective  function  (11)  minimizes  the  maximum  distance  between  a  demand  node 
and  its  closest  open  facility.  Constraints  (12)  and  (13)  designate  how  many  facilities 
will  be  located  and  ensure  that  each  demand  node  is  assigned  to  exactly  one  facility, 
respectively.  Constraint  (14)  allows  node  assignments  to  open  facilities  only.  Con¬ 
straint  (15)  sets  a  lower  bound  for  the  maximum  distance  that  is  being  minimized. 
Constraint  (16)  restricts  the  decision  variable  to  binary  values.  Constraint  (17)  re¬ 
quires  the  demand  at  a  node  to  be  asigned  to  one  facility  only  [23]. 

The  PCP  was  developed  and  solved  by  Hakimi  [22],  The  PCP  is  similar  to  the  SCLP 
and  the  MCLP  in  that  they  all  have  covering  criteria  that  drive  the  objective  function. 
The  PCP  and  the  SCLP  both  require  that  all  demands  are  covered  by  a  facility,  but 
only  the  PCP  requires  a  fixed  number  of  facilities.  The  PCP  and  the  MCLP  both 
require  a  fixed  number  of  facilities,  but  only  the  PCP  requires  that  all  demands  are 
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covered  by  a  facility,  albeit  with  a  more  permissible  definition  of  coverage. 

Mumphrey  et  al.  [32]  advanced  Hakimi’s  work  in  this  area,  but  called  it  “Locating 
Controversial  Facilities.”  As  an  example,  consider  a  ground  distribution  trucking 
company  that  wants  to  minimize  the  longest  distance  any  of  its  trucks  travels.  The 
company  can  locate  only  p  facilities,  and  it  wants  to  place  them  such  that  the  furthest 
distance  any  delivery  truck  travels  is  as  short  as  possible.  The  PCP  can  be  thought 
of  as  an  “equity”  approach  to  facility  location  problems,  meaning  that  the  decision 
maker  desires  that  no  truck  travels  farther  than  another  truck,  and  will  locate  facili¬ 
ties  to  minimize  the  greatest  distance  any  individual  vehicle  travels. 

The  weighted  p-center  problem  (WPCP)  uses  a  weight  hij  that  allows  the  decision 
maker  to  give  preference  or  penalty  to  demands.  The  notation  for  the  WPCP  and 
the  constraints  are  the  same  as  the  SCLP,  but  constraint  (15)  would  change  as  shown 
below: 


W  -  ^2  hijdijUij  >0,  Vi  G  I.  (18) 

ieJ 

The  capacitated  p-center  problem  (CPCP)  supposes  that  for  a  given  facility  at  j , 
it  can  support  at  most  n  demand  locations.  Accordingly,  the  CPCP  would  add  the 
following  constraint: 


^2vij<n,  Vj  G  J.  (19) 

iei 
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2.2.4  p-Median  Problem. 


The  PMP  requires  a  fixed  number  of  facilites  and  minimizes  the  average  distance 
between  the  facilities  and  the  demands.  The  notation  for  the  PMP  is  identical  to 
that  introduced  for  the  PCP,  and  the  PMP  formulation  is  shown  below. 


Minimize  djji/ij, 

(20) 

iei  jeJ 

such  that  Xj  =  p, 

(21) 

16  J 

Vij  ~  1  ’ 

Vi  e  /, 

(22) 

16  J 

o' 

VI 

■c-i 

B 

1 

Vi  e  i,j  e  J, 

(23) 

Xj  e  {0,1}, 

Vj  e  J, 

(24) 

yn  e  (0,  l}, 

Vi  g  IV j  e  J. 

(25) 

The  objective  function  (20)  minimizes  the  average  distance  between  demands  and 
assigned  facilities.  Constraints  (21)  and  (22)  designate  how  many  facilities  will  be 
located  and  ensure  that  each  demand  node  is  assigned  to  exactly  one  facility,  respec¬ 
tively.  Constraint  (23)  allows  node  assignments  only  to  open  facilities.  Constraint 
(24)  restricts  the  decision  variable  to  binary  values.  Constraint  (25)  requires  the  de¬ 
mand  at  a  node  be  assigned  to  only  one  facility  [23] . 

Current  et  al.  [11]  provide  an  application  of  a  PMP  using  20  demand  nodes  labeled 
A,  B,...,T  with  10  of  those  nodes  considered  possible  intermediate  locations,  namely 
A,  B,...,J,  as  shown  in  Figure  7  [11].  In  this  problem  one  of  the  constraints  called  for 
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exactly  two  facilities  to  be  built  such  that  when  all  demands  were  routed  to  one  of 
the  two  facilities,  the  sum  of  the  distances  was  minimized. 

The  optimal  solution  is  to  place  nodes  at  B  and  J  to  minimize  the  total  length  of 
the  arcs.  The  solution  can  be  understood  graphically  by  examining  Figure  7.  The 
optimal  solution  was  5133  units  when  nodes  B  and  J  were  selected.  [Note:  [11]  did 
not  provide  units  to  his  research  so  the  generic  term  “units”  is  used.]  The  pairing 
with  the  closest  total  distance  traveled  to  nodes  B  and  J  was  B  and  H  at  a  distance 
of  6151  units.  Placing  the  facilities  at  nodes  F  and  J  resulted  in  a  distance  of  116,410 
units  which  was  the  maximum  distance  for  this  network. 


Optimal  EOL  Solution:  P=2 

Opportunity  Loss  =  5133 


|~ □  Node  Name  Assignments  to  a  Open  Facility 

Figure  7.  p-Median  Example 

The  weighted  p-median  problem  (WPMP)  uses  a  weight  hl3  that  allows  the  decision 
maker  to  give  preference  or  penalty  to  desired  demands.  The  notation  for  the  WPMP 
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and  the  objective  function  would  change  as  shown  below: 


Minimize  EE  hidijUij  (26) 

iei  j&J 

The  capacitated  p-median  problem  (CPMP)  supposes  that,  for  a  given  facility  at  j, 
it  can  support  at  most  n  demand  locations.  The  CPMP  would  add  the  following 
constraint: 


Vj  G  J.  (27) 

iei 

2.3  Network  Flow  Models 

In  1954  T.E.  Harris  and  F.S.  Ross  formulated  one  of  the  first  known  flow  problems. 
In  1956  Dantzig,  G.B,,  L.  R.  Ford,  and  D.R.  Fulkerson  furthered  this  work  by  de¬ 
veloping  the  first  algorithm  to  compute  the  maximum  flow  through  a  network  [13]. 
This  initial  work  opened  exciting  opportunities  in  the  discipline  of  network  modeling 
and  was  followed  by  a  more  than  a  decade  of  advancements  in  this  held.  Three  of 
the  more  notable  advancements  were  by  Bellman  [5]  with  the  SPP,  Dijkstra  [17]  with 
another  SPP  algorithm,  and  Dantzig’s  work  on  minimum  cost  how  models  [12], 

Network  how  models  are  used  to  find  the  most  efficient  path  for  moving  a  commodity 
from  a  source  node  to  a  sink  node  through  a  group  of  interconnected  transshipment 
nodes.  The  goal  is  to  move  these  goods  across  the  network  to  satisfy  the  most  demand 
at  the  cheapest  cost,  the  shortest  path,  or  the  maximum  how.  A  network  has  at  least 
one  source  node  (supply  node)  and  at  least  one  sink  node  (demand  node).  Goods 
how  through  the  network  from  the  source  node  to  the  sink  node.  At  the  source  node 
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the  amount  of  flow  out  of  the  node  is  always  greater  than  the  amount  of  flow  into 
the  node.  The  opposite  is  true  at  the  sink  node.  Any  other  node  that  may  exist 
in  the  network  is  a  transshipment  node.  The  amount  of  flow  into  a  transshipment 
node  is  equal  to  the  amount  of  flow  out  of  a  transshipment  node  [24],  Networks  exist 
in  the  fields  of  transportation,  communications,  manufacturing,  power  distribution, 
resource  management,  financial  planning,  and  many  others. 

The  study  of  network  flow  models  in  effect  addresses  three  questions  [1]. 

1.  “What  is  the  best  way  to  traverse  a  network  to  get  from  one  point  to  another 
as  cheaply  as  possible?” 

2.  “If  a  network  has  capacities  on  arc  flows,  how  can  we  send  as  much  flow  as  pos¬ 
sible  between  two  points  in  the  network  while  honoring  the  arc  flow  capacities?” 

3.  “If  we  incur  a  cost  per  unit  flow  on  a  network  with  arc  capacities  and  we  need  to 
send  units  of  a  good  that  reside  at  one  or  more  points  in  the  network  to  one  or 
more  other  points,  how  can  we  send  the  material  at  a  minimum  possible  cost?” 

These  questions  motivate  three  of  the  most  common  types  of  network  flow  models:  the 
shortest  path  problem  (SPP),  the  maximum  flow  problem  (MFP),  and  the  minimum 
cost  flow  problem  (MCFP). 

2.3.1  Shortest  Path  Problem. 

In  this  section,  we  present  and  discuss  the  formulation  for  the  SPP.  Although  several 
other  versions  of  the  SPP  exist,  this  model  is  a  simple  network  with  one  source  and 
one  sink.  More  complex  models  with  multiple  sources  and  multiple  sinks  are  not 
covered. 

Model  Development 
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The  following  notation  is  used. 

Set  Notation 

•  A  =  the  set  of  m  directed  arcs. 

•  N  =  the  set  of  n  nodes. 

•  G  =  (N,  A)  the  underlying  network. 

•  /  =  the  set  of  nodes  indexed  by  i. 

•  J  =  the  set  of  candidate  facility  loctaions  indexed  by  j. 

Decision  Variables 

•  Xij  =  a  binary  decision  variable  that  reflects  whether  arc  (i,j)  is  identified  as 
being  on  the  shortest  path  from  s  to  t. 

Parameters 

•  Cjj  =  the  arc  length  or  arc  cost  from  node  %  to  node  j. 

•  s  —  the  source  node. 

•  t  =  the  sink  node. 

With  the  given  set  notation,  decision  variables,  and  parameters  the  SPP  formulation 
is  shown  below. 


Minimize  E  Cij  Xij , 


(28) 


such  that  5]  xa~ 

j:(i,j)eA  j-UA&A 


1  for  i  =  s 
(  0  for  all  i  E  N  \  {s,t} 
—  1  for  i  —  t 


(29) 
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=  {o.l}.  V(i,j)eA 


(30) 


The  objective  function  (28)  minimizes  the  summation  of  arc  flows  x%1  times  their  costs 
Cij.  The  parameter  cl3  can  represent  either  a  distance  or  a  monetary  value.  Constraint 
(29)  is  commonly  referred  to  as  flow  conservation  or  nodal  balance  equations.  This 
constraint  require  that  one  unit  of  flow  begins  at  s  and  ends  at  t.  Constraint  (30)  is 
a  binary  variable  that  requires  Xij  to  be  a  1  or  0. 

The  most  basic  application  of  the  SPP  is  using  a  road  map  to  move  from  a  start 
location  to  an  end  location.  Although  many  s  —  t  paths  may  exist,  there  is  only  one 
shortest  length  for  an  s  —  t  path  (such  a  shortest  length  path  may  correspond  to  more 
than  one  s  —  t  path,  as  in  the  case  of  alternative  optima).  Beasley  &  Christofides  [4] 
used  an  integer  programming  formulation  to  solve  a  “resource-constrained  SPP.”  A 
traveler  was  given  a  budget  of  various  resources  and  was  required  to  reach  a  given 
destination  as  quickly  as  possible.  The  various  resources  were  consumed  along  the 
path.  He  had  to  reach  his  destination  within  the  resource  constraints  of  his  budget. 
Beasley  &  Christofides  [4]  used  a  lagrangean  relaxation  to  find  the  shortest  path  and 
satisfy  the  budget  constraints. 

2.3.2  Maximum  Flow  Problem. 

In  this  section,  we  present  and  discuss  the  formulation  for  the  MFP.  The  MFP  sends 
as  much  flow  as  possible  through  a  network  between  a  source  node  s  and  a  sink  node 
t  without  exceeding  the  capacity  of  the  arcs. 

Model  Development 
The  following  notation  is  used. 

Set  Notation  As  defined  for  the  SPP. 
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Decision  Variables 


•  v  =  a  scalar  variable  representing  the  value  of  the  maximum  s  —  t  flow. 

•  Xij  =  the  flow  from  i  to  j. 

Parameters 

•  s  and  t  remain  as  defined  in  the  SPP. 

•  Uij  =  the  non- negativity  flow  capacity  for  arc  (i,j). 

With  the  given  set  notation,  decision  variables,  and  parameters  the  MFP  formulation 
is  shown  below. 


Maximize  v 


(31) 


such  that 


%ji  \ 


v  for  i  —  s 
0  for  all  i  £  N  \  {s,  t} 
—v  for  i  —  t 


(32) 


0  <  xi:j  <  u^,  V(i,j)eA 


(33) 


The  objective  function  (31)  maximizes  the  amount  of  flow  v  through  the  network.  At 
the  start  node  s,  constraint  (32)  ensures  that  the  amount  of  flow  out  of  the  node  is 
greater  than  the  amount  of  flow  into  the  node.  Constraint  (32)  also  ensures  that  at 
the  end  node  the  amount  of  flow  in  to  the  node  is  greater  than  the  amount  of  flow 
out  of  the  node.  This  constraint  also  ensures  that  at  every  transshipment  node  in 
the  network  the  flow  into  the  node  is  equal  to  the  amount  of  flow  out  of  the  node. 
Constraint  (33)  ensures  that  the  flow  xtJ  will  always  be  non-negative. 


The  MFP  is  used  to  model  the  flow  of  petroleum  products  through  a  pipe  network, 
cars  in  a  road  network,  messages  in  a  telecommunication  network,  and  electricity  in 
an  electrical  network.  The  MFP  models  arc  capacities  but  not  arc  costs,  whereas  the 
SPP  models  arc  costs  but  not  arc  capacities;  they  are  complementary.  Together,  the 
properties  of  the  MFP  and  SPP  make  np  the  basis  for  network  flow  analysis  and  lead 
to  the  minimum  cost  flow  model  described  in  the  next  section  [1], 

The  Ford-Fulkerson  method  (1956)  was  the  first  algorithm  recorded  for  computing 
maximum  flow.  In  1970  Jack  Edmonds  and  Richard  M.  Karp  developed  a  more  ef¬ 
ficient  maximum  flow  algorithm.  The  Edmonds-Karp  algorithm  is  identical  to  the 
Ford-Fulkerson  algorithm  with  one  exception.  Edmonds  and  Karp  pointed  out  that 
“an  improper  choice  of  flow  augmenting  paths  can  lead  to  severe  computational  dif¬ 
ficulties.”  Edmonds  and  Karp  showed  that  “if  each  flow  augmentation  is  made  along 
an  augmenting  path  having  a  minimum  number  of  arcs,  then  a  maximum  flow  in  an 
n-node  network”  can  be  obtained  more  quickly  [18]. 

2.3.3  Minimum  Cost  Flow  Problem. 

The  SPP  is  a  special  case  of  the  MCFP,  and  the  MCFP  is  closely  related.  The  MCFP 
seeks  to  “determine  a  least  cost  shipment  of  a  commodity  through  a  network  in  order 
to  satisfy  demands  at  certain  nodes  from  available  supplies  at  other  nodes”  [1]. 

Model  Development 

The  following  notation  is  used. 

Set  Notation  As  defined  in  the  SPP. 

Decision  Variables 

•  Xij  =  the  flow  from  i  to  j. 

Parameters 
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•  b,  —  a  variable  indicating  the  amount  of  supply  at  node  i,  where  6,;  >  0  indi¬ 
cates  a  supply  node,  bt  <  0  indicates  a  demand  node,  and  6*  =  0  indicates  a 

transshipment  node. 

•  Cij  =  the  cost  of  a  unit  of  flow  from  i  to  j. 

•  Uij  =  the  upper  bound  on  the  flow  from  i  to  j. 

With  the  given  set  notation,  decision  variables,  and  parameters  the  MFP  formulation 
is  shown  below. 


Minimize  CijXij 

(34) 

h,i)eA 

such  that  Xij  —  xi3  = 

Vi  G  N, 

(35) 

j:(i,j)eA  j'(j,i)£A 

0  ^  Xij  ^  Uij , 

V(i,  j)  e  A. 

(36) 

The  objective  function  (34)  minimizes  the  cost  of  the  flow  through  the  network.  Con¬ 
straint  (35)  ensures  a  balance  of  flow  at  each  node.  Constraint  (36)  ensures  that  the 
flow  Xij  will  always  be  non-negative  and  bounded. 

The  MCFP  is  used  to  model  product  distribution  from  the  location  of  product  cre¬ 
ation  to  a  storage  facility  such  as  a  warehouse;  distribution  from  warehouses  to  re¬ 
tailers;  the  flow  of  raw  materials  used  in  an  industrial  process  as  they  move  through 
various  stations  in  a  production  line;  and  the  flow  of  cars  through  a  network  of  city 
streets.  The  MCFP  algorithms  are  not  as  efficient  as  the  SPP  and  MFP  algorithms, 
but  the  MCFP  algorithms  are  more  versatile  because  they  address  both  flow  costs 
and  capacities.  Fulkerson  &  Harding  [19]  used  a  MCFP  algorithm  to  maximize  the 
minimum  source-sink  path  subject  to  a  budget  constraint.  With  a  fixed  budget  for 
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each  expenditure  on  arcs  of  the  network,  they  allocated  the  budget  among  the  arcs 
so  as  to  maximize  the  length  of  the  shortest  path  from  source  to  sink”  [19]. 

Of  note,  there  are  three  special  cases  of  the  MCFP  (other  than  the  SPP)  commonly 
found  in  the  literature:  the  general  transportation  problem,  the  assignment  problem, 
and  the  transshipment  problem. 

The  transportation  problem  is  identical  to  the  MCFP  with  the  following  exceptions. 
First,  the  set  of  nodes  N  is  divided  into  two  sets  of  nodes  (possibly  unequal),  Ni  and 
N2,  n  1  is  a  set  of  supply  nodes  and  N2  is  a  set  of  demand  nodes.  No  transshipment 
nodes  are  allowed.  Hillier  &  Lieberman  [24]  states  that  “all  arcs  are  directed  from  a 
supply  node  to  a  demand  node,  where  distributing  xtj  units  from  source  i  to  destina¬ 
tion  j  corresponds  to  a  flow  of  xtj  through  arc  i  — y  j” .  The  transportation  problem 
does  not  impose  an  upper  bound  on  xVJ ,  therefore  Uij  =  00.  The  distribution  of  goods 
from  warehouses  to  customers  is  an  example  of  a  transportation  problem  where  N\ 
represents  the  warehouses,  N2  represents  the  customers,  and  an  arc  (i,j)  represents 
a  path  from  warehouse  i  to  customer  j  [1] . 

Like  the  transportation  problem,  the  assignment  problem  is  identical  to  the  MCFP 
with  some  exceptions.  The  assignment  problem  divides  the  set  of  nodes  N  into  two 
sets  of  nodes,  N\  and  N2,  which  represent  supply  and  demand  nodes,  respectively. 
However,  the  cardinality  of  the  two  new  sets  must  be  equal.  Also  each  supply  node 
has  non-negative  flow  and  each  demand  node  has  a  negative  flow.  The  assignment 
problem  seeks  to  create  at  a  minimum  cost  a  collection  of  pairs,  one  from  each  set. 
The  pairings  must  be  one-to-one.  The  most  common  examples  of  assignment  prob¬ 
lems  are  assigning  jobs  to  machines,  students  to  a  computer,  runners  to  a  lane  at  a 
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track  meet,  and  aircraft  to  a  runway  [24]  and  [1]  . 


The  last  special  case  of  the  MCFP  is  the  transshipment  problem.  One  difference  be¬ 
tween  the  transshipment  problem  and  the  MCFP  is  that  the  transshipment  problem 
does  not  have  finite  arc  capacities.  For  the  transshipment  problem,  an  upper  bound 
on  the  amount  of  goods  that  can  flow  across  an  arc  is  limitless;  each  arc  can  carry  any 
desired  amount  of  flow  required  to  meet  demands  [24],  Like  both  the  transportation 
and  assignment  problems,  the  transshipment  problem  has  two  sets  of  nodes,  N\  and 
N2,  which  represent  supply  and  demand  nodes,  respectively.  Any  node  that  is  not 
within  N\  or  N2  is  called  a  transshipment  node.  A  transshipment  node  is  an  inter¬ 
mediate  node  used  as  a  temporary  stopping  location  for  goods  flowing  from  a  source 
node  to  a  sink  node.  Goods  may  flow  to  a  transshipment  node  instead  of  directly 
from  a  source  node  to  a  sink  node  because  it  is  less  costly  or  out  of  necessity.  For 
example,  consider  a  set  of  goods  leaving  a  source  node  bound  for  three  locations  via 
train.  The  first  location  is  on  an  island.  The  second  location  is  in  a  remote  region 
that  does  not  have  a  rail  system.  The  third  location  has  a  rail  system,  but  needs 
only  a  relatively  small  amount  of  goods.  To  move  goods  to  the  first  location,  either 
a  boat  or  an  airplane  is  required.  To  move  goods  to  the  second  location  an  airplane 
could  be  used,  but  it  is  most  likely  cheaper  to  use  ground  transportation.  Goods  to 
the  third  location  could  be  delivered  via  train,  but  it  may  be  cheaper  to  use  ground 
transportation  due  to  the  relatively  small  demand.  To  meet  the  demands  of  all  three 
locations,  goods  should  flow  into  a  transshipment  node  first  and  then  switch  to  the 
most  economical  mode  of  transportation. 
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2.4  Location  and  Network  Flow  Problem 


The  models  previously  described  are  well  defined  and  can  be  formulated  to  fit  many 
situations.  However,  there  are  scenarios  where  both  facility  location  and  flow  must 
be  considered  in  the  same  network.  Consider  a  solid  waste  collection  transfer  where 
trash  is  collected  and  then  distributed.  The  flow  into  the  substation  from  the  various 
pick-up  routes  will  have  a  maximum  capacity  and  so  will  the  transfer  station  [1], 
In  these  situations  using  the  facility  location  problem  models  and  network  location 
models  to  derive  a  solution  is  insufficient.  Specifically,  representing  the  actual  cost  of 
this  operation  by  modeling  a  single  round  trip  from  a  facility  to  a  customer  would  be 
inadequate  [42], 

Location  and  network  flow  problems  are  also  called  location  routing  problems  (LRP). 
This  name  is  a  misnomer  because  the  LRP  is  not  a  well  defined  problem  with  distinct 
parameters.  An  LRP  seeks  to  solve  both  facility  location  problems  and  network  flow 
problems.  The  LRP  accounts  for  the  individual  parameters  of  the  facility  location 
problem,  the  individual  parameters  of  the  network  location  problem  and  the  relation¬ 
ship  between  the  two.  The  LRP  is  also  able  to  consider  multiple  stops  throughout 
the  network  [33]. 

The  most  pertinent  field  of  study  to  this  research  is  the  LRP.  In  2014,  Prodhon  and 
Prins  [36]  provided  a  thorough  review  of  LRP  advancements.  Much  of  the  earlier 
research  endeavors  such  as  the  warehouse  location- routing  problem  (see,  e.g.,  [2,  35]) 
are  limited  in  scope  and  not  immediately  relevant  to  our  problem;  they  only  address 
the  location  of  one  echelon  of  sustainment  and  the  routing  of  supplies  to  meet  demand. 
More  closely  related  is  research  on  two-echelon  location-and-routing  problems,  as 
initially  examined  by  Jacobsen  and  Madsen  [27]  and  Madsen  [30].  A  general  form 
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of  their  two-echelon  location-routing  problem  (denoted  ‘LRP-2E’)  has  been  studied 
more  extensively  in  recent  years  (e.g.,  see  [6,  10,  29,  34]).  Although  the  LRP-2E 
model  is  the  closest  model  in  the  literature  to  our  application,  it  still  differs  in  several 
key  aspects.  The  LRP-2E  model  seeks  to  locate  intermediate  distribution  centers 
wherein  only  the  highest  and  lowest  echelon  locations  are  fixed,  whereas  we  also  seek 
to  locate  the  highest  level  echelon  (i.e.,  the  A/SPODs).  Moreover,  the  LRP-2E  model 
routes  the  delivery  of  supplies  via  tours.  This  constraint  is  appropriate  for  certain 
applications  such  as  newspaper  or  grocery  delivery,  but  the  delivery  of  supplies  on  a 
road  network  for  our  context  is  not  similarly  restricted. 
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III.  Methodology 


“Clearly,  logistics  is  the  hard  part  of  fighting  a  war.” 

Lt.  Gen.  E.  T.  Cook,  USMC  [25] 

This  chapter  explains  the  methodology  that  we  use  to  solve  this  problem.  First,  in 
Section  3.1  we  recall  the  need  for  a  mathematical  model  to  provide  support  to  TRAC- 
LEE’s  logistics  studies.  Next,  in  Section  3.2  we  set  forth  a  formal  presentation  of 
the  math  programming  model.  Finally,  Section  3.3  discusses  the  software  used  to 
implement  the  model  and  solve  the  problem. 

3.1  Overview 

Recall  that  TRAC-LEE  uses  LBC  (a  dynamic  stochastic  simulation)  to  assess  the 
effectiveness  of  locations  for  A/SPODs  and  intermediate  logistics  hubs.  LBC  requires 
initial  solutions  which  currently  are  provided  by  subject  matter  experts  (SMEs).  This 
research  seeks  to  provide  TRAC-LEE  with  analytic  starting  solutions  for  logistics  hub 
locations  (APODs,  SPODs,  and  DCs). 

United  States  Army  Africa  (USARAF)  and  TRAC-LEE  are  collaborating  to  develop 
supply  chain  models  for  contingency  operations  within  the  United  States  Africa  Com¬ 
mand  (AFRICOM)  theater  of  operations.  TRAC-LEE  uses  LBC  to  provide  USARAF 
with  supply  chain  analysis.  This  research  is  conducted  to  assist  TRAC-LEE  in  pro¬ 
viding  improved  support  to  USARAF.  The  model  will  use  Nigeria  as  a  sample  area 
of  operations  (AO).  Once  the  model  is  validated,  it  can  be  applied  to  any  country 
within  AFRICOM. 
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3.2  Model  Development 


This  research  models  the  movement  of  supplies  through  a  transportation  network  as 
a  mixed-integer  linear  program  (MILP).  Within  the  network,  various  units  transport 
the  military  classes  of  supply  from  CONUS  through  A/SPODs  into  the  AO,  over 
air  and  ground  lines  of  communication  (A/GLOCs)  through  logistic  hubs  managed 
by  regional  warehouse  distribution  centers  (DCs),  and  subsequently  over  GLOCs  to 
units  at  forward  operating  bases  (FOBs).  Our  model  seeks  to  locate  a  fixed  number  of 
SPODs,  APODs,  and  DCs  in  order  to  meet  the  weekly  demands  for  supplies  at  FOBs 
and  logistic  hubs  while  minimizing  one  of  several  considered  objectives.  As  a  baseline, 
we  consider  the  objective  of  minimizing  the  total  risk  of  transporting  supplies,  but  we 
also  consider  minimizing  the  total  transportation  distances  for  the  supplies  as  well  as 
minimizing  the  imbalance  of  per  unit  capita  responsibility  for  supply  transportation 
distances. 

We  represent  the  A/SPOD-to-FOB  distribution  infrastructure  as  a  directed  network, 
wherein  seaports,  airports,  FOBs,  and  road  junctions  correspond  to  nodes,  and  roads 
and  intra-theater  air  transportation  routes  correspond  to  directed  arcs.  Over  this  net¬ 
work,  we  aggregate  the  point-demand  for  supplies  on  a  weekly  basis,  and  we  model 
the  flow  of  supplies  specific  to  one  of  two  echelon  levels:  between  the  A/SPOD  and 
the  DCs,  and  between  the  DCs  and  the  FOBs. 

Assumptions 

Underlying  our  model  are  five  key  assumptions.  First,  we  will  consider  locating 
a  facility  only  at  a  node  within  the  network.  This  means  that  a  facility  cannot  be 
placed  along  an  arc  or  off  of  the  network  itself.  This  assumption  limits  the  complexity 
of  the  formulation  and  is  valid  for  most  objective  functions  considered  [14],  Second, 
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a  DC  node  can  be  located  at  any  node  in  the  network.  This  assumption  ensures  the 
widest  latitude  for  the  decision  maker  for  emplacing  DCs  within  the  limits  of  the  first 
assumption.  Third,  all  arcs  are  directed.  This  ensures  that  we  can  characterize  the 
flow  over  an  arc  by  its  direction  of  movement.  Fourth,  there  is  an  unbounded  supply  of 
material  from  CONUS.  This  assumption  is  used  because  we  are  concerned  both  with 
the  ability  of  in-theater  infrastructure  and  logistics  nodes  to  meet  demands.  Finally, 
we  will  ship  only  the  amount  of  supplies  needed  to  meet  the  demands.  This  implies 
that  we  will  not  ‘warehouse’  supplies.  This  research  models  the  sustained  conditions 
after  initial  deployment  into  the  area  of  operations  when  steady-state  conditions  have 
been  achieved.  By  limiting  flow  to  exactly  meet  demands,  we  can  obtain  a  more 
accurate  assessment  of  how  supplies  flow  from  an  A/SPOD  to  a  FOB  instead  of  from 
a  warehouse  in  the  network  to  a  FOB. 

Set  Notation 

•  N  :  the  set  of  nodes  in  the  network  indexed  by  i. 

—  Na  C  N  :  the  subset  of  nodes  in  the  network  that  can  serve  as  an  APOD. 

—  Ns  C  N  :  the  subset  of  nodes  in  the  network  that  can  serve  as  an  SPOD. 

•  A  :  the  set  of  arcs  in  the  network  indexed  by 

•  G[N,A]  :  the  network. 

•  E  :  the  echelon  of  commodity  flow  (indexed  by  e)  in  the  network,  with  E  = 
{1,2}  corresponding  to  flow  between  the  A/SPOD  and  intermediate  distribution 
centers,  and  to  flow  between  the  intermediate  distribution  centers  and  the  BSBs, 
respectively. 

Parameters 
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•  b,  :  the  amount  of  weekly  demand  for  supplies  at  node  i  E  N,  based  on  the 
geographic  locations  of  BCTs. 

•  bdc  :  the  generic  weekly  demand  for  an  intermediate  distribution  center  (DC). 

•  cfj  :  the  per  unit  risk  incurred  to  transport  supplies  of  echelon  level  e  on  arc 

(i,j)  e  A. 

•  dij  :  the  distance  from  node  i  to  node  j,  defined  for  (i,j)  €  A. 

•  na  :  the  maximum  number  of  APODs  being  designated  for  CONUS-into-theater 
flow. 

•  ns  :  the  maximum  number  of  SPODs  being  designated  for  CONUS-into-theater 
flow. 

•  ne  :  the  number  of  medium  truck  companies  assigned  to  transport  supplies  at 
echelon  e  G  E. 

•  ride  '■  the  maximum  number  of  intermediate  distribution  centers  (DCs)  being 
located. 

•  :  weekly  upper  bound  on  throughput  for  an  APOD  located  at  node  i  G  Na. 
This  allows  the  model  to  account  for  limitations  on  maximum  on  ground  (MOG) 
capacity  of  potential  APODs.  MOGs  data  can  be  determined  by  Air  Mobility 
Command  (AMC)  surveys  of  potential  APODs.  For  the  purposes  of  keeping  this 
work  unclassified  U-1  :  values  were  developed  by  subject  matter  expert  input. 

•  Uf  :  weekly  upper  bound  on  throughput  for  an  SPOD  located  at  node  i  € 
Ns.  This  allows  the  model  to  account  for  limitations  on  berthing  capacities  of 
potential  SPODs.  Berthing  capacities  can  be  determined  by  (SDDC)  surveys  of 
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potential  SPODs  for  berthing  capacities.  For  the  purposes  of  keeping  this  work 
unclassified  Uf  :  values  were  developed  by  subject  matter  expert  input. 

•  Ufc  :  weekly  generic  upper  bound  on  throughput  for  a  DC  located  at  node 
i  G  N.  This  allows  the  model  to  account  for  limitations  on  repackaging  and 
managing  supplies  at  intermediate  distribution  centers  based  on  the  availability 
of  material  handling  equipment  (MHE). 

•  Uij  :  upper  bound  on  weekly  flow  on  arc  (i,j)  G  A. 

•  w i,w2,w3  :  the  non-negative  weights  on  the  objective  function  components  qi, 
Q2-,  and  qs,  respectively. 

Decision  Variables 

1.  Location  Related 

•  yf  :  1  if  an  APOD  is  located  at  node  i  G  Na,  and  0  otherwise. 

•  y\  :  1  if  an  SPOD  is  located  at  node  i  G  Ns,  and  0  otherwise. 

•  y^c  :  1  if  a  DC  is  located  at  node  i  G  N,  and  0  otherwise. 

2.  Flow  Related 

•  x?j  :  the  weekly  flow  of  supplies  at  echelon  e  G  E  on  arc  (i,j)  G  A. 

•  zf  :  the  amount  of  weekly  flow  passing  through  node  %  G  N  that  is  converted 
from  Echelon  1  to  Echelon  2.  This  ‘conversion’  represents  the  receipt, 
management,  and  repackaging  of  bulk  supplies  into  tailored  packages  for 
delivery  to  BSBs. 

•  zf  :  the  amount  of  weekly  flow  passing  through  node  i  G  JV  that  is  not 
converted  from  Echelon  1  to  Echelon  2. 
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3.  Intermediate  Decision  Variables 


•  qi  :  the  total  risk  accumulated  to  supplies  flowing  through  the  network. 

•  q-2  :  the  maximum  transportation  distance  through  the  network. 

•  g3  :  the  maximum  per  capita  (i.e.,  per  medium  truck  company)  workload 
supported  by  transportation  assets  at  a  given  echelon  e  E  E. 

Problem  Formulation 

Within  this  framework,  we  propose  the  following  formulation  for  the  Military  Theater 
Distribution  Network  Design  Problem  (MTDNDP): 


3 

Minimize  E  wkqk, 

k= 1 

Subject  to  q\  =  E  E  ^ijXij  i 

( i,j)eA  eeE 
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e  (i,J)eA 

Y  yt  -  Us’ 

i£Na 
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i£Na 

Y  yfc  ^  ndc} 
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yt  +  y‘  +  yf  <  i 
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(37) 

(38) 
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V  e  G  E,  (40) 
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(43) 

V  i  e  N,  (44) 

V  i  e  N,  (45) 
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V  4<z? 
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(48) 

y  y<z: 

+  W, 

V  i  e  Na\{Nan Ns}, 

(49) 
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\/ie{NanNs}, 

(50) 
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V  i  e  N, 

(51) 

Xtj  <  Uij, 

e£E 

V  (i,j)  e  A, 

(52) 

x]j  >  0, 

V(i,i)  e  A, 

(53) 

<2: 

IV 

o 

V(i,i)  e  A, 

(54) 

24s  e  {0,1}, 

Vt  G  Ns, 

(55) 

Vi  e  {0,1}, 

V*  e  Na, 

(56) 

vT  e  {0, 1}, 

Vi  e  N. 

(57) 

The  objective  function  (37)  minimizes  a  non-negatively  weighted  sum  of  three  mea¬ 
sures  of  performance:  the  total  risk  encountered  by  transporting  supplies,  the  total 
distance  traveled  by  supplies,  and  the  maximum  per  capita  workload  supported  by 
transportation  assets  at  a  given  echelon.  Constraints  (38)  and  (39)  calculate  the  total 
accumulated  risk  and  traveled  distance  by  units  of  supply  flowing  through  the  net¬ 
work,  respectively.  Constraint  (40)  bounds  the  average  distance  traveled  per  truck 
company  transporting  supplies  for  each  echelon  in  the  network.  Constraints  (41),  (42), 
and  (43)  enforce  the  respective  upper  bounds  on  the  number  of  SPODs,  APODs,  and 
DCs  that  may  be  emplaced,  and  constraint  (44)  restricts  location  decisions  to  prevent 
a  co-location  of  an  SPOD,  APOD,  and/or  DC  at  the  same  node.  Although  such  a 
disposition  of  sustainment  nodes  may  be  possible,  it  is  trivial  within  our  context,  and 
so  we  utilize  this  constraint  to  set  aside  the  possibility.  Should  a  reader  wish  to  allow 


41 


such  a  co-location,  one  need  only  eliminate  the  appropriate  binary  variables  from 
constraint  (44),  or  eliminated  the  constraint  in  entirely.  Constraints  (45)-(51)  en¬ 
force  conservation  of  flow  at  each  node  while  accounting  for  the  transition  of  supplies 
between  echelons  at  distribution  centers.  Constraint  (45)  requires  that  all  Echelon  1 
flow  coming  into  a  node  remains  as  Echelon  1  flow  or  is  converted  into  Echelon  2  flow. 
Constraint  (46)  ensures  that  flow  may  only  be  converted  from  Echelon  1  to  Echelon 
2  at  a  node  if  an  intermediate  distrnbntion  center  node  is  emplaced  there  and,  for 
such  a  case,  it  bounds  the  corresponding  weekly  converted  flow  by  Udc.  For  nodes  at 
which  neither  an  APOD  nor  an  SPOD  may  be  located,  constraint  (47)  requires  that 
the  amount  of  Echelon  1  flow  into  this  node  that  was  not  converted  into  Echelon  2 
flow  will  depart  the  node  as  Echelon  1  flow.  For  nodes  at  which  an  SPOD  but  not 
an  APOD  may  be  located,  constraint  (48)  ensures  that  the  amount  of  Echelon  1  flow 
departing  a  node  is  bounded  above  by  the  sum  of  the  non-converted  Echelon  1  flow 
that  entered  the  node  and  the  flow  that  entered  from  CONUS  if  that  node  is  desig¬ 
nated  as  SPOD,  where  Uf  represents  the  limit  on  weekly  throughput  into  the  port. 
Using  a  parallel  construct,  constraint  (49)  enforces  a  similar  constraint  for  nodes  that 
may  be  selected  to  serve  as  an  APOD  but  not  an  SPOD,  and  constraint  (50)  does 
the  same  for  nodes  that  may  be  selected  to  serve  as  either  an  SPOD  or  an  APOD. 
Given  the  restriction  imposed  by  constraint  (44),  it  may  appear  reasonable  to  replace 
constraints  (48)-(49)  with  constraint  (50)  indexed  over  i  €  {lVa  U  Ns };  however,  we 
maintain  them  within  this  model  for  the  purpose  of  extensibility,  in  the  event  that  a 
user  does  not  assume  that  distribution  centers  cannot  be  co-located  with  A/SPODs. 
For  every  node,  constraint  (51)  requires  that  the  sum  of  Echelon  2  flow  and  converted 
Echelon  1  flow  coming  into  a  node  is  equal  to  the  sum  of  the  Echelon  2  flow  coming 
out  of  the  node,  the  demand  at  that  node,  and  any  demand  that  is  induced  by  lo¬ 
cating  a  distribution  center  at  that  node.  Constraint  (52)  imposes  an  upper  limit  on 
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the  weekly  total  of  Echelon  1  and  Echelon  2  flow  for  each  arc  in  the  network,  con¬ 
straints  (53)  and  (54)  ensure  that  both  Echelon  1  and  2  flows  are  non-negative,  and 
constraints  (55)-(57)  enforce  binary  logical  constraints  on  facility  location  decisions. 


3.3  Software  Used  in  Modeling 

Two  types  of  software  are  used  in  this  research:  Microsoft  Excel  (Excel)  and  IBM 
ILOG  CPLEX  Optimization  Studio  (CPLEX).  Excel  is  a  spreadsheet  application  that 
allows  calculation,  graphing,  and  a  macro  programming  language  called  Visual  Basic 
for  Applications  (VBA).  CPLEX  is  a  commercial  solver  for  linear,  integer,  and  mixed 
integer  linear  programming  problems. 

We  utilize  Microsoft  Excel  to  store  the  sets  and  parameters  underlying  a  given  in¬ 
stance,  and  we  use  the  embedded  VBA  capability  to  represent  our  MILP  for  the 
instance.  Also  using  VBA,  we  invoke  the  commercial  solver  CPLEX  to  solve  our 
MILP  and  present  our  solution  to  the  decision  maker. 
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IV.  Analysis  and  Results 


“The  line  between  disorder  and  order  lies  in  logistics...” 

Sun  Tzu  [25] 

This  chapter  explains  the  results  yielded  in  this  analysis.  First,  in  Section  4.1  we 
present  the  scenario  which  was  used  to  test  the  modeling  tool.  Next,  Section  4.2 
discusses  the  results  using  seven  weighting  variations.  Finally,  in  Section  4.3  we 
include  a  sensitivity  analysis  using  one  of  the  three  components  of  the  objective 
function. 

4.1  Scenario  Description 

To  test  our  model,  we  developed  an  unclassified  representative  humanitarian  assis¬ 
tance  (HA)  scenario  in  the  country  of  Nigeria  described  in  Appendix  A.  Although 
any  country  could  be  used,  we  chose  Nigeria  because  it  provides  the  opportunity  for 
selection  among  multiple  A/SPODs,  and  there  was  sufficient  available  data  to  keep 
this  research  unclassified. 

USARAF  is  responsible  for  strategic,  operational,  and  tactical  operations  across  the 
full  range  of  military  operations  for  the  African  continent.  Moving  military  supplies 
into  an  AO  to  meet  war  fighter  demands  is  critical  to  USARAF’s  success.  The  scenario 
provides  the  background  and  rationale  for  deploying  troops  to  Nigeria,  and  includes 
a  regional  overview  of  neighboring  countries  and  Nigeria,  a  terrain  analysis,  a  human 
terrain  analysis,  and  a  time  line  of  significant  events  leading  to  a  deployment  of  US 
troops.  With  a  mission  to  provide  humanitarian  aid  and  enable  civilian  authority, 
this  research  concerns  the  placement  of  logistic  hubs  to  support  this  mission. 
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As  illustrated  in  Figure  8,  this  scenario  requires  the  design  of  a  theater  distribution 
network  to  sustain  nine  BCT-sized  units  at  geographically  dispersed  FOBs,  and  we 
restricted  supplies  to  flow  through  at  most  one  of  five  possible  APODs  and  at  most 
one  of  two  possible  SPODs.  We  further  restricted  the  scenario  to  use  at  most  two 
intermediate  distribution  centers.  The  choice  of  ndc  —  2  for  use  within  constraint 
(43)  was  based  on  the  number  and  disposition  of  units,  which  correlate  to  a  mission 
run  by  a  joint  task  force  with  two  subordinate  division  headquarters,  each  of  which 
would  have  an  aligned  sustainment  brigade  operating  a  regional  distribution  center. 
The  induced  network,  G[N,A],  was  developed  using  a  more  detailed  road  network 
map  from  the  World  Food  Programme  [43].  Considering  highways,  primary  roads, 
and  secondary  roads  for  the  transport  of  military  supplies,  the  resulting  network  has 
|iVj  =  135  and  |A|  =  414.  Within  this  structure,  we  selected  ne  =  2  medium  truck 
companies  for  e  =  1,2  based  on  subject  matter  expert  input  [21]. 


I  4  I  Possible  SPODs 

1.  Lagos 

2.  PortHarcourt 

Ilk  |  Possible  APODs 

1.  Abuja 

2.  Enugu 

3.  Kano 

4.  Lagos 

5.  PortHarcourt 

I  ><  I  FOBs  (~BCT-sized) 

1.  Abuja  (JTF  HQ  +  Fires  BDE  HQ) 

2.  Oshogbo  (MEB) 

3.  Makurdi  (BFSB) 

4.  Gombe  (IBCT/SBCT) 

5.  Mai duguri  (IBCT/SBCT) 

6.  Bauchi(CAB) 

7.  Kano  (IBCT/SBCT) 

8.  Birnin  Kebbi  (IBCT/SBCT) 

9.  Kaduna(CAB) 

FOBs  to  be  located  (n  ot  sh  own ) 

1 .  SPOD  (TSC  &  ESC  elements) 

2.  APOD  (TSC  &  ESC  elements) 

3.  Distribution  Center  #1 

4.  Distribution  Center #2 


Figure  8.  Scenario  -based  Disposition  of  Units  and  Possible  A/SPODs  in  Nigeria  [31] 


Figure  9  displays  the  undirected  road  network  used  in  the  scenario.  The  network  is 
only  a  portion  of  the  road  system  across  Nigeria.  Each  of  its  135  nodes  in  the  network 
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which  are  depicted  as  blue  circles  with  a  corresponding  number  inside,  is  located  at 
a  road  intersection  and  represent  the  only  locations  on  the  map  where  facilities  can 
be  located.  Each  arc  is  depicted  as  a  heavy  black  line  connecting  nodes.  As  stated  in 
Chapter  3,  we  assume  a  directed  network  that  ensures  that  we  can  characterize  the 
flow  over  an  arc  by  its  direction  of  movement. 


Figure  9.  The  Undirected  Road  Network,  G[N,A]  [31] 

We  estimated  the  demand  at  each  FOB  (6j)  using  a  theater  sustainment  planning 
reference  [9],  yielding  approximately  285  tons  of  weekly  demand  for  food,  fuel,  spare 
parts,  mail,  and  selected  other  classes  of  supply  at  each  such  location.  Weekly  demand 
at  to-be-located  DCs  ( bdc )  was  estimated  at  9  tons  of  similar  supplies.  In  the  absence 
of  specific  information  about  material  handling  equipment  (MHE)  restrictions,  we  set 
the  upper  bounds  on  throughput  at  APODs  (Cf),  SPODs  (C/),  and  DCs  (Ufc)  to  be 
1000,  3000,  and  1800,  respectively,  so  that  the  scenario  would  require  the  use  of  an 
APOD,  and  SPOD,  and  both  DCs.  We  estimated  distances  (djj)  to  the  nearest  five 
meters  using  a  road  network  map  [43],  and  we  assumed  upper  bounds  {utj)  of  1600, 
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800,  and  200  tons  for  highways,  primary  roads,  and  secondary  roads,  respectively.  To 
determine  the  risk  parameters  for  each  arc,  we  quantified  the  qualitative  geographic 
risk  assessments  within  Nigeria  by  the  United  Kingdom’s  Foreign  and  Commonwealth 
Office  (FCO)  [39].  We  used  the  FCO’s  2012  assessment  in  lieu  of  a  more  recent 
assessment  in  order  to  maintain  the  notional  nature  of  the  scenario,  and  we  weighted 
the  FCO’s  categories  of  “advise  against  all  travel”,  “advise  against  all  but  essential 
travel”,  and  “see  our  travel  advice  before  traveling”  with  linearly  scaled  values  of 
c\j  =  4  =  3,2,  and  1,  respectively.  Of  note,  as  shown  in  Figure  10  each  arc’s  risk 
classification  was  determined  by  the  highest  risk  category  for  any  region  through 
which  it  traverses,  and  we  did  not  assume  a  different  category  of  risk  for  supplies 
being  transported  at  different  echelons  of  flow. 
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4.2  Results 


In  this  section  we  analyze  results  using  seven  instances  with  different  weighting  vari¬ 
ations  corresponding  to  a  letter  (A,  B,  C,  D,  E,  F,  and  G).  We  analyze  both  the  ob¬ 
jective  function  components  and  the  decision  variables  with  respect  to  each  instance. 
As  described  in  Section  3.2,  the  objective  function  (37)  consists  of  three  measures 
of  performance;  the  accumulated  risk  to  supplies  flowing  through  the  network  (gi); 
the  total  distance  traversed  by  supplies  (g2);  and  the  maximum  supply-distance  load 
among  the  two  defined  echelons  of  transportation  (^3).  However,  it  is  true  that  ex¬ 
amining  the  entire  objective  function  value  at  optimality  is  not  particularly  useful. 
The  reason  is  because  the  units  on  qx ,  q2,  and  q3  are  not  the  same,  and  those  com¬ 
ponents  might  not  be  well  scaled.  Moreover,  a  component-wise  optimization  of  the 
objective  function  will  likely  indicate  that  they  are  in  tension.  These  reasons  indicate 
a  component- wise  analysis  of  the  objective  function  is  in  order,  which  we  accomplish 
by  examining  the  efficient  Pareto  frontier.  The  terms  risk,  distance ,  and  balance  are 
used  interchangeably  for  the  calculations  qi,  q2,  and  q3,  respectively. 

Because  the  multiple  objectives  have  different  units  of  measure  and  different  scales,  we 
examined  seven  different  instances  for  possible  corresponding  weights,  as  indicated 
in  Table  1,  each  of  which  yields  a  unique  non-dominated  solution  on  the  efficient 
Pareto  frontier.  To  identify  distinguishable  non-dominated  solutions  for  our  three- 
dimensional  Pareto  frontier,  we  first  consider  three  instances  (A-C)  for  the  given 
network,  where  we  weight  only  one  objective  function  component  at  a  time.  Next  we 
consider  two  objective  function  component  instances.  The  weights  for  Instances  D-F, 
a  scaled  weighting  for  the  combinations  of  (;))  objectives,  with  w±,w2,  and  w3  are  cal¬ 
culated  using  the  weighted  sum  method  based  on  the  optimal  solutions  to  Instances 
A-C.  Finally,  we  determine  the  weights  for  all  three  components  in  a  similar  manner 
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to  generate  and  solve  the  three  objective  function  component,  Instance  G  [3]. 

To  determine  the  two  weights,  w.t  and  Wj,  we  solve  the  set  of  equations: 

(q*)wi  =  ( q*)wj .  (58) 

Wi  +  Wj  =  1.  (59) 

where  q*  is  the  value  from  the  instance  (A,  B,  or  C)  where  Wi  =  1,  and  likewise 
for  q*.  This  weighted  sum  method  ensures  the  respective  weights  better  address 
the  relative  scaling  of  the  objective  function  components.  In  this  manner,  we  seek 
to  determine  appropriate  weights  that  account  for  possible  poor  scaling  among  the 
objective  function  components. 

Table  1.  w- values  for  Instances  A-G 


Instance 

w  i 

w2 

w3 

A 

1 

0 

0 

B 

0 

1 

0 

C 

0 

0 

1 

D 

0.98 

0.02 

0 

E 

0.94 

0 

0.06 

F 

0 

0.23 

0.77 

G 

0.93 

0.02 

0.05 

Given  comma-delimited  hies  for  the  node-  and  arc-specific  data  for  the  scenario’s 
network  shown  in  Figure  10  (or  any  other  instance),  we  encoded  in  Visual  Basic  a 
program  to  construct  the  formulation  within  Microsoft  Excel  and,  given  user-input 
weights  for  the  respective  objectives,  invoke  the  commercial  solver  CPLEX  (Version 
12.5)  to  solve  Problem  MTDNDP.  The  average  required  computational  effort  for  all 
seven  instances  was  less  than  two  minutes  on  a  computer  having  an  2.70  GHz  AMD 
Athlon  II  X2  processor  and  4.00  GB  of  RAM,  indicating  the  efficacy  of  our  formula¬ 
tion  and  solution  method  to  address  larger  instances.  Reported  in  Table  2  are  selected 
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elements  from  the  optimal  solutions  for  Instances  A-G;  the  first  three  columns  tabu¬ 
late  the  optimal  values  for  the  respective  objectives  corresponding  to  the  accumulated 
risk,  the  distance,  and  the  maximum  per  capita  workload  supported  by  transportation 
assets  at  a  given  echelon  (balance),  and  the  final  three  columns  indicate  the  number 
of  the  node  selected  for  use  as  an  APOD,  an  SPOD,  and  DCs.  For  the  last  row  of 
Table  2  we  denote  RPDj  to  be  the  relative  percentage  deviation  between  the  worst 
and  best  optimal  values  for  q% ,  i  =1,  2,  3,  over  Instances  A-G. 


Table  2.  Selected  Optimal  Decision  Variable  Values,  Instances  A-G 


Instance 

<1\ 

Q2 

<h 

APODs 

SPODs 

DCs 

A 

21,666 

1,536,885 

678,868 

110 

1 

2,  107 

B 

28,457 

1,213,825 

543,838 

110 

1 

2,  111 

C 

44,400 

1,486,964 

371,741 

110 

10 

40,  111 

D 

22,968 

1,400,515 

563,183 

110 

1 

2,  114 

E 

26,136 

1,524,649 

381,162 

110 

1 

28,  92 

F 

42,457 

1,416,260 

384,700 

110 

10 

40,  111 

G 

25,568 

1,450,775 

400,425 

110 

1 

28,  92 

RPDj 

105% 

27% 

83% 

Intuitively,  each  of  the  respective  objectives  attained  its  lowest  value  when  only  it  was 
weighted,  as  is  the  case  for  Instances  A-C.  For  example,  the  weighting  for  Instance 
A  (wi,w2,w3)  =  (1,0,0)  returns  the  lowest  value  for  q\  for  all  seven  instances.  With 
this  weighting,  only  the  risk-related  objective  is  forced  to  minimize.  The  distance- 
related  and  balance-related  objectives  are  allowed  to  fluctuate  in  order  to  minimize 
the  risk-related  objective  and  are  not  the  lowest  values  for  the  seven  instances.  The 
weighting  for  Instances  B  and  C  return  the  lowest  values  for  distance-related  and 
balance-related  objectives,  respectively,  in  the  same  manner. 

Despite  the  common  use  of  the  APOD  at  Kano  (i.e.,  Node  110)  for  all  instances,  the 
optimal  solution  for  each  instance  was  unique  with  regard  to  other  sustainment  node 
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emplacements  and/or  supply  flow  routing.  Although  other  non-dominated  solutions 
may  exist  on  the  efficient  Pareto  frontier,  we  did  not  seek  to  find  them  via  an  ex¬ 
haustive  search;  we  intend  for  this  model  merely  to  provide  several  distinguishable 
alternative  (nondominated)  solutions  for  a  decision  maker  to  consider.  Of  note,  the 
RPD  between  the  worst  and  best  optimal  values  for  q*Xl  q*2  and  q\  over  Instances  A-G 
were  respectively  105%,  27%,  and  83%.  This  result  indicates  a  higher  sensitivity  of 
both  the  risk-related  and  the  workload-balancing-related  objectives  to  variations  in 
the  objective  weights  for  this  scenario. 

Although  the  RPD  for  the  distance-related  objective  was  only  27%,  this  result  is 
not  insignificant;  it  corresponds  to  a  possible  increase  in  road  mileage,  wear-and-tear 
on  vehicles,  and  fuel  usage  of  up  to  27%,  depending  on  the  non  dominated  solution 
selected  for  implementation.  The  RPD  for  the  balance-related  objective  was  83% 
between  the  highest  and  lowest  solutions  which  corresponds  to  a  significant  imbal¬ 
ance  between  the  work  load  of  units  designated  to  haul  Echelon  1  and  Echelon  2 
flow.  With  such  a  strong  discrepancy  in  work  load,  a  decision  maker  may  consider 
changes  to  the  task  organization  to  prevent  both  personnel  and  vehicles  being  over 
utilized.  The  RDP  for  the  risk-related  objective  was  105%.  This  corresponds  to  units 
hauling  supplies  being  as  much  as  twice  as  vulnerable  to  a  threat  depending  on  the 
non  dominated  solution  selected. 

Figure  11  illustrates  the  seven  non-dominated  points  on  the  Pareto  frontier  corre¬ 
sponding  to  Instances  A-G.  the  Pareto  frontier  consists  of  non-dominated  solutions. 
The  continuous  surface  that  includes  these  points  looks  somewhat  like  an  irregular- 
shaped  bowl  made  of  facets.  We  seek  to  minimize  the  objective  function  components. 
The  points  closer  to  the  origin  (points  ‘below’  or  ‘outside’  the  bowl)  are  not  feasible. 
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In  contrast  and  in  the  context  of  Figure  11,  feasible  points  exist  that  are  closer  to  the 
reader  (‘above’  or  ‘inside’  the  bowl).  However,  these  points  are  dominated  by  points 
on  the  the  Pareto  frontier  because  there  exists  a  point  on  the  frontier  that  has  the 
risk,  distance,  and  balance  values  at  least  as  low  as  a  dominated  solution,  with  one 
of  the  component  values  being  strictly  lower. 


Figure  11.  Pareto  Frontier  Representation  of  the  Optimal  Solutions  for  Instances  A-G 

Beyond  a  comparison  of  the  values  for  the  objective  function  components,  there  is 
merit  to  examining  the  other  optimal  objective  function  values  for  the  seven  instances. 
Such  a  comparison  seeks  to  identify  whether  the  optimal  solutions  are  distinguishable 
[40].  Moreover,  it  may  identify  commonalities  among  the  optimal  solutions  as  well, 
thereby  generating  insights  on  solution  characteristics  that  are  sensitive  to  relative 
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component  weightings.  Thus,  while  Instances  A-G  are  all  on  the  Pareto  frontier, 
we  show  that  Instances  A,  B,  and  C  provide  the  best  solutions  when  considering 
trade-offs  between  the  risk-related,  distance-related  and  balance-related  objectives. 
Instances  D,  E,  F,  or  G  do  not  allow  for  the  minimum  value  of  any  combination 
of  measure  of  performance  for  the  the  risk-related,  distance-related,  and  balance- 
related  objectives.  For  example,  when  comparing  Instances  D  and  A,  more  balance 
is  achieved  between  Echelon  1  and  Echelon  2  flow  in  Instance  D,  as  well  as  a  lower 
distance  traveled  in  Instance  D.  However,  Instance  A  still  returns  a  smaller  risk  value. 

To  illustrate  elements  of  the  characteristics  of  our  solutions,  we  depict  in  Figures  12 
and  13,  the  optimal  solutions  to  Instances  A  and  C,  overlaid  on  a  map  that  depicts  the 
geographic  risk  assessments,  respectively  [39].  Within  each  figure,  we  depict  APODs, 
SPODs,  and  DCs  as  red,  orange,  and  yellow  nodes,  with  all  other  nodes  being  blue 
in  color;  and  we  represent  Echelon  1  and  2  flows  by  coloring  the  arcs  on  which  they 
are  shipped  as  green  and  purple,  respectively,  with  some  arcs  supporting  flow  at 
both  echelons.  Using  the  graphic  solution  for  each  instance  along  with  Table  2,  the 
relationships  among  the  instances  can  be  explored. 
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Node 

^MEchelon  1 
Echelon  2 
O  SPOD 
•  APOD 
O  DC 


Advise  against  all  travel 

Advise  against  all  but  essential  travel 

See  our  travel  advice  before  travelling 


Figure  12.  Optimal  Solution  for  Instance  A  [39] 


Node 
Echelon  1 
Echelon  2 
O  SPOD 
•  APOD 
O  DC 


Advise  against  all  travel 

Advise  against  all  but  essential  travel 

See  our  travel  advice  before  travelling 


Figure  13.  Optimal  Solution  for  Instance  C  [39] 


Visible  in  Figure  12  is  the  most  risk-averse  solution;  it  utilizes  Lagos  (i.e.,  Node  1)  as 
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the  SPOD  in  lieu  of  Port  Harcourt  (i.e. ,  Node  10),  and  it  uses  very  few  of  the  roads 
in  the  high  threat  areas  of  the  country.  However,  it  also  emplaces  the  distribution 
centers  very  close  to  the  selected  APOD  and  SPOD,  yielding  an  imbalanced  load 
between  the  two  echelons.  The  supply-distances  supported  by  the  second  echelon 
are  7.5  times  as  great  as  the  first  echelon,  indicating  that  a  change  in  transportation 
resource  allocation  will  be  necessary  to  implement  this  solution. 

In  contrast,  the  optimal  solution  for  Instance  C  perfectly  balances  the  workload  be¬ 
tween  the  two  echelons  by  emplacing  one  of  the  the  distribution  centers  in  the  more 
centrally-located  town  of  Lokoja  (i.e.,  Node  40).  However,  the  greatest  risk-based 
objective  is  attained  via  the  solution  to  Instance  C,  wherein  Port  Harcourt  is  utilized 
as  an  SPOD  and  flows  are  routed  throughout  the  high  threat  areas  without  regard 
to  the  risk  incurred  to  convoys.  This  graphic  finding  is  supported  by  Table  2  where 
the  highest  value  for  q\  was  attained  when  (wi,w2,w3)  =  (0,0,1)  (i.e.,  Instance  C), 
and  the  highest  value  for  q was  attained  when  (ity,  w2,  w^)  =  (1,0,0)  (i.e.,  Instance 
A),  indicating  the  tension  between  these  two  objectives  for  this  scenario.  Moreover, 
both  Instances  A  and  C  do  not  address  the  distance-based  metric,  resulting  in  the 
respective  (ordinally  ranked)  first  and  third  worst  values  for  q*2. 

Should  a  decision  maker  be  concerned  only  with  minimizing  the  total  distance  tra¬ 
versed  by  all  supplies,  consider  the  solution  to  Instance  B  as  depicted  in  Figure  14 
[39].  In  order  to  attain  the  shortest  supply  routings,  this  solution  incurs  a  31%  in¬ 
crease  in  the  risk-based  metric  compared  to  Instance  A  by  routing  supplies  through 
higher-risk  regions,  and  it  allows  for  a  second-echelon  workload  that  is  8.6  times  the 
first-echelon  workload  by  locating  the  two  distribution  centers  close  to  the  selected 
APOD  and  SPOD. 
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Node 
Echelon  1 
Echelon  2 
O  SPOD 
•  APOD 
O  DC 


Advise  against  all  travel 

Advise  against  all  but  essential  travel 

See  our  travel  advice  before  travelling 


Figure  14.  Optimal  Solution  for  Instance  B  [39] 


4.3  Sensitivity  Analysis 

Acknowledging  that  our  quantification  of  risk  via  a  linear  metric  is  not  the  only 
suitable  method,  we  examined  the  sensitivity  of  our  solutions  to  the  imposition  of 
a  quadratic  risk-based  metric,  with  values  of  {1,4,9}  in  lieu  of  {1,2,3}.  Using  the 
same  weighted-sum  method  to  determine  the  objective  weights  when  more  than  one 
objective  is  considered,  we  calculated  the  new  weights  for  all  seven  instances  as  shown 
in  Table  3. 

We  then  used  the  same  comma-delimited  files  for  the  node-  and  arc-specific  data, 
Visual  Basic  program,  and  CPLEX  (Version  12.5)  solver  to  obtain  optimal  solutions 
for  all  seven  instances  as  shown  in  Table  4. 

Similar  to  the  linearly  weighted  optimal  solutions,  the  quadratic  risk-based  optimal 
solutions  attain  their  lowest  value  when  only  each  objective  is  individually  weighted, 
as  in  the  cases  for  Instances  A’-C’.  Also,  the  placement  of  APODs,  SPODs,  and  DCs 
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Table  3.  w-values  for  Instances  A’-G’  with  Quadratic  Weighting 


Instance 

Wi 

w2 

w3 

A’ 

1 

0 

0 

B’ 

0 

1 

0 

C’ 

0 

0 

1 

D’ 

0.98 

0.02 

0 

E’ 

0.93 

0 

0.07 

F’ 

0 

0.23 

0.77 

G’ 

0.91 

0.02 

0.07 

Table  4.  Selected  Optimal  Decision  Variable  Values  with  Quadratic  Weighting,  In¬ 
stances  A’-G’ 


Instance 

Qi 

92 

<h 

APODs 

SPODs 

;  DCs 

A’ 

27,374 

1,579,285 

700,068 

110 

1 

2,  107 

B’ 

49,625 

1,213,825 

543,838 

110 

1 

2,  111 

C’ 

87,907 

1,486,964 

371,741 

110 

10 

40,  111 

D’ 

29,157 

1,411,490 

616,170 

110 

1 

2,  107 

E’ 

35,495 

1,529,231 

382,308 

110 

1 

28,  92 

F’ 

84,119 

1,416,260 

384,700 

110 

10 

40,  111 

G’ 

35,778 

1,475,800 

390,250 

110 

1 

28,  92 

RPD* 

221% 

30% 

88% 

is  the  same  between  the  linear  and  quadratic  optimal  solutions  with  the  exception  of 
Instance  D  and  D’.  Differing  notably  between  the  linear  and  quadratic  risk-based  op¬ 
timal  solutions  are  the  RPD  values  for  each  of  the  objectives.  The  RPD  between  the 
worst  and  best  optimal  values  for  q\ ,  q2  and  q 3  over  the  instances  with  the  alternative 
risk  quantification  were  221%,  30%,  and  88%,  respectively,  an  increase  of  116%,  3%, 
and  5%  over  the  linear  risk-based  optimal  solutions  displayed  in  Table  2.  A  compari¬ 
son  of  the  optimal  objective  function  components  for  the  linear  and  quadratic  models 
is  shown  in  Table  5.  The  third  column  displays  the  percent  change  (PC)  between  the 
linear  and  quadratic  objectives  for  each  instance 


The  risk-related  objective  experienced  in  some  cases  almost  100%  PC  between  the 
linear  and  quadratic  models,  indicating  that  the  risk-related  objective  is  highly  sen- 
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Table  5.  Linear  and  Quadratic  Objective  Function  Component  Comparison 


Instances 

Qi 

Lin  Quad 

PC  (%) 

<72 

Lin  Quad 

PC  (%) 

13 

Lin  Quad 

PC  (%) 

A,  A' 

21,666 

27,374 

26 

1,536,885 

1,579,285 

2.8 

678,868 

700,068 

3.12 

B,  B’ 

28,457 

49,625 

74 

1,213,825 

1,213,825 

0.0 

543,838 

543,838 

0.00 

C,  C’ 

44,400 

87,907 

98 

1,486,964 

1,486,964 

0.0 

371,741 

371,741 

0.00 

D,  D’ 

22,968 

29,157 

27 

1,400,515 

1,411,490 

0.8 

563,183 

616,170 

9.41 

E,  E' 

26,136 

35,495 

36 

1,524,649 

1,529,231 

0.3 

381,162 

382,308 

0.30 

F,  F’ 

42,457 

84,119 

98 

1,416,260 

1,416,260 

0.0 

384,700 

384,700 

0.00 

G,  G’ 

25,568 

35,778 

40 

1,450,775 

1,475,800 

1.7 

400,425 

390,250 

-2.54 

APCj 

57 

0.8 

1.5 

sitive  to  change.  The  largest  average  percent  change  (APC*)  occurred  in  the  risk- 
related  objective  indicating  that  risk  is  the  most  sensitive  objective.  Regarding  the 
effect  of  the  quadratic  risk  scale,  we  illustrate  in  Figure  15  the  optimal  solution  when 
(wi,W2,ws)  =  (1,0,0)  and  note,  compared  to  the  solution  depicted  in  Figure  12,  the 
rerouting  of  the  flow  of  supplies  to  reduce  transit  through  the  high  threat  areas. 


Node 

^MEchelon  1 
Echelon  2 

•  SPOD 

•  APOD 
O  DC 


Advise  against  all  travel 

Advise  against  all  but  essential  travel 

See  our  travel  advice  before  travelling 


Figure  15.  Optimal  Solution  for  Instance  A  with  Quadratic  Risk  Quantification  [39] 


The  distance-  and  balance-related  objectives  experienced  0.8%  and  1.5%  difference 
between  the  linear  and  quadratic  models,  respectively.  This  result  further  supports 
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the  identification  of  the  insensitivity  of  the  distance-based  and  workload-balancing 
objectives  to  the  alternative  quantification  of  risk. 
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V.  Conclusion 


Background 

“Logistics  ...  as  vital  to  military  success  as  daily  food  is  to  daily  work.” 

Capt  Alfred  Thayer  Mahan,  USN  [25] 

This  chapter  provides  the  research  conclusion  from  combining  our  developed  model 
and  a  representative  scenario.  First,  in  Section  5.1  we  provide  a  summary  of  the 
problem  statement  and  methodology  to  solve  the  problem.  Next,  in  Section  5.2  we 
discuss  possible  further  research. 

5.1  Conclusions 

This  work  examined  the  problem  of  designing  a  military  theater  distribution  network, 
wherein  a  decision  maker  selects  Air  and  Sea  Ports  of  Debarkation  and  intermediate 
logistical  distribution  centers,  through  which  we  route  military  supplies  over  a  directed 
transportation  network  to  meet  demands.  Considering  alternative  objectives  to  min¬ 
imize  the  total  risk  encountered  by  transporting  supplies,  the  total  distance  traveled 
by  supplies,  and  the  maximum  per  capita  workload  supported  by  transportation  as¬ 
sets  at  a  given  echelon  (i.e. ,  port-to-distribution  center  versus  distribution  center  to 
demands),  we  set  forth  a  multi-objective,  mixed- integer  linear  programming  (MILP) 
formulation  to  solve  this  problem.  For  a  given  instance,  we  developed  an  automated 
tool  that  generates  the  MILP  formulation  and  solves  it  using  the  commercial  solver 
CPLEX.  To  demonstrate  the  efficacy  of  the  model,  we  applied  it  to  a  representative 
instance  and  generated  seven  nondominated  solutions  on  the  efficient  Pareto  frontier. 
Computational  tests  indicated  the  relative  variation  of  the  three  objectives  among 
the  instances  examined  and  identified  two  of  the  objectives  as  being  in  tension  for 
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the  given  scenario.  Moreover,  a  successive  analysis  indicated  that  the  distance-based 
and  load-balancing  objectives  were  relatively  insensitive  to  an  alternative  metric  for 
quantifying  risk. 


5.2  Future  Research 

In  Section  1.3  we  outlined  four  objectives  for  this  research.  The  first  research  objective 
was  met  by  developing  an  MILP  for  the  movement  of  supplies  through  a  transporta¬ 
tion  network.  The  second  research  objective  was  met  by  encoding  in  Visual  Basic  a 
program  to  construct  the  formulation  within  Microsoft  Excel  and  invoke  a  commer¬ 
cial  solver  CPLEX  (Version  12.5)  to  solve  the  Problem  MTDNDP.  The  third  research 
objective  was  met  by  using  a  scenario-based  case  study  to  demonstrate  the  model. 
The  fourth  research  objective  discussed  examining  extensions  to  the  analysis.  The 
remainder  of  this  section  discusses  several  areas  for  model  improvement  and  areas  for 
further  research. 

There  are  at  least  four  areas  for  model  improvement.  First,  an  upper  bound  should 
be  imposed  on  the  total  distance-based  workload.  This  will  ensure  that,  if  a  decision 
maker  seeks  to  balance  the  workload  between  echelons  of  sustainment  transport,  it 
does  not  come  at  the  cost  of  irrationally  increasing  the  total  distances  traveled.  Sec¬ 
ond,  overall  testing  could  be  improved  by  using  actual  node  throughput  restrictions 
at  APODs,  SPODs,  and  DCs  based  on  the  availability  of  MHE.  The  availability  of 
MHE  varies  on  the  task  organization  and  number  of  CSSBs  and  this  research  made 
assumptions  about  the  task  organization.  Model  users  will  determine  assets  available 
for  mission  requirements  and  have  a  refined  understanding  of  MHE  available.  Third, 
improve  the  intuitive  understanding  of  q2  and  q 3  by  converting  them  from  supply- 
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distances  and  supply-distances  per  truck  company  to  distance  per  unit  of  supply  and 
average  distance  per  truck  company,  respectively.  Fourth,  the  model  and  solution 
methodology  should  be  applied  to  additional  scenarios  to  further  assess  its  efficacy. 

One  suggested  area  for  further  research  is  the  consideration  of  the  reverse  supply 
chain  for  the  evacuation  of  casualties  and  materials.  In  the  context  of  steady-state 
operations  casualties,  whether  due  to  hostile  fire  or  non-combat  injuries,  will  occur. 
While  some  local  medical  facilities  will  have  the  capacity  to  treat  casualties,  some 
casualties  will  require  evacuation  to  receive  higher  levels  of  care.  This  model  could 
be  modified  to  support  the  movement  of  casualties  within  the  country  or  out  of  the 
country.  Assuming  there  is  no  intent  for  the  United  States  Government  to  maintain 
a  long  term  presence  in  a  country,  all  material  deployed  into  the  country  will  have  to 
be  redeployed  to  its  home  station.  In  the  context  of  a  redeploying  force,  the  model 
could  be  modified  to  support  a  reverse  supply  chain. 

Another  area  worthy  of  further  research  is  to  allow  direct  flow  of  supplies  from  a 
located  POD  to  a  FOB  which  could  be  either  through  ground  or  air  lines  of  commu¬ 
nication.  Reasons  for  allowing  direct  flow  from  a  POD  to  a  FOB  include  but  are  not 
limited  to  an  operational  consideration  that  necessitates  a  rapid  flow  of  supplies  into 
a  FOB  or  a  located  POD  closer  (time  and/or  distance)  to  a  set  of  FOBS  than  it  is  to 
a  located  DC.  In  the  case  of  direct  flow  of  supplies  via  air  from  a  POD  to  a  FOB  two 
additional  reasons  are  to  forgo  the  risk  of  ground  high  threat  areas  and  non-existent 
ground  lines  of  communication.  Changes  to  the  model  to  accommodate  direct  flow  of 
supplies  from  a  located  POD  to  a  FOB  would  be  minimal.  A  third  echelon  of  flow, 
which  describes  direct  flow  from  a  located  POD  to  a  FOB,  would  be  introduced.  The 
model  would  have  to  allow  Echelon  3  flow  out  of  a  POD  and  into  a  FOB.  The  amount 
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of  Echelon  3  flow  out  of  a  POD  would  be  bounded  the  summation  of  the  demand  of 


the  FOBs  that  it  are  closer  to  the  POD  than  a  located  DC.  The  amount  of  Echelon  3 
flow  into  a  FOB  would  be  bounded  by  the  demand  that  FOB  requires.  A  derivative 
to  the  direct  flow  of  supplies  from  a  located  POD  to  a  FOB  is  adding  additional  arcs 
to  the  network  to  represent  intra-theater  flow  between  airports,  or,  as  appropriate, 
transportation  via  railway. 

We  used  the  Army’s  Command  and  General  Staff  College’s  Theater  Sustainment  Bat¬ 
tle  Book  and  easily  accessible  maps  with  varying  degreees  of  fidelity  to  determine  arc 
parameters  for  the  network.  These  arc  parameters  were  valid  for  our  proof  of  prin¬ 
ciple,  however  use  of  this  model  in  the  future  should  involve  better  network  data  for 
the  country  or  region  that  is  analyzed.  Some  network  data  would  necessarily  classify 
this  analysis,  however  there  is  available  open  source  data  that  would  make  the  results 
valuable  to  a  decision  maker.  Finally,  a  more  general  network  flow  model  should  be 
examined,  representing  the  possible  loss  of  supplies  either  along  arcs  or  at  nodes  due 
to  pilferage. 
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Appendix  A.  Notional  Nigeria  Tactical  Scenario 


This  map  provides  an  overview  of  the  Nigerian  international  borders,  coast  line,  major 
cities,  and  major  roads.  The  road  network  is  the  primary  method  for  transporting 
supplies  from  A/SPODs  to  FOBs.  The  country  of  Nigeria  in  West  Africa  is  the  most 
populace  African  nation  and  contains  an  abundance  of  some  of  the  worlds  riches.  Its 
neighbors  are  the  Republic  of  Benin  in  the  west,  Chad  and  Cameroon  in  the  east,  and 
Niger  in  the  north.  Its  coast  in  the  south  lies  on  the  Gulf  of  Guinea  in  the  Atlantic 
Ocean  which  gives  it  access  to  extensive  trade  opportunities. 


UNCLASSIFIED  MODELING  SCENARIO 

CJTF-Relief  CONOPS  AAtfl' 


1  The  AFIT  of  Today  is  the  Air  Force  of  Tomorrow. 
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The  description  below  gives  a  summary  of  Nigeria’s  neighbors.  Nigeria’s  stability 
is  affected  by  events  in  neighboring  countries.  Nigerias  neighbors  vary  in  levels  of 
economic,  social,  and  political  stability.  As  a  nation  formed  by  colonial  powers  in 
the  18th  century,  Nigerias  borders  are  at  times  non-existent  to  tribal  communities 
that  extend  beyond  borders  into  neighboring  countries.  Along  the  border  regions, 
and  often  times  in  the  interior,  loyalty  to  tribe  supersedes  loyalty  to  the  Nigerian 
government. 


UNCLASSIFIED  MODELING  SCENARIO 

Regional  Overview 

Neighboring  Countries 


The  AFIT  of  Today  is  the  Air  Force  of  Tomorrow. 

Benin 

Gained  independence  from  France  on  1  August  1960.  Relatively  stable  republic  government.  Has  three  branches  of 
government  with  political  power  resting  with  a  prime  minister.  A  majority  of  the  people  live  along  the  southern 
coastline.  Half  of  the  population  is  split  between  Roman  Catholic  and  Muslim.  The  economy  is  dependent  on 
subsistence  agriculture,  cotton  production,  and  regional  trade. 


^FIT 

Air  Foroc  Institute  or  Tcohnoloiiy 


Niger 

Gained  independence  from  France  on  1  August  1960.  Since  independence  the  government  was  experienced 
numerous  coups.  A  presidential  election  was  help  in  2011  with  no  clear  winner.  Consequently,  the  current  president 
rules  with  little  authority.  More  than  80%  of  the  country  is  covered  by  the  Sahara  desert  (central  and  north)  and  most 
people  live  in  the  south  or  west  of  the  country.  The  population  is  mostly  80%  Muslim.  It  is  a  landlocked  nation  whose 
economy  is  concentrated  around  subsistence  farming.  It  is  a  developing  country  and  consistently  ranks  low  on  the 
United  nations’  Human  Development  Index. 


Chad 

Gained  independence  from  France  on  11  August  1960  and  endured  three  decades  of  civil  war  and  an  invasion  by 
Lybia  before  peace  was  restored  in  1990.  By  the  constitution  the  president  is  elected  every  5  years  and  contains 
almost  all  power.  More  than  54%  of  the  country  is  Muslim  and  34%  Christian.  By  some  indexes  it  is  considered  a 
failed  state  and  consistently  ranks  as  one  of  the  poorest  countries  on  the  planet.  Although  there  are  oil  reserves  in  the 
country,  due  to  corruption  over  80%  of  the  country  relies  on  subsistence  farming. 


Cameroon 

Gained  independence  from  France  on  1  January  1960.  It  is  considered  a  dominant-party  president  republic  with  an 
authoritarian  president  selected  every  seven  years  and  enjoys  social  and  political  stability.  Cameroon  has  developed 
roads,  railways,  agriculture  systems,  and  large  petroleum  and  timber  industries.  Indigenous  and  Christian  beliefs 
each  make  up  40%  of  the  population  with  the  remaining  20%  being  Muslim.  Agriculture  and  natural  resources 
dominate  the  economy  and  Cameroon  is  an  exporter  of  its  goods  to  industrialized  nations.  Still,  30%  of  the 
population  is  unemployed  and  one-third  of  the  population  lives  below  international  poverty  standards. 

Air  University:  The  Intellectual  and  Leadership  Center  of  the  Air  Force 

Aim  High. ..Fly  ■  Fight  -  Win  s 


65 


The  description  below  describes  Nigeria’s  geography  and  topography  in  a  military 
context.  Nigeria  is  about  twice  the  size  of  California  and  varies  in  terrain  from 
tropical  jungle  in  the  south  to  a  dry  high  plateau  in  the  north.  Lines  of  communication 
throughout  Nigeria  are  marginal  and  even  poor  in  rural  areas. 

UNCLASSIFIED  MODELING  SCENARIO 

Regional  Overview  ^^FIT 

Terrain  «,« «...  ...  r..«.ou.„ 

The  AFIT  of  Today  is  the  Air  Force  of  Tomorrow. 


•  Overview 

•  The  total  area  is  about  twice  the  size  of  California  and  is  the  world’s  32nd  largest  country.  Nigeria  has  five 
major  geographic  regions:  a  low  coastal  zone  along  the  Gulf  of  Guinea;  hills  and  low  plateaus  north  of  the 
coastal  zone;  the  Niger-Benue  river  valley;  a  broad  stepped  plateau  stretching  to  the  northern  border  with 
elevations  exceeding  4000  feet;  and  a  mountainous  zone  along  the  eastern  border,  which  includes  the 
country's  highest  point,  Chappal  Waddi  (8000  feet).  The  Niger  River  (from  west)  and  Benue  River  (from 
east)  meet  in  the  center  of  Nigeria  and  flow  southward  into  the  Gulf  of  Guinea. 

•  Observation  and  Fields  of  Fire: 

•  Limited  in  the  rain  forests  of  the  south  and  in  all  urban  areas  (Lagos  - 1 1M,  Abuja  -  2M,  and  Kano  -  3.5M) 

•  In  the  central  savannah  and  northern  highlands  there  is  less  vegetation  and  visible  distances  are  measured 
in  miles. 

•  Obstacles: 

•  In  the  south  vehicular  traffic  is  limited  to  roads  due  to  rain  forest  and  mountains. 

•  In  the  central  and  northern  sections  of  the  country  off-road  is  possible.  In  the  central  area  the  convergence 
of  the  Niger  and  Benue  Rivers  restrict  some  off-road  movement.  In  the  north  there  are  steep  valleys  that 
can  restrict  some  cross  country  movement. 

•  Cover  and  Concealment: 

•  Excellent  in  the  south.  Entire  villages  and  small  settlements  that  are  not  documented  by  the  government  are 
plentiful. 

•  The  heavily  vegetated  areas  of  the  central  savannah  along  the  major  river  provide  excellent  hiding  locations. 

•  The  northern  areas  offer  little  vegetation,  however  wadis  and  broken/uneven  terrain  are  plentiful. 

•  Avenues  of  Approach: 

•  Almost  67,000  miles  of  surfaced  road  exist,  but  rarely  maintained  and  can  be  treacherous  at  high  speeds. 

All  major  cities  are  connected  with  some  4-lane  highways.  The  railway  system  runs  mostly  north-south  and 
is  in  poor  repair.  There  are  two  major  airports  and  5  international  airports. 

Air  University:  The  Intellectual  and  Leadership  Center  of  the  Air  Force 
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The  description  below  describes  Nigerias  human  terrain  in  a  military  context.  Nigeria 
is  split  along  two  lines;  the  northern  dominated  Muslim  population  and  the  south¬ 
ern  dominated  Christian  population.  There  are  four  major  tribes  that  align  with 
geographic  and  religious  lines.  Competition  for  land  and  other  natural  resources  be¬ 
tween  these  groups  has  a  history  of  violence,  including  a  Civil  War  in  1967. 


UNCLASSIFIED  MODELING  SCENARIO  _ 

W  AO  OE  &  PMESII  +  PT  ^FIT 

F....  IH.TITUTI  ..  TlOnNOLOBY 


The  AFIT  of  Today  is  the  Air  Force  of  Tomorrow. 


OE: 

•  Militias  are  hidden  within  the 
population  but  are  comprised  of 
local  ledership. 

•  Militias  are  dispersed  but  mainly 
in  the  northern  region. 

•  To  survive  a  majority  of  the 
population  has  chosen  sides. 


Political 

Split  along  tribal,  regional,  and  religious  lines.  Muslim 
north  are  from  two  tribes  while  the  south  has  three 
tribes.  Southern  Nigerians  have  neglected  northern 
concerns 

Military 

Northerners  see  the  military  as  heavy  handed  against 
Muslims  and  do  not  support  it. 

Economic 

Revenue  from  rich  natural  resources  in  the  south  have 
been  squandered  due  to  corruption. 

Social 

Hausa  and  Fulani  in  the  north.  Yoruba  in  the  south¬ 
west  Ijaws  in  the  south  and  Igbo  in  the  south-east 
Nigerian  Civil  War  (1967  - 1970)  brought  secessionist 
Igbo's  back  into  Nigeria  Over  1 M  killed 

Infrastructure 

Numerous  road  networks  in  poor  condition  due  to 
neglect.  Rail  system  in  marginal  condition  and  runs 
north  to  south. 

Information 

Muslim  militias  (predominantly  Boko  Haram)  in  the 
north  benefit  from  world-wide  Muslim  support  Some 
foreign  fighters  have  bee  found  in  Nigeria 

Physical 

Terrain 

Complex  terrain  in  the  south  limits  ISR  capabilities. 
Central  and  northern  regions  are  more  open 

Time 

Muslim  militias  know  that  the  government  will  not  curb 
violence  by  Christian  militias.  Time  is  on  the  side  of  the 
Boko  Haram  as  southern  patience  is  waning 

15  August  2014 
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The  description  below  provides  a  historical  accounting  of  how  US  forces  became  in¬ 
volved  in  the  Nigerian  conflict.  The  deployment  of  US  forces  to  Nigeria  is  traced 
directly  to  unrest  starting  in  2014.  The  roots  of  this  unrest  dates  back  decades  to 
Nigerias  independence  (1960)  and  even  centuries  to  the  colonial  period  (1800-1960). 


UNCLASSIFIED  MODELING  SCENARIO 

\  J  Road  to  War 

Air  Fcrhe  Ihbti.ute  or  Teshulqik 

The  AFIT  of  Today  is  the  Air  Force  of  Tomorrow. 

I960:  Gained  independence  from  the  United  Kingdom.  Underlying  tensions  between 
south  (Christian,  educated,  more  affluent,  natural  resource  rich)  and  north  (Muslim,  less 
educated,  few  natural  resources). 

1999:  After  16  years  of  military  rule  a  new  constitution  and  civilian  government  put  in 
place.  Brutal  military  dictatorships  were  harsh  on  the  north. 

2002:  Militant  Islamic  group  Boko  Haram  founded  by  Mohammed  Yusuf. 

2003:  The  first  violent  Boko  Haram  attacks  were  against  police  stations. 

2004:  More  attacks  against  government.  Hundreds  of  Boko  Haram  members  rounded  up 
and  jailed.  Christian  militias  kill  over  600  Muslims  in  north. 

2004  -  2008:  Boko  Haram  limits  their  attacks  and  builds  strength. 

2009  -  1010:  Boko  Haram  is  able  to  attack  prisons  and  free  over  1000  jailed  members. 

Conduct  significant  attacks  all  across  the  north.  Christian  militias  organize  in  the  north. 

201 1 :  Boko  Haram  kill  over  200  people  and  begin  to  move  their  attacks  southward. 

2012:  Over  500  deaths  attributed  to  Boko  Haram  mostly  in  the  north  and  central  areas. 
President  Goodluck  Jonathan  calls  for  talks  with  Boko  Haram  which  quickly  fall  apart. 

2013:  Boko  Haram  limits  their  attacks  and  builds  strength.  Kidnappings  are  rampant. 

2014:  To  date  over  1200  people  have  been  killed  in  Boko  Haram  attacks.  The  two  most 
notable  were  the  attacks  on  Bamboru  and  Ngala  which  left  over  300  dead  and  the 
abduction  of  over  300  young  girls  from  Chibok.  No  data  on  Christian  caused  deaths. 

2014:  Thousands  of  Nigerians  in  the  north  and  central  have  fled  their  homes  and  headed 
to  neighboring  countries  or  to  the  south.  Reception  in  the  south  has  been  met  with 
violence  as  overcrowding  has  increased  tensions. 


15  August  2014 
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The  description  below  displays  all  US  military  forces  that  are  used  in  the  operation. 
The  wire  diagram  further  shows  the  US  Army  forces  that  participate  in  the  opera¬ 
tion.  US  military  operations  in  Nigeria  are  conducted  under  a  JTF.  Although  the 
operation  is  conducted  by  US  combat  units,  the  mission  is  a  humanitarian  assistance 
mission. 


UNCLASSIFIED  MODELING  SCENARIO 

CJTF-Relief  Task  Organization 
and  Mission 

The  AFIT  of  Today  is  the  Air  Force  of  Tomorrow. 


Am  Force  Institute  or  Teomnoloiiy 


3rd  ID  Mission;  On  order  31 D  deploys  to 
secure  key  infrastructure  within  the 
Nigerian  AO,  provide  humanitarian  aid, 
and  enable  civilian  authority  IOT  to 
prevent  humanitarian  crisis  from 
overwhelming  host  nation  capabilities  and 
becoming  regional  conflict.  On  order 
transition  control  of  humanitarian 
operations  to  host  nation  forces  and 
NGOs  and  redeploy. 


rnl 

2  x  SBCT 

1  xTSC 

2  x  IBCT 

1  xESC 

1  x  FIRES 

2  x  Sust  Bde 

1  xMEB 

1  x  CSSB 

2  x  CAB 

4  x  Med  Truck  Co  Cargo  (22T) 

1  x  BFSB 

2  x  Med  Trcuk  Co  POL 

1  x  HET  Co 
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The  description  below  provides  a  general  timeline  and  description  for  US  military 
operations  in  Nigeria.  The  times  are  for  planning  purposes  only  and  may  vary  de¬ 
pending  on  operational  considerations  on  the  ground.  The  six  phase  operation  is 
planned  for  at  least  180  days. 


UNCLASSIFIED  MODELING  SCENARIO 

Scenario  Operational  Timeline 

CJTF  Relief  Phases  of  Operation 
The  AFIT  of  Today  is  ihe  Air  Force  of  Tomorrow. 


Phase  0- 

Shape  political 
and  social 

conditions  in 
Nigeria 

Phase  1  - 

Deter  further 
dispersion  and 
unrest  with  civilian 
population. 

Phase  II  - 

Seize  the  initiative 
through  rapid 
deployment  of 
forces. 

Phase  III  - 

Dominate  the 
countryside  to 
provide  aid  and 
maintain  stability. 

Phase  IV- 

Stabilize  unsettled 
regions  along  with 
host  nation 

government 
forces  and 
agencies. 

Phase  V  - 

Enable  civil 
authority  to 
operate 
independent  of 
coalition  forces. 

now  to  E-Day 

E+lto  D-l 

D-Day  to  D+60 

D+61  to  D+180 

D+181  to  D+365 

D+  180  -  UTC 

Alia  rouse  INSTITUTE  OT  TtGHNOLOST 


Phase  0:  Strategic  and  political  posturing  by  senior  leaders. 

Phase  I:  Military  forces  are  alerted  for  pending  deployment.  Mission  orders  are 
issued  and  deliberate  planning  begins. 

Phase  II:  Military  forces  deploy  and  develop  the  infrastructure  for  follow  on  forces. 
Coordinate  with  coalition  partners  and  host  nation  for  support.  Protect  the  force. 
Phase  III:  Conduct  Humanitarian  operation.  Protect  the  force. 

Phase  IV:  Enable  host  nation  forces  to  conduct  humanitarian  operations 
Phase  V:  Oversight  of  humanitarian  operations  and  re-deploy. 


10  September 2014 
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The  map  below  provides  an  overview  of  the  nine  BCTs  used  to  conduct  the  humani¬ 
tarian  mission.  It  also  shows  the  location  for  possible  A/SPODs. 


UNCLASSIFIED  MODELING  SCENARIO 

W  CJTF-Relief  CONOPS  4SFIT 

Am  Fo acc  Institute  or  rcaHNOuior 


The  AFIT  of  Today  is  the  Air  Force  of  Tomorrow. 


m  Possible  SPODs 

1.  Lagos 

2.  Port  Harcourt 

|  |  Possible  APODs 

1.  Abuja 

2.  Enugu 

3.  Kano 

4.  Lagos 

5.  Port  Harcourt 

^1  FOBS  (~BCT-sized) 

1 .  Abuja  (JTF  HQ  +  Fires  BDE  HQ) 

2.  Oshogbo  (MEB) 

3.  Makurdi  (BFSB) 

4.  Gombe  (IBCT/SBCT) 

5.  Maiduguri  (IBCT/SBCT) 

6.  Bauchi  (CAB) 

7.  Kano  (IBCT/SBCT) 

8.  Bimin  Kebbi  (IBCT/SBCT) 

9.  Kaduna  (CAB) 

FOBs  to  be  located  (not  shown) 

1.  SPOD  (TSC  &  ESC  elements) 

2.  APOD  (TSC  &  ESC  elements) 

3.  SustBde#1 

4.  Sust  Bde  #2 


Air  University:  The  Intellectual  and  Leadership  Center  of  the  Air  Force 
Aim  Hiah...Flv  -  Fiaht  -  Win 
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Appendix  B.  VBA  Code 


The  following  screen  shots  of  the  VBA  code  are  provided  to  document  the  code  used 
to  run  the  model  developed  in  this  research.  They  are  not  intended  for  the  reader 
to  use  to  reproduce  the  code  as  a  copy  of  the  code  was  provided  to  the  sponsor;  the 
copy  in  this  annex  is  provided  as  a  back-up. 

The  format  for  displaying  the  code  is  based  on  the  five  modules:  master,  node  data, 
arc_data,  constraint -data,  and  cplexvba.  Each  screen  shot  is  labeled  at  the  top  with 
the  appropriate  module  name.  Additionally,  each  module  name  is  highlighted  in  blue 
in  the  ”  Project- VBAProject”  window  in  the  upper  left  side  in  each  of  the  screen 
shots. 
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master  (1  of  1) 
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^  F  111  J 

Project  -  VBAProject 

1  3  Q 


£s  J 


m  a 

X 


10  J  ^ 

4=1  ^  =, 
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0  Jgl  VBAProject  (20150202_input  data 

E  •i(§3  Microsoft  Excel  Objects 
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Sheet  10  (weighting) 
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PQ  Sheet  13  (Cl) 

Bj  Sheet  14  (Dl) 

Bj  SheetlS  (El) 

Bj  Sheet  16  (FI) 

®  Sheet  17  (Gl) 

Sheet2  (A) 

Bj  Sheet3  (B) 

®  Sheet4(F) 

B]  Sheets  (C) 

Bj  Sheet6  (D) 

Bj  Sheet7  (E) 

Bj  Sheet8  (G) 

Sheets  (flow) 

ThisWorkbook 
S' Modules 
j  arc_data 

;■  constraint_data 

i  cplexvba 
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^  node_data 


(General) 
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< 


j 


Properties  -  master 
|  master  Module 

Alphabetic  I  Categorized  | 


X J 
"3 


master 


Public  numNodes  As  Integer 
Public  lastRow  As  Long 
Public  lastRow2  As  Long 
Public  lastCol  As  Long 
Public  lastCol2  As  Long 
Public  firstRow  As  Long 
Public  numArcs  As  Integer 
Public  dij Range  As  Range 
Public  cijlRange  As  Range 
Public  xijlRange  As  Range 
Public  uij Range  As  Range 
Public  cij2Range  As  Range 
Public  xij2Range  As  Range 
Public  const47  As  Integer 
Public  const48  As  Integer 
Public  const49  As  Integer 
Public  const50  As  Integer 


Sub  master  () 

'  Stops  the  screen  from  updating  i.e.  will  speed  up  code. 
’Application. ScreenUpdating  =  False 

•  Clears  the  worksheet  below  Row  11 

Rows ("11: 65536") .Delete 

•  Clears  the  OF,  DV,  and  constraints  out  of  CPLEX 

CPXclear 

'  Centers  the  entire  worksheet 
Cells . Select 
With  Selection 

. HorizontalAlignment  =  xlCenter 
•VerticalAlignment  =  xlBottom 
•Wraplext  =  False 
•Orientation  =  0 
•Addlndent  =  False 
. IndentLevel  =  0 
•ShrinkToFit  =  False 
. ReadingOrder  =  xlContext 
•MergeCells  =  False 
End  With 

Call  nodeData 

Call  arcData 

Call  constraints 

MsgBox  "Data  is  now  uploaded." 

CPXsolve 

’  Allows  the  screen  to  update . 

Application. ScreenUpdating  =  True 

End  Sub 
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node.data  (1  of  5) 
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S]  Sheetl6  (FI) 
g)  Sheetl7(Gl) 
g)  Sheet2  (A) 

■)Sheet3(B) 

I  @Sheet4(F) 

a  Sheets  (C) 

®  Sheett  (D) 

@  Sheet7  (E) 

®  Sheet8  (G) 

H]  Sheet9  (flow) 

30  ThisWorkbook 
Modules 
:  arc_data 

constraint_data 
Mt  cplexvba 
««  master 

«gg  I 


j 


Properties  -  node_data 
|node data  Module 

Alphabetic  I  Categorized  | 


node_data 


*J 


Option  Explicit: 


Sub  nodeData() 

Dim  myFile  As  String 

Dim  nodeData  As  String 

Dim  lastComma  As  Integer 

Dim  i  As  Integer 

Dim  j  As  Integer 

Dim  q  As  Integer 

Dim  row_number  As  Integer 

Dim  Jc  As  Variant 

Dim  eachitem  As  Variant 

Dim  nodeCount  As  Integer 

'  Makes  the  gap  between  'Individual  Values'  and  'node-specific  data'  YELLOW. 
Rows ( "11 : 11") .Select 

With  Selection. Interior 
. Pattern  =  xlSolid 
. PatternColorlndex  =  xlAutomatic 
.Color  =  65535 
.TintAndShade  =  0 
. PatternTintAndShade  =  0 
End  With 

'  Makes  the  'Node-specific  data'  row  RED. 

Rows ("12:12") .Select 

With  Selection. Interior 
.Pattern  =  xlSolid 
. PatternColorlndex  =  xlAutomatic 
.Color  =  255 
.TintAndShade  =  0 
.PatternTintAndShade  =  0 
End  With 

'  Writes  the  column  heading  "Node-Specific  Data" . 

Sheets ("Sheetl") .Select 
Range ("bl2") .Select 

ActiveCell . FormulaRICl  =  "Node-Specific  Data" 


"Node-Specific  Data" 


Writes  the  1st  column  heading 

Sheets ("Sheetl") .Select 
Range ("A13") .Select 
ActiveCell . FormulaRICl  =  "Node" 

Range ("bl3") .Select 
ActiveCell . FormulaRICl  =  "b_i" 

Range ("cl3") .Select 

ActiveCell . FormulaRICl  =  "U_iAa" 

Range ("dl3") .Select 

ActiveCell . FormulaRICl  =  "D_iAs" 

Range ("el3") .Select 

ActiveCell . FormulaRICl  =  "U_i^dc" 

Shades  the  column  heading  for  "Node_Specific  Data"  of  first  matrix 
Range ("al3:el3") .Select 
With  Selection. Interior 
.Pattern  =  xlSolid 
.PatternColorlndex  =  xlAutomatic 
.ThemeColor  =  xlThemeColorDarkl 
.TintAndShade  =  -0.349986266670736 
.PatternTintAndShade  =  0 
End  With 
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3  Microsoft  Excel  Objects 
m  Sheet  1  (Sheet  1) 

«  Sheet  10  (weighting) 
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m  Sheet  13  (Cl) 

<5]  Sheetl4(Dl) 
m  Sheet  15  (El) 
m  Sheet  16  (FI) 

^  Sheet  17  (Gl) 
a  Sheet2  (A) 
m  Sheet3 (B) 

I  B  Sheet4(F) 

a  Sheets  (C) 
a  Sheett  (D) 

@  Sheet7 (E) 

(P)  Sheet8  (G) 

@j  Sheet9  (flow) 
jT)  This  Workbook 
“3  Modules 

arc_data 
constraint_data 
cplexvba 
master 

«%l  i 


Properties  -  node_data 


|node data  Module 

Alphabetic  I  Categorized  ] 


*1 

U 


node  data 


|  (General) 


'  Writes  the  1st  column  heading  "Node-Specific  Data". 

Sheets ("Sheetl") .Select 

Range ("gl3") .Select 

ActiveCell . FormulaRlCl  =  "y_iAa" 

Range ("hl3") .Select 

ActiveCell . FormulaRlCl  =  "y_i^s" 

Range ("il3") .Select 

ActiveCell . FormulaRlCl  =  "y_i~dc" 

Range ("jl3") .Select 

ActiveCell . FormulaRlCl  =  "z_i~c" 

Range ("kl3") .Select 

ActiveCell . FormulaRlCl  =  "z_iAnc" 

Shades  the  column  heading  for  the  decision  variable  (second)  matrix. 

Range  ("gl3 : )cl3")  .Select 
With  Selection. Interior 
.Pattern  =  xlSolid 
. PatternColorlndex  =  xlAutomatic 
.ThemeColor  =  xlThemeColorDarlcl 
. TintAndShade  =  -0.349986266670736 
. PatternTintAndShade  =  0 
End  With 

Places  title  "Decision  Variable  Matrix"  above  the  decision  variable  matrix. 
Range ("il2") .Select 

ActiveCell . FormulaRlCl  =  "Decision  Variable  Matrix" 


1  Reads  the  node  data  from  a  .txt  file 

'myFile  =  "I:\My  Documents\Thesis\Chapter  4\VBA  Interface\node_datal.txt" 
myFile  =  "I:\My  Documents\Thesis\Chapter  4\VBA  Interface\node_data2.txt" 
Open  myFile  For  Input  As  #1 


'  Converts  the  data  in  the 
nodeCount  =  0 
Cells (1,  1). Select 
row_n umber  =  13 
Do  Until  EOF (1) 

Line  Input  #1,  Jc 
eachitem  =  Split (k. 


.txt  file  into  Excel-cell  format 


ActiveCell .Off set (row_number, 
ActiveCell .Off set (row_number, 
ActiveCell .Offset (row_number, 
ActiveCell .Offset (row_number, 
ActiveCell .Offset (row_number, 
row_n umber  =  row_number  +  1 
nodeCount  =  nodeCount  +  1 


0) .Value 

1)  .Value 

2)  .Value 

3)  .Value 

4)  .Value 


eachitem(O) 
eachitem (1) 
eachitem (2 ) 
eachitem (3) 
eachitem (4) 


Loop 


’  Finds  the  position  of  the  last  row 

lastRow  =  ActiveSheet .Cells (Rows .Count,  "a") .End(xlUp) .Row 
Close  #1 

'  Assigns  the  number  of  nodes  to  ' numNodes ' 
numNodes  =  nodeCount 
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node.data  (3  of  5) 


File  Edit  View  Insert  Format  Debug  Run  Jools  Add-Ins  Window  Help 

li'd  <  ►  u  a  K  $ 


!  .-=  1  lig: 


Ml  B 

a  3  a _ 

ASAPUtilitiesRibbon  (ASAP_Utilrtie~ 

B  ^  VBAProject  (20150202 Jnput  data 

B--<§3  Microsoft  Excel  Objects 
i  ®  Sheetl  (Sheetl) 

!  IfQ  SheetlO  (weighting) 

(P|  Sheetll(Al) 

@  Sheet  12  (Bl) 

(P|  Sheet  13  (Cl) 

@  Sheetl4(Dl) 

Sheetl5  (El) 

]  Sheet  16  (FI) 

}  Sheet  17  (Gl) 

]  Sheet2  (A) 
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@Sheet7(E) 
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:  cplexvba 
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Properties  -  node_data 


node_data  Module 
Alphabetic  |  Categorized  | 
node_data 


XJ 


▼  |  JnodeData 


Makes  a  grid  system  for  the  'Node-Specific  Data'  matrix 
Range ("A13:E"  £  lastRow) .Select 
With  Selection. Borders 

.LineStyle  =  xlContinuous 
.Colorlndex  =  0 
.TintAndShade  =  0 
.Weight  =  xlThin 
End  With 

Makes  a  grid  system  for  the  'Decision  Variable  Matrix' 

Range ("gl3 : k"  £  lastRow) . Select 
With  Selection. Borders 

•LineStyle  =  xlContinuous 
.Colorlndex  =  0 
.TintAndShade  =  0 
.Weight  =  xlThin 
End  With 

Places  "0"  in  all  of  the  cells  of  the  "Decision  Variable  Matrix' 
Range (Cells (14,  7),  Cells (14  +  numNodes  -  1,  11))  =0 


of  first  matrix 


Shades  the  "Node"  row  for  "Node-Specific  Data" 

Range ("al3:a"  £  lastRow) .Select 
With  Selection. Interior 
.Pattern  =  xlSolid 
. PatternColorlndex  =  xlAutomatic 
.ThemeColor  =  xlThemeColorDarkl 
.TintAndShade  =  -0.349986266670736 
. PatternTintAndShade  =  0 
End  With 

Shades  the  APOD  decision  variables  (y_i~a)  blue  that  are  possible  solutions  in  the  DV  matrix 
For  i  =  lastRow  -  numNodes  +  1  To  lastRow 
If  Cells (i,  3)  =  0  Then 
Else 


CPXaddVariable  Variable :=Cells (i,  7),  Lb:=0,  Binary:=True  'Add  CPLEX  decision  variables  everytime  cell  turns  blue 
Cells (i,  7). Select 
With  Selection. Interior 
.Pattern  =  xlSolid 
. PatternColorlndex  =  xlAutomatic 
.Color  =  15773696 
.TintAndShade  =  0 
.PatternTintAndShade  =  0 
End  With 
End  If 
Next  i 
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node.data  (4  of  5) 
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node  Data 


Shades  the  APOD  decision  variables  (y_iAa)  blue  chat  are  possible  solutions  in  the  DV  matrix 
For  i  =  lastRow  -  numNodes  +  1  To  lastRow 
If  Cells (i,  3)  =  0  Then 
Else 

CPXaddVariable  Variable :=Cells (i,  7),  Lb:=0,  Binary:=True  ‘Add  CPLEX  decision  variables  everytime  cell  turns  blue 
Cells(i,  7) .Select 
With  Selection. Interior 
. Pattern  =  xlSolid 
. PatternColorlndex  =  xlAutomatic 
.Color  =  15773696 
.TintAndShade  =  0 
. PatternTintAndShade  =  0 
End  With 
End  If 
Next  i 

Shades  the  SPOD  decision  variables  (y_iAs)  blue  that  are  possible  solutions  in  the  DV  matrix 
For  i  =  lastRow  -  numNodes  +  1  To  lastRow 
If  Cells (i,  4)  =  0  Then 
Else 

Cells (i,  8). Select 

CPXaddVariable  Variable :=Cells (i,  8),  Lb:=0,  Binary:=True  ‘Add  CPLEX  decision  variables  everytime  cell  turns  blue 
With  Selection. Interior 
. Pattern  =  xlSolid 
.PatternColorlndex  =  xlAutomatic 
•Color  =  15773696 
.TintAndShade  =  0 
.PatternTintAndShade  =  0 
End  With 
End  If 
Next  i 


Shades  the  DC  decision  variables  (y_i*dc)  blue  that  are  possible  solutions  in  the  DV  matrix 
For  i  =  lastRow  -  numNodes  +  1  To  lastRow 
If  Cells (i,  5)  =  0  Then 
Else 

CPXaddVariable  Variable :=Cells (i,  9),  Lb:=0,  Binary:=True  ’Add  CPLEX  decision  variables  everytime  cell  turns  blue 
Cells (i,  9). Select 

With  Selection. Interior 
.Pattern  =  xlSolid 
•PatternColorlndex  =  xlAutomatic 
.Color  =  15773696 
.TintAndShade  =  0 
.PatternTintAndShade  =  0 
End  With 
End  If 
Next  i 
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g  Sheet4  (F) 
g  Sheet5  (C) 
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g  ThisWorkbook 
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Properties  -  node_data 
|  node  data  Module 


x| 
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j  (General) 


With  Selection. Interior 
. Pattern  =  xlSolid 
. PatternColorlndex  =  xlAutomatic 
.Color  =  15773696 
.TintAndShade  =  0 
. PatternTintAndShade  =  0 
End  With 
End  If 
Next  i 

Shades  the  DC  decision  variables  (y_i*dc)  blue  that  are  possible  solutions  in  the  DV  matrix 
For  i  =  lastRow  -  numNodes  +  1  To  lastRow 
If  Cells (i,  5)  =  0  Then 
Else 

CPXaddVariable  Variable :=Cells (i,  9),  Lb:=0,  Binary:=True  'Add  CPLEX  decision  variables  everytime  cell  turns  blue 
Cells (i,  9). Select 

With  Selection. Interior 
.Pattern  =  xlSolid 
-PatternColorlndex  =  xlAutomatic 
.Color  =  15773696 
•TintAndShade  =  0 
.PatternTintAndShade  =  0 
End  With 
End  If 


This  section  shades  DC  decision  variables  (y_i"dc)  blue  that  are  possible  solutions  in  the  DV  matrix  _ 
and  reads  them  into  CPLEX  as  DVs. 

Dim  ydv  As  Range 

Set  ydv  =  Range (Cells (14,  9),  Cells (14  +  numNode3  -  1,  9)) 

CPXaddVariable  Variable :=ydv,  Lb:=0,  Binary :=True 
ydv. Select 

With  Selection. Interior 

.Pattern  =  xlSolid 
.PatternColorlndex  =  xlAutomatic 
.Color  =  15773696 
.TintAndShade  =  0 
.PatternTintAndShade  =  0 
End  With 

This  section  shades  DVs  (z_i~c,  z_i~nc)  blue  if  they  return  a  non-zero  value  and  reads  them  into  CPLEX  as  DVs. 
Dim  zdv  As  Range 

Set  zdv  =  Range ("j 14 :K"  &  lastRow) 

CPXaddVariable  Variable :=zdv,  Lb:=0 

Range ("j 14 :K"  &  lastRow) .Select 
With  Selection. Interior 
.Pattern  =  xlSolid 
. PatternColorlndex  =  xlAutomatic 
.Color  =  15773696 
.TintAndShade  =  0 
.PatternTintAndShade  =  0 
End  With 
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■3  Sheet  16  (PI) 

Die  f  As  integer 

C  Sheet  17  (Gl) 

Die  1  As  Integer 

13  Sheet2  (A) 

Die  j  As  Integer 

(]Sheet3(B) 
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Die  Jc  As  Integer 

■3  Sheet.  (F) 

Die  p  As  Integer 

13  Sheets  (C) 

Die  t  As  Integer 

■3  Sheets  (D) 

Die  v  As  Integer 

•]Sheet7(E) 

Die  z  As  Integer 

1]  Sheets  (G) 
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■T)  Sheet9  (flow) 

Die  1  As  Variant 

%J  ThsWorfcbook 

Die  eachitem  As  Variant 

S  £3  Modules 

Die  startPos  As  Integer 

~i 

coastrant_data 

oplexvba 

'Die  firstRov  As  Integer 

master 

firstRov  -  lascRow  +  3  '  This  provides  the  position  for  the  grid. 

node_data 

- 

«  m  , 

startPos  ■  nuEilodes  ♦  14 

•  Makes  the  gap  between  'Node-specific  data*  and  'Arc-specific  data'  YELLOW. 
Rows (startPos) .Select 

Properties  -  arc.data 

*J 

arc  data  MocUe  -1 

with  Selection. Interior 

Alphabetic  |  Categorized  | 

.Pattern  ■  xlSolid 
.  Patterr.Color Index  ■  xlAutomatic 

I  arc.data 

.Color  -  6S53S 

. Tin t AndShade  ■  0 

•PatternTintAndShade  -  0 

End  With 

'  Makes  the  'Arc-specific  data'  row  RED. 

Rows (startPos  ♦  1). Select 

With  Selection. Interior 

.Pattern  -  xlSolid 
.  Patterr.Color  Index  -  xlAutomatic 
.Color  -  2SS 
.Tint AndShade  ■  0 

.PatternTintAndShade  -  0 

End  With 

'  Writes  the  column  heading  "Arc-Specific  Data". 

Cells (startPos  ♦  1,  2). Select 

AetlveCell . ForuulaRICl  -  "Arc-Specific  Data" 
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constraint_data 
|  cplexvba 

;••••<$£  master 
:  node_data 

'  1 


®S 

0T 


(General) 


'  Writes  the  1st  column  heading  "Arc-Specific  Data" . 

Cells (startPos  +  2,  1). Select 
ActiveCell . FormulaRICl  =  "From" 

Cells  (startPos  +  2,  2). Select 
ActiveCell . FormulaRICl  =  "To" 

Cells (startPos  +  2,  3) .Select 
ActiveCell . FormulaRICl  =  "d_ij" 

Cells  (startPos  +  2,  4). Select 
ActiveCell . FormulaRICl  =  "c_ijAl" 

Cells  (startPos  +  2,  5). Select 
ActiveCell . FormulaRICl  =  "c_ij~2" 

Cells (startPos  +  2,  6). Select 
ActiveCell . FormulaRICl  =  "u_ij" 

•  Shades  the  first  row  of  the  "Arc_Specific  Data"  matrix 
Range (Cells (firstRow,  1),  Cells (firstRow,  6)). Select 
With  Selection. Interior 
.Pattern  =  xlSolid 
. PatternColorlndex  =  xlAutomatic 
.ThemeColor  =  xlThemeColorDarJcl 
.TintAndShade  =  -0.349986266670736 
. PatternTintAndShade  =  0 
End  With 

■  Reads  the  arc  data  from  a  .txt  file 

'myFile  =  "I:\My  Documents\Thesis\Chapter  4\VBA  Interface\arc_datal.txt" 

myFile  =  "I:\My  Documents\Thesis\Chapter  4\VBA  Interface\arc_data2.txt" 

Open  myFile  For  Input  As  #1 


'  Put  the  arc  data  into  appropriate  format  for  use  in  CPLEX 
Cells  (1,  1). Select 
row_number  =  startPos  +  2 
numArcs  =  0 
Do  Until  EOF (1) 

Line  Input  #1,  1 
eachitem  =  Split  (1,  ",") 


ActiveCell .Offset (row_number, 
ActiveCell .Offset (row_number, 
ActiveCell .Offset (row_number, 
ActiveCell . Offset ( row_number , 
ActiveCell . Of f set ( row_number , 
ActiveCell .Offset (row_number, 
row_number  =  row_n umber  +  1 
numArcs  =  numArcs  +  1 


0) .Value 

1)  .Value 

2)  .Value 

3)  .Value 

4)  .Value 

5)  .Value 


eachitem (0) 
eachitem (1) 
eachitem (2) 
eachitem  (3) 
eachitem  (4) 
eachitem (5) 


Loop 

'  P  is  the  starting  point  for  w;  used  to  read  arc  data  into  d_ij  matrix, 
p  =  lastRow  +  4 

********  Finds  the  position  of  the  last  row***"*** 
lastRow  =  ActiveSheet. Cells (Rows. Count,  "a") .End(xlUp) .Row 
'  MsgBox  lastRow 

Close  *1 
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►  ii  *  „  c*  *!  in  m 

6* 

SJ  -  .  A  -  -  . 

Project  -  VBAProject  X 

|  (General) 

EE)'"4$  ASAPUtilities  Ribbon  (ASAP Utilrtie 

A 

Range (Cells (firstRow,  1),  Cells (lastRow,  6)). Select 

E  ^  VBAProject  (20150202_input  data 

With  Selection. Borders 

IE1 -{Sj  Microsoft  Excel  Objects 

.LineStyle  =  xlContinuous 

Sheetl  (Sheetl) 

•Colorlndex  =  0 

if)  Sheet  10  (weighting) 

.TintAndShade  =  0 

if)  Sheet  11  (Al) 

.Weight  =  xlThin 

■3  Sheet  12  (Bl) 

End  With 

@  Sheet  13  (Cl) 

If]  Sheet  14  (Dl) 

'  Shades  the  first  column  of  the  "Arc  Specific  Data"  matrix. 

Sheetl5  (El) 

Range (Cells (firstRow,  1),  Cells ( (firstRow  +  numArcs),  2)). Select 

Hj  Sheet  16  (FI) 

With  Selection. Interior 

B]  Sheet  17  (Gl) 

.Pattern  =  xlSolid 

SQ  Sheet2  (A) 

. PatternColorlndex  =  xlAutomatic 

Sheet3  (B) 

.ThemeColor  =  xlThemeColorDarkl 

B]  Sheet4  (F) 

.TintAndShade  =  -0.349986266670736 

@Sheet5(C) 

. PatternTintAndShade  =  0 

■f)  Sheet6  (D) 

End  With 

SP)Sheet7(E) 

:  S]Sheet8(G) 

■  Writes  the  column  heading  for  the  x  ij~l  and  x  ijA2  decision  variable  matrix. 

B}  Sheets  (flow) 

Cells (startPos  +  2,  8). Select 

L  ThisWorkbook 

ActiveCell .  FormulaRICl  =  "x  ij/vl" 

Modules 

Cells (startPos  +  2,  9). Select 

«l - 1 

ActiveCell .  FormulaRICl  =  "x  ij~2" 

«$£  constraint_data 

cplexvba 

'  Shades  the  heading  for  the  x  ij^l  and  x  ij~2  decision  variable  matrix. 

1  master 

Range (Cells (startPos  +  2,  8),  Cells (startPos  +  2,  9)). Select 

node_data 

▼ 

With  Selection. Interior 

.Pattern  =  xlSolid 

Properties  -  arc_data 

x 

.ThemeColor  =  xlThemeColorDarkl 

arc  data  Module 

.TintAndShade  =  -0.349986266670736 

Alphabetic  |  Categorized  | 

End  With 

SSjarcjdat^^^ 

Makes  the  grid  system  for  the  x  ij^l  and  x  ijA2  decision  variable  matrix. 
Range (Cells (startPos  +2,  8),  Cells (lastRow,  9)). Select 
With  Selection. Borders 

.LineStyle  =  xlContinuous 
.Colorlndex  =  0 
.TintAndShade  =  0 
.Weight  =  xlThin 
End  With 


d_ij,  u_ij,  c_ij*l,  c_ij~2,  x_ij~l,  and  x_ij~2  matrix 


firstRow  =  lastRow  +  4 

i  i  i  i  i  ■  i  i  i  i  i  i  . . .  i  i  i  i  i  i  i  i  i  i  i  i  i  Begin  d_i j  matrix 

•  Places  the  d_ij  heading 

Cells (firstRow,  1) .Select 
ActiveCell . FormulaRICl  =  "d_ij" 
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!  ®  Sheets  (C) 

|~  ®  Sheet6  (D) 

J  ®  Sheet7  (E) 

®  Sheets  (G) 

Sheet9  (flow) 

ThisWorkbook 
£-S  Modules 

«g!°!  i 

I  constraint_data 

!■•••«$£  cplexvba 
|  master 

node_data 

G 


Properties  -  arc_data 


|arc_data  Module 


* 


*1 

“II 


|  (General) 


Places  the  row  designation  for  d_ij  matrix 
j  =  1 

For  i  =  (firstRow  +  1)  To  (firstRow  +  numNodes) 

Cells (i,  2)  =  j 
3=3+1 
Next  i 

Places  the  column  designation  for  d_ij 
k  =  1 

For  i  =  3  To  numNodes  +  2 
Cells (firstRow,  i)  =  k 
k  =  k  +  1 
Next  i 

Makes  a  grid  system  for  the  d_ij  matrix 

Range (Cells (firstRow,  2),  Cells (firstRow  +  numNodes,  numNodes  +  2)). Select 
With  Selection. Borders 

.LineStyle  =  xlContinuous 
.Colorlndex  =  0 
.TintAndShade  =  0 
.Weight  =  xlThin 
End  With 

Shades  the  first  row  of  the  "d_ij"  matrix 

Range (Cells (lastRow  +  4,  2),  Cells (lastRow  +  4,  numNodes  +  2)). Select 
With  Selection. Interior 
.Pattern  =  xlSolid 
. PatternColorlndex  =  xlAutomatic 
.ThemeColor  =  xlThemeColorDarkl 
.TintAndShade  =  -0.349986266670736 
. PatternTintAndShade  =  0 
End  With 


numNodes,  2)). Select 


Shades  the  first  column  of  the  "d_ij"  matrix 
Range (Cells (lastRow  +4,  2),  Cells (lastRow 
With  Selection. Interior 
.Pattern  =  xlSolid 
. PatternColorlndex  =  xlAutomatic 
.ThemeColor  =  xlThemeColorDarkl 
.TintAndShade  =  -0.349986266670736 
.PatternTintAndShade  =  0 
End  With 


'  Places  "0"  in  all  of  the  cells  of  d_ij 

Range (Cells (firstRow  +1,  3),  Cells (firstRow  +  numNodes,  numNodes 


2))  =  0 


Shades  d_ij  brown  if  a  non-zero  value 

Range (Cells (firstRow  +1,  3),  Cells (firstRow  +  numNodes,  numNodes  +  2)). Select 
Selection. FormatConditions .Add  Type:=xlCellValue,  Operator :=xlNotEqual,  _ 
Formulal :="=0" 

Selection. FormatConditions (Selection. FormatConditions . Count) . SetFirstPriority 
With  Selection. FormatConditions (1) .Interior 
.PatternColorlndex  =  xlAutomatic 
.Color  =  3368601 
.TintAndShade  =  0 
End  With 

Selection . FormatConditions (1) . StopIfTrue  =  False 
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Q  ThisWorkbook 

^3  Modules 

M  1 

!  constraint_data 

*«£  cplexvba 

|  master 

:  node_data 

- 

Properties  -  arc_data 

X 

arc  data  Module 

- 3 

Alphabetic  |  Categorized  | 

arc_data 

|  (General) 


▼  [arcData 


.  the  appropriate  cells 


'  Places  the  d_ij '  s  ii 
For  w  =  p  To  lastRow 
f  =  Cells (w,  1) 
t  =  Cells (w,  2) 
dij  =  Cells  (w,  3) 

For  i  =  (firstRow  +  1)  To  (numNodes  +  firstRow)  ’  This  checks  all  the  row  values. 
For  j  =  3  To  (numNodes  +2)  '  This  checks  all  the  column  values. 

If  (Cells  (i,  2)  =  f  And  Cells (firstRow,  j)  =  t)  Then 
Cells (i,  j)  =  dij 

End  If 
Next  j 


1  Creates  d_ij  matrix  as  a  range  for  use  in  constraint  calculations 

Set  dijRange  =  Sheets ("Sheetl") .Range (Cells (firstRow  +1,  3),  Cells (firstRow  +  numNodes,  numNodes  • 


End  d_ij  matrix 


'  ’  ' ’ ' ' ’ ' 1 ' 1 ’ ’ ' ' 1 ' ' ’ ’ ' 1 ’ ’ ’ ’ ' ' ’ ' 1 ’ 1 'Begin  c_ij~l  matrix  . . .  '  '  '  '  '  '  '  '  '  '  '  '  '  '  '  '  '  '  '  '  '  ' 

'  Finds  the  position  of  the  last  col 

Finds  the  position  of  the  last  column 

lastRow  =  ActiveSheet. Cells (Rows. Count,  "b") .End(xlUp) .Row 
lastCol  =  Cells (lastRow,  Columns. Count) .End (xlToLeft) .Column 
Places  the  d_ij  heading 

Cells (firstRow,  lastCol  +  2) .Select 
ActiveCell . FormulaRICl  =  "c_ij~l" 

Places  the  row  designation  for  c_ij'‘l  matrix 
3=1 

For  i  =  (firstRow  +1)  To  (firstRow  +  numNodes) 

Cells  (i,  lastCol  +  3)  =  j 
j  =  3  +  1 
Next  i 

Places  the  column  designation  for  c_ij  ~1  matrix 

k  =  1 

For  i  =  (lastCol  +  4)  To  lastCol  +  3  +  numNodes 
Cells (firstRow,  i)  =  k 
k  =  k  +  1 
Next  i 

Makes  a  grid  system  for  the  c_ij~l  matrix 

Range (Cells (firstRow,  lastCol  +  3),  Cells (firstRow  +  numNodes,  lastCol  +  3 
With  Selection. Borders 

.LineStyle  =  xlContinuous 
.Colorlndex  =  0 
. TintAndShade  =  0 
.Weight  =  xlThin 
End  With 


numNodes) ) .Select 
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®Sbeet4(F) 

®Sheet5(C) 

®  Sheets  (D) 

®Sheet7(E) 

|  ®  Sheet8  (G) 

®  Sheet9  (flow) 

$5*1  ThisWorkbook 
B  ■  ®  Modules 

arc_data 

i  constramt_data 

j  tflj  cplexvba 
|  master 

node_data 


Properties  -  arc_data 


arc_data  Module 
Alphabetic  |  Categorized  | 


arc Data 


Shades  the  first  row  of  the  matrix 

Range (Cells (lastRow  -  numNodes,  lastCol  +  3),  Cells (lastRow  -  numNodes,  lastCol  +  3  +  numNodes) ) .Select 
With  Selection. Interior 
.Pattern  =  xlSolid 
. PatternColorlndex  =  xlAutomatic 
.ThemeColor  =  xlThemeColorDarlcl 
.TintAndShade  -  -0.349986266670736 
.PatternTintAndShade  =  0 
End  With 

Shades  the  first  column  of  the  "c_i;j''l"  matrix 

Range (Cells (lastRow  -  numNodes,  lastCol  +  3),  Cells (lastRow,  lastCol  +  3)). Select 
With  Selection. Interior 
.Pattern  =  xlSolid 
. PatternColorlndex  -  xlAutomatic 
.ThemeColor  =  xlThemeColorDarlcl 
.TintAndShade  =  -0.349986266670736 
.PatternTintAndShade  =  0 
End  With 


'  Places  "0"  in  all  of  the  cells  of 
Range (Cells (fir stRow  +  1,  lastCol  + 


:_ij“l  matrix 

4),  Cells (fir stRow  h 


numNodes,  lastCol  +  3  +  numNodes))  : 


Shades  c_ij'‘l  brown  if  a  non-zero  value 

Range (Cells (firstRow  +  1,  lastCol  +4),  Cells (firstRow  +  numNodes,  lastCol  +  3  +  numNodes) ) .Select 
Selection. FormatConditions . Add  Type:=xlCellValue,  Operator :=xlNotEqual,  _ 

Formulal 

Selection . FormatConditions (Selection . FormatConditions . Count) . SetFirstPriority 
With  Selection. FormatConditions (1) .Interior 
. PatternColorlndex  =  xlAutomatic 
.Color  =  3368601 
.TintAndShade  =  0 
End  With 

Selection. FormatConditions (1) .Stoplf True  =  False 


Ln  the  appropriate  cells| 
-  numNodes  -  4) 


*  Places  the  c_ij "1 1 a  : 

For  w  =  p  To  (lastRow  - 
f  =  Cells (w,  1) 
t  =  Cells (w,  2) 
cijl  =  Cells (w,  4) 

For  i  =  (firstRow  +  1)  To  (numNodes  +  firstRow)  1  This  checks  all  the  row  values. 

For  j  =  (lastCol  +  4)  To  (lastCol  +  numNodes  +  4)  ■  This  checks  all  the  column  values. 

If  (Cells (i,  2)  =  f  And  Cells (firstRow,  j)  =  t)  Then 
Cells (i,  j)  =  ci31 

End  If 
Next  3 


•  Creates  c  ij'l  matrix  as  a  range  for  use  in  constraint  calculations 

Set  ci31Range  =  Sheets ("Sheetl”) .Range (Cells (firstRow  +  1,  lastCol  +  4),  Cells (firstRow  +  numNodes,  lastCol  - 


'End 


i3"l  matrix 
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©  Sheet8  (G) 
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;  constraint_data 

i  cplexvba 

v*:J  master 
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'  Places  the  u_ij  heading 
firstRow  =  lastRow  +  2 

Cells (firstRow,  1). Select 
ActiveCell . FormulaRICl  =  "u_ij" 

'  Places  the  row  designation  for  u_ij  matrix 

j  =  1 

For  i  =  (firstRow  +  1)  To  (firstRow  +  numNodes) 

Cells (i,  2)  =  j 
j  =  j  +  1 
Next  i 

'  Places  the  column  designation  for  u_ij 
k  =  1 

For  i  =  3  To  numNodes  +  2 
Cells (firstRow,  i)  =  k 
Jc  =  Jc  +  1 
Next  i 

'  Makes  a  grid  system  for  the  u_ij  matrix 

Range (Cells (firstRow,  2),  Cells (firstRow  +  numNodes,  numNodes  +  2)). Select 
With  Selection. Borders 

.LineStyle  =  xlContinuous 
.Color Index  =  0 
.TintAndShade  =  0 
.Weight  =  xlThin 
End  With 

'  Places  "0"  in  all  of  the  cells  of  u_ij 

Range (Cells (firstRow  +1,  3),  Cells (firstRow  +  numNodes,  numNodes  +2))  =0 


'  Shades  c_ij~l  brown  if  a  non-zero  value 

Range (Cells (firstRow  +1,  3),  Cells (firstRow  +  numNodes,  numNodes  +  2)). Select 
Selection. FormatCondit ions .Add  Type :=xlCellValue,  Operator :=xlNotEqual,  _ 
Formulal : ="=0" 

Selection . FormatConditions ( Selection . FormatConditions . Count) . Set First Priority 
With  Selection. FormatConditions (1) .Interior 
. PatternColorlndex  =  xlAutomatic 
.Color  =  3368601 
.TintAndShade  =  0 
End  With 

Selection. FormatConditions (1) . StopIfTrue  =  False 

'  Shades  the  first  row  of  the  "u_ij"  matrix 

Range (Cells (firstRow,  2),  Cells (firstRow,  numNodes  +  2)). Select 
With  Selection. Interior 
.Pattern  =  xlSolid 
. PatternColorlndex  =  xlAutomatic 
.ThemeColor  =  xlThemeColorDarkl 
.TintAndShade  =  -0.349986266670736 
. PatternTintAndShade  =  0 
End  With 

'  Shades  the  first  column  of  the  "u_ij"  matrix 

Range (Cells (firstRow,  2),  Cells (firstRow  +  numNodes,  2)). Select 
With  Selection. Interior 
.Pattern  =  xlSolid 
. PatternColorlndex  =  xlAutomatic 
.ThemeColor  =  xlThemeColorDarkl 
.TintAndShade  =  -0.349986266670736 
. PatternTintAndShade  =  0 
End  With 
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®Sheet7(E) 
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~3  F 


'  Places  the  u_ij 1 s  in  the  appropriate  cells 
For  w  =  p  To  (lastRow  -  numNodes  -  4) 
f  =  Cells (w,  1) 
t  =  Cells (w,  2) 
uij  =  Cells (w,  6) 

For  i  =  (firstRow  +  1)  To  (numNodes  +  firstRow)  •  This  checks  all  the  row  values. 

For  j  =  3  To  (numNodes  +  2)  1  This  checks  all  the  column  values. 

If  (Cells (i,  2)  =  f  And  Cells (firstRow,  j)  =  t)  Then 
Cells (i,  j)  =  uij 

End  If 
Next  j 
Next  i 

Next  w 

1  Creates  u_ij  matrix  as  a  range  for  use  in  constraint  calculations 

Set  uij Range  =  Sheets ("Sheetl") .Range (Cells (firstRow  +1,  3),  Cells (firstRow  +  numNodes,  numNodes  +  2)) 
uij  Range . Select 

. . i  i  i  i  . . .  . . . . End  u_ij  matrix  . . . . . . . .  1  1 


. . . . . . . . . . . . . . . . Begin  c_ij~2  matrix 

Places  the  c_ij_2  heading 

Cells (firstRow,  lastCol  +  2) .Select] 

ActiveCell . FormulaRICl  =  "c_ij“2" 

Places  the  row  designation  for  c_ij"2  matrix 
j  =  1 

For  i  =  (firstRow  +1)  To  (firstRow  +  numNodes) 
Cells (i,  lastCol  +  3)  =  j 
j  =  j  +  1 


Next  i 

'  Places  the  column  designation  for  c_ij~2  matrix 
k  =  1 

For  i  =  (lastCol  +4)  To  lastCol  +  3  +  numNodes 
Cells (firstRow,  i)  =  k 
k  =  k  +  1 
Next  i 

•  Makes  a  grid  system  for  the  c_ij/'2  matrix 

Range (Cells (firstRow,  lastCol  +  3),  Cells (firstRow  +  numNodes,  lastCol  +  3  +  numNodes) ). Select 
With  Selection. Borders 

.LineStyle  =  xlContinuous 
.Colorlndex  =  0 
.TintAndShade  =  0 
.Weight  =  xlThin 
End  With 

’  Places  "0"  in  all  of  the  cells  of  c_ij~2  matrix 

Range (Cells (firstRow  +  1,  lastCol  +  4),  Cells (firstRow  +  numNodes,  lastCol  +  3  +  numNodes))  =  0 
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El  Microsoft  Excel  Objects 
®  Sheet  1  (Sheet  1) 

®  SheetlO  (weighting) 

I  Iff)  Sheet  11  (Al) 

®Sheetl2(Bl) 

®Sheetl3(Cl) 

®Sheetl4(Dl) 

:  ®  SheetlS  (El) 

@Sheetl6(Fl) 
i  ®  Sheetl7  (Gl) 

@Sheet2(A) 

@Sheet3(B) 

I  ®Sheet4(F) 

©  Sheets  (C) 

I  fi  Sheet6  (D) 

©Sheet7(E) 

©  Sheets  (G) 

©  Sheets  (flow) 

®  ThisWorkbook 
El  '3  Modules 


|  (General) 


~3  I" 


Shades  c_i3~2  brown  if  a  non-zero  value 
Range (Cells (firstRow  +  1,  lastCol  +  4),  Cells (firstRow  +  numNodes,  lastCol  +  3 
Selection. Forma ^Conditions .Add  Type:=xlCellValue,  Operator :=xlNotEqual,  _ 
Formulal 

Selection . FormatConditions ( Selection . FormatConditions . Count) . SetFirstPriority 
With  Selection. FormatConditions (1) .Interior 
.PatternColorlndex  =  xlAutomatic 
.Color  =  3368601 
. TintAndShade  =  0 
End  With 

Selection. FormatConditions (1) .StopIfTrue  =  False 


Shades  the  first  row  of  the  "c_ij''2"  matrix 

Range (Cells (firstRow,  lastCol  +  3),  Cells (firstRow,  lastCol  +  3  +  numNodes) ). Select 
With  Selection. Interior 
.Pattern  ”  xlSolid 
.PatternColorlndex  =  xlAutomatic 
.ThemeColor  =  xlThemeColorDarkl 
.TintAndShade  =  -0.349986266670736 
. PatternTintAndShade  =  0 
End  With 


numNodes) ) .Select 


Properties  -  arc_data 


|arc data  Module 


■  Shades  the  first  column  of  the  "c_ij/'2"  matrix 

Range (Cells (firstRow,  lastCol  +  3),  Cells (firstRow  +  numNodes,  lastCol  +  3)). Select 
With  Selection. Interior 
.Pattern  -  xlSolid 
.PatternColorlndex  =  xlAutomatic 
.ThemeColor  =  xlThemeColorDarkl 
. TintAndShade  =  -0.349986266670736 
.PatternTintAndShade  =  0 
End  With 

■  Places  the  c_±l‘'2,s  in  the  appropriate  cells 

For  w  =  p  To  (lastRow  -  numNodes  -  4) 
f  =  Cells (w,  1) 
t  =  Cells (w,  2) 
cij2  =  Cells (w,  5) 

For  i  =  (firstRow  +  1)  To  (numNodes  +  firstRow)  *  This  checks  all  the  row  values. 

For  j  =  (lastCol  +  4)  To  (lastCol  +  numNodes  +  4)  ‘  This  checks  all  the  column  values. 

If  (Cells (i,  2)  =  f  And  Cells (firstRow,  j)  =  t)  Then 
Cells (i,  j)  -  cij2 

End  If 
Next  j 
Next  i 

Next  w 

’  Creates  c_i3”'2  matrix  as  a  range  for  use  in  constraint  calculations 

Set  cij2Range  -  Sheets ("Sheetl”) .Range (Cells (firstRow  +  1,  lastCol  +  4),  Cells (firstRow  +  numNodes,  lastCol  +  3  +  numNodes)) 

cij2Range . Select 


:_i j  ^2  matrix  1 
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End  Mich 

Cells (firscRow  -  4  -  numNodes  -  2  -  numArcs  +  z,  9).Sele 
Mich  Seleccion. Incerior 

•PaccernColorlndex  =  xlAucomacic 
•Color  =  15773696 
•TincAndShade  =  0 

Cells  (i,  j)  =  "="  l  Cells  (firscRow  -  4  -  numNod 
RowAbsoluce : =False ,  ColumnAbsoluce : =False ) 
End  Mich 


'  Colors  che  x_ij"2  cell  in  Che  “Arc-specific  daca“  seccion  blue 


,  9) .Address (ReferenceScyle:=xlAl, 


o  che  blue  cell  in  che  “Arc-specific  Daca  Seccion". 


the  code  below  exisced  before  my  accesses  co  shade  x_ig“2 ' 3  blue  if  che  u_l 
Shades  x_ij"2  blue  if  ic  recurns  a  non-zero  value 
Range (Cells (lascRow  +  3,  lascCol  +  lascCol  +  5),  Cells (lascRow  +  numNodes 
Seleccion. FormacCondicions. Add  Type:=xlCellValue,  Operacor :=xlNocEqual,  _ 
Formulal  :="=0" 

Seleccion . FormacCondicions (Seleccion . FormacCondicions . Count ) . SecFitscPrior 
Mich  Seleccion . FormacCondicions ( 1 ) . Incerior 
.PaccernColorlndex  =  xlAucomacic 
.Color  =  15773696 
■TincAndShade  =  0 


0.  I  will  use  ic  if  I  c 
,  lascCol  +  lascCol  +  4  -i 


in  gee  BVa  co  call  CPLEX. 
numNodes) ) .Selecc 


92 


constraint  .data  (1  of  11) 


File  Edit  View  Insert  Format  Debug  Run  Tools  Add-Ins  Window  Help 

lElB'H  j>  -j  &  *1  P  >  li  a  ^  g  tir  ® _ 

i  |£  [►  oi  a  j  4i  ti  ^  i  i  i  w  &  ;  t.  %  a*  i  «J 
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Project  -  VBAProject 
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i  ^  ASAPUtilitiesRibbon  (ASAP  Utilitie  * 

El  ^gS  VBAProject  (20150202_input  data 

□•••Sj  Miaosoft  Excel  Objects 
BP]  Sheet  1  (Sheet  1) 

§§]  Sheet  10  (weighting) 
m  Sheet  11  (A  1) 

■]  Sheetl2  (Bl) 

H]  Sheetl3  (Cl) 

■3  Sheetl4(Dl) 
m  SheetlS  (El) 

Sf)  Sheet  16  (FI) 
d  Sheetl7(Gl) 
d  Sheet2  (A) 
d  Sheets  (B) 
d  SheeH  (F) 
d  Sheets  (C) 
d  Sheet6  (D) 
d  Sheet7  (E) 
d  Sheets  (G) 

Sheet9  (flow) 

ThisWorkbook 
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;  arc_data 
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«££  cplexvba 
Ht  master 

node_data 


•L 


Properties  -  constraint_data 


|constraint_data  Module 
Alphabetic  I  Categorized  | 


constraint  data 


*1 

"3 


|  (General) 


Option  Explicit 


Counter  for  constraints 
Counter  for  constraints 


(47), 

(47), 


Sub  constraints ( ) 

Dim  i  As  Integer 
Dim  1c  As  Integer  ' 
Dim  c  As  Integer  ' 
Dim  j  As  Integer 
Dim  w  As  Integer 
Dim  p  As  Integer 
Dim  f  As  Integer 
Dim  t  As  Integer 
Dim  uij  As  Integer 


i *******  Finds  the  position  of  the  last  row* 
lastRow  =  Act iveSheet. Cells (Rows. Count, 

'  Malces  the  gap  between  'Arc-specific  data' 
Rows (lastRow  +  3). Select 
With  Selection. Interior 
.Pattern  =  xlSolid 
. PatternColorlndex  =  xlAutomatic 
.Color  =  65535 
.TintAndShade  =  0 
. PatternTintAndShade  =  0 
End  With 

'  Malces  the  'Arc-specific  data'  row  RED. 

Rows(lastRow  +  4). Select 

With  Selection. Interior 
.Pattern  =  xlSolid 
.PatternColorlndex  =  xlAutomatic 
.Color  =  255 
.TintAndShade  =  0 
. PatternTintAndShade  =  0 
End  With 


Writes  the  column  heading  for  the  "Set  Matrix" 
Sheets ("Sheetl") .Select 
Range ("nl2") .Select 

ActiveCell. ForraulaRICl  =  "Constraint  Number" 

Range ("ml3") .Select 

ActiveCell. FormulaRICl  =  "'(47)" 

Range ("nl3") .Select 
ActiveCell. FormulaRICl  =  "'(48)" 

Range ( "ol3") .Select 

ActiveCell . FormulaRICl  =  "'(49)" 

Range ("pl3") .Select 
ActiveCell. FormulaRICl  =  "'(50)" 

Makes  a  grid  system  for  the  Set  matrix 

Range ("ml3 :p"  &  13  +  numNodes) . Select 
With  Selection. Borders 

.LineStyle  =  xlContinuous 
.Colorlndex  =  0 
.TintAndShade  =  0 
.Weight  =  xlThin 
End  With 


'  Places  "0"  in  all  of  the  cells  of  d_ij 
Range (Cells (14,  13),  Cells (13  +  numNodes,  16))  =0 


(48), 

(48), 


(49), 
(49)  , 


and 

and 


(50) 

(50) 


'b") . End (xlUp) .Row 


and  'Constraints'  YELLOW. 
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®Sheet7(E) 
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S]  Sheets  (flow) 

ThisWorkbook 
Modules 
arc_data 
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cplexvba 

master 

A  node_data 


Properties  -  constraint_data 
|constraint data  Module 

Alphabetic  |  Categorized  | 
constraint_data 


T3 


Shades  Che  fir sc  row  for  Sec  macrix. 

Range ("ml3:pl3") .Selecc 
Wich  Seleccion. Incerior 
. Paccern  =  xlSolid 
. PaccernColorlndex  =  xlAucomacic 
.ThemeColor  =  xlThemeColorDarkl 
.TinCAndShade  =  -0.349986266670736 
. PaccernTincAndShade  =  0 
End  Wich 

Wrices  che  column  heading  "Conscraincs" . 

Cells (lascRow  +  4,  2) .Selecc 
AcciveCell . FormulaRICl  =  "Conscraincs" 

Cells (lascRow  +  6,  2) .Selecc 

AcciveCell . FormulaRICl  =  "Conscrainc  Number" 


Creaces  Che  macrix  for  counring  conscraincs  47,  48,  49,  50 
For  i  =  14  To  13  +  numNodes 
■  check  conscrainc  47 

If  (Cells (i,  3) .Value  +  Cells (i,  4) .Value)  =  0  Then 
Cells (i,  13) .Value  =  1 
Else 

'  check  conscrainc  48 

If  (Cells (i,  3) .Value)  =  0  And  (Cells (i,  4) .Value)  >  0  Then 
Cells (i,  14) .Value  =  1 
Else 
End  If 

'check  conscrainc  49 

If  (Cells (i,  3) .Value)  >  0  And  (Cells (i,  4) .Value)  =  0  Then 
Cells (i,  15) .Value  =  1 
Else 
End  If 

'check  conscrainc  50 

If  (Cells (i,  3) .Value)  >  0  And  (Cells(i,  4) .Value)  >  0  Then 
Cells (i,  16) .Value  =  1 
End  If 


Nexc  i 

Sum  for  conscraincs  47,  48,  49,  50  (Sec  macrix) 

Cells  (14  +  numNodes,  13)  =  WorksheecFunccion. Sum (Range (Cells (14, 
Cells (14  +  numNodes,  14)  =  WorksheecFunccion. Sum (Range (Cells (14, 
Cells (14  +  numNodes,  15)  =  WorksheecFunccion. Sum (Range (Cells (14, 
Cells (14  +  numNodes,  16)  =  WorksheecFunccion . Sum (Range (Cells (14, 
consc47  =  Cells (14  +  numNodes,  13) 
consc48  =  Cells (14  +  numNodes,  14) 
consC49  =  Cells (14  +  numNodes,  15) 
consc50  =  Cells (14  +  numNodes,  16) 


Cells (14  +  numNodes, 

Cells (14  +  numNodes, 

15) ,  Cells (14  +  numNodes,  15))) 

16) ,  Cells (14  +  numNodes,  16))) 


13) , 

14) , 


13) )) 

14) )) 


■  ■  Begin  conscrainc  (38)  »•»»•••••■•••••»••■•••••••» 

Wrices  (38)  as  conscrainc  idencifier 
Cells (lascRow  +  7,  3) .Value  =  "'(38)" 

Wrices  Che  LHS  of  conscrainc  (38) . 

Cells (lascRow  +  7,  4)  =  "="  &  Cells(4,  8) .Address (ReferenceScyle :=xlAl, 
RowAbsoluce : =False ,  ColumnAbsoluce : =False ) 


94 


constraint  .data  (3  of  11) 


File  Edit  View  Insert  Fermat  Debug  gun  Iools  Add-Ins  Window  Help 
!£■  ►  JJ  J  «£  ^5  *£  tar  @ 


*1 


i 


!§  ASAPUtilrtiesRibbon  (ASAP_Utilitie 
VBAProject  (20150202_input  data 

B  S3  MiCTOsoft  Excel  Objects 
e  Sheetl  (Sheetl) 

B)  SheetlO  (weighting) 
fi  Sheetll  (Al) 

I  ®Sheetl2(Bl) 

I  ®  Sheet  13  (Cl) 
fi)  Sheet  14  (Dl) 

®  SheetlS  (El) 
fiSheetl6(Fl) 

(Q  Sheet  17  (Gl) 

6]  Sheet2  (A) 

®  Sheet3  (B) 

®  Sheet*!  (F) 

B  Sheets  (C) 

B)  Sheet6  (D) 

@Sheet7(E) 

B)  Sheets  (G) 

6]  Sheet9  (flow) 
t- jtl  ThisWorkbook 
0  £3  Modules 
;■■■■«%  arc_data 
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“3 


Writes  the  RHS  of  constraint  (38) 

■  This  code  [ (cijlRange .Address (ReferenceStyle :=xlAl,  _  RowAbsolute :=False,  ColumnAbsolute :=False) ) ]  turns  ci  j 1 Range  into  its  Excel  n< 
Cells (lastRow  +  7,  6)  =  "=+SOMPRODUCT ( "  £  (cij IRange .Address (ReferenceStyle :=xlAl,  _ 

RowAbsolute: “False,  ColumnAbsolute :=False) )  £  &  (xijlRange. Address (ReferenceStyle :=xlAl,  _ 

RowAbsolute: “False,  ColumnAbsolute:=False) )  £  ")  +SOMPRODOCT ( "  £  (cij2Range. Address ( ReferenceStyle :=xlAl, 
RowAbsolute : “False,  ColumnAbsolute : “False ) )  £  ",n  £  (xij2Range. Address ( ReferenceStyle :—xlAl,  _ 
RowAbsolute: “False,  ColumnAbsolute : “False ) )  £  ■) " 


e.g.  R55:AC66 


1  End  constraint  (38) 


1 '  Begin  constraint  (39)  . . . . 

;  identifier 
-  8,  3) .Value  =  ■■ (39) " 

1  Writes  the  LHS  of  constraint  (39) . 

Cells (lastRow  +  8,  4)  =  "="  £  Cells (4,  9) .Address ( ReferenceStyle :=xlAl, 
RowAbsolute : “False ,  ColumnAbsolute : “False) 

•Cells (lastRow  +  8,  4) .Value  =  Cells(4,  9) .Value 


'  Writes  the  RHS  of  constraint  (39) 

Cells (lastRow  +  8,  6)  =  "“SUMPRODUCT ( "  £  (dijRange. Address ( ReferenceStyle :=xlAl,  _ 

RowAbsolute :=False,  ColumnAbsolute: “False) )  £",("£_ 

(xijlRange. Address (ReferenceStyle :=xlAl,  _ 

RowAbsolute : “False,  ColumnAbsolute: “False) )  £  "+"  £  _ 

(xij2Range. Address (ReferenceStyle :=xlAl,  RowAbsolute :=Fal3e,  ColumnAbsolute : “False) )  £  ' 


1  End  constraint  (39) 


. ,  i . . . .  1 1  1 1  . . .  1 1  1 1 . . .  1 1  1 1 1 . .  Begin  constraint  (40*1) 

•  Writes  (40*1)  as  constraint  identifier 

Cells  (lastRow  +  7,  9)  .Value  =  ’"(40*1)" 

1  Writes  the  LHS  of  constraint  (40*1)  . 

Cells (lastRow  +  7,  10)  =  "="  £  Cell3(4,  10) .Address (ReferenceStyle:=xlAl, 
RowAbsolute : “False ,  ColumnAbsolute : “False ) 

•  Writes  the  “  for  the  constraint  (40*1) 

Cells (lastRow  +  7,  11) .Value  =  "■=>=" 


’  Writes  the  RHS  of  constraint  (40*1) 

Cells (lastRow  +  7,  12)  =  "= (1/B9) *SOMPRODOCT (-  _ 

£  (dijRange. Address (ReferenceStyle :=xlAl,  RowAbsolute : “False,  ColumnAbsolute :=False) )  _ 

£  ”,"  £  (xi j IRange. Address (ReferenceStyle :=xlAl,  RowAbsolute :=False,  ColumnAbsolute : “False ) ) 
£  ")" 


End  constraint  (40*1) 


. . . .  Begin  constraint  (40*2) 

■  Writes  (40*2)  as  constraint  identifier 
Cells (lastRow  +  8,  9) .Value  =  "•(40*2)" 
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®Sheet2(A) 
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®  Sheet4(F) 
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®Sheet7(E) 
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(P)  Sheet9  (flow) 
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4t 
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’  Writes  the  LHS  of  constraint  (40*2)  . 

Cells (lastRow  +  8,  10)  =  "="  £  Cells (4,  10) .Address (ReferenceStyle:=xlAl, 
RowAbsolute : =False,  ColumnAbsolute : “False) 

■Cells (lastRow  +  8,  10) .Value  “  Cells(4,  10) .Value 


’  Writes  the  RHS  of  constraint  (40*2) 

Cells (lastRow  +  8,  12)  =  "“(1/B10) "SOMPRODUCT ("  _ 

£  (dijRange. Address (ReferenceStyle:=xlAl,  RowAbsolute: “False,  ColumnAbsolute : “False ) )  _ 

£  £  (xij2Range. Address (ReferenceStyle :=xlAl,  RowAbsolute :=False,  ColunmAbsolute:=False) ) 

£  ")" 


1  End  constraint  (40*2) 


■  ■  ■  i  ■  •  ■  ■  ■  ■  . . Begin  constraint  (41)  y_i*s  . """" . . . 

1  Writes  (41)  as  constraint  identifier 

Cells (lastRow  +  7,  15) .Value  =  ",(41)" 

'  Writes  the  LHS  of  constraint  (41) . 

■Range (Cells (14,  8),  Cells(13  +  numNodes,  8)). Select  'Used  to  find  the  column  in  question. 
Cells (lastRow  +  7,  16)  =  "=SumC  £  Range (Cells (14,  8),  Cells (13  +  numNodes,  8))  _ 

.Address (ReferenceStyle:=xlAl,  RowAbsolute:=False,  ColumnAbsolute: “False)  £  ") " 


Properties  -  constraint.data  X| 

|  constraint.data  Module  33 


’  Writes  the  RHS  of  constraint  (41) 

Cells (lastRow  +  7,  18)  =  "="  £  Cells(4,  2) .Address (ReferenceStyle :=xlAl, 
RowAbsolute : “False,  ColumnAbsolute : “False) 

■  ■  ■  •  ■  a  ■  •  •  ■  a  ■  ■  ■  ■  i  ■  ■  •  •  ■  a . End  constraint  (41)  . . .  ■  ■  ■  ■  ■ 


■■■■■■' . a  i  a  a  a  t  i  a  a  a  a  a  a  a  t  a  Begin  Constraint  (42)  y  i*a  . . a  a  a  a  a  a  a  a  a  a  a  a  a  a  a  a 

■  Writes  (42)  as  constraint  identifier 

Cells  (lastRow  +  8,  15)  .Value  = 

’  Writes  the  LHS  of  constraint  (42) . 

Cells (lastRow  +  8,  16)  =  "=Sum("  £  Range (Cells (14,  7),  Cells (13  +  numNodes,  7))  _ 
.Address (ReferenceStyle:“xlAl,  RowAbsolute: “False,  ColumnAbsolute :“False)  £  ") " 

’  Writes  the  =  for  the  constraint  (42) 

Cells (lastRow  +  8,  17) .Value  =  "■<=" 

■  Writes  the  RHS  of  constraint  (42) 

Cells (lastRow  +  8,  18)  -  £  Cells (5,  2) .Address (ReferenceStyle :-xlAl,  _ 

RowAbsolute :=False,  ColumnAbsolute :=False) 

'■'■■'■' . a  a  a  a  a  a  a  a  a  a  a  En(J  constraint  (42)  . . a  a  a  a  a  a  a  a  a  a  a  a  a  a  a  a  a  a  a 


’  . . .  . . .  ■  Begin  constraint  (43)  y_i*dc  "■■■ . a  a  a  a  a  a  a  a  a  a  a  a  a  I 

'  Writes  (43)  as  constraint  identifier 

Cells  (lastRow  +  9,  15)  .Value  =  ’"HS)" 

■  Writes  the  LHS  of  constraint  (43) . 

Cells (lastRow  +  9,  16)  =  "“SumC  £  Range (Cells (14,  9),  Cells (13  +  numNodes,  9))  _ 
.Address (ReferenceStyle :=xlAl,  RowAbsolute :=False,  ColumnAbsolute : “False)  £  ")" 


96 


constraint  .data  (5  of  11) 
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I-®  ASAPUtilrtiesRibbon  (ASAPJJtilitie  - 
B  ^  VBAProject  (20150202 .input  data 

El  Microsoft  Excel  Objects 
\m  Sheet  1  (Sheet  1) 

B§3  Sheet  10  (weighting) 
m  Sheetll  (Al) 

13  Sheetl2  (Bl) 

13  Sheet  13  (Cl) 

1  D  Sheetl4  (Dl) 

13  SheetlS  (El) 

13  Sheetl6  (FI) 
he  Sheetl7(Gl) 

!  fflj  Sheet2  (A) 

!•  ®Sheet3(B) 

!  13  Sheet4  (F) 

!•  13  Sheets  (C) 

!  13  Sheet6  (D) 
i  @Sheet7(E) 

I  13  Sheets  (G) 

;  13  Sheet9  (flow) 

^3  ThisWorkbook 
El  Modules 

arc_data 

<gr  i 

I  *«J  cplexvba 
i  master 
J  node data 


3roperties  -  constraint.data 

2< 

constraint_data  Module 

Alphabetic  |  Categorized  | 

'■  30  constraint_data 

|  (General) 


T3  F 


Writes  the  =  for  the  constraint  (43) 

Cells (lastRow  +  9,  17) .Value  =  "■<=" 

Writes  the  RHS  of  constraint  (43) 

Cell3 (lastRow  +  9,  18)  =  "="  £  Cells(6,  2) .Address (ReferenceStyle :=xlAl, 
RowAbsolute: “False,  ColumnAbsolute : “False) 


■Cells (lastRow  +  9,  18)  =  Cells (6,  2) 

I  I  ■  ■  ■  ■  •  ■  i  I  •  ■  ■  i  •  ■  i  ■  ■  ■  ■  ■  •  ■  i  ■  Encl  constraint  (43)  . . . funiiiinnuuniini 

1 1 ■ 1 1 1 1 1 1 ■ • ■ 1 1 1 ■ i > • ■ i > •  ■  i  ■  Begin  constraint  (44)  . . .  ■  ■  ■  ■ 

For  i  =  1  To  numNodes 

Cells (lastRow  +  9  +  i,  15)  =  (44_"  &  i  6  ")"  1  Writes  (44)  as  a  heading. 

Cells (lastRow  +  9  +  i,  16)  =  "=Sum("  £  Range (Cells (13  +  i,  7),  Cells (13  +  i,  9) 
.Address (Reference5tyle:=xlAl,  RowAbsolute :=Fal3e,  ColumnAbsolute : “False)  _ 
&  ")"  '  Writes  the  LHS  of  constraint  (44). 

Cells (lastRow  +  9  +  i,  17) .Value  =  '  Writes  <=  of  constraint  (44) 

Cells (lastRow  +  9  +  i,  18) .Value  =  1  '  Writes  the  RHS  of  constraint  (44) 

Next  i 

1 . . . .  • ' ' 1 1 ' '  End  constraint  (44)  . . . . . . 


>  i  ■  i  >  i  ■  •  > i i i > 1 1 • i 1 1 • i 1 1 • i 1 1  Begin  constraint  (45) 

Writes  (45)  as  constraint  identifier 
For  i  =  1  To  numNodes 

Cells (lastRow  +  8  +  i,  3)  =  "■ (45_"  £  i  £  ")"  '  Writes  (45)  as  a  heading. 

Cells (lastRow  +  8  +  i,  4)  =  "=Sum("  £  Range (Cells (lastRow  -  2  *  numNodes  -  1,  2  *  lastCol  +  4  +  i) ,  _ 
Cells (lastRow  -  numNodes  -  2,  2  *  lastCol  +  4  +  i) ) .Address (ReferenceStyle :=xlAl,  _ 

RowAbsolute : “False,  ColumnAbsolute :=Fal3e)  £  ")"  1  Writes  the  LHS  of  constraint  (45). 

Cells (lastRow  +  8  +  i,  5) .Value  =  1  Writes  the  =  of  constraint  (45) 

Cells (lastRow  +  8  +  i,  6)  =  "="  £  Cells(13  +  i,  10) .Address (Ref erenceStyle:=xlAl,  _ 

RowAbsolute :=False,  ColumnAbsolute :=False)  £  "+"  £  Cells(13  +  i,  11) .Address (ReferenceStyle :=xlAl,  _ 
RowAbsolute : “False,  ColumnAbsolute :=False)  '  Writes  the  RHS  of  constraint  (45) 

Next  i 

1  1  ' ' 1  ' '  ■  1  1  ' ' 1  ' '  ■  1  1  ' ' 1  ' ' '  1  1  '  End  constraint  (45)  . . . . .  ■  ■  ■  ■ 

i i ■ • t 1 1 • i 1 1 1 > i i • i 1 1 • t  1 1 • t 1 1  Begin  constraint  (46)  • 

For  i  =  1  To  numNodes 

Cells (lastRow  +  9  +  numNodes  +  i,  15)  =  (46_"  £  i  £  ")"  1  Writes  (46)  as  a  heading. 

Cells (lastRow  +  9  +  numNodes  +  i,  16)  =  "="  £  Cells (13  +  i,  10) .Address (ReferenceStyle :=xlAl,  _ 

RowAbsolute:=False,  ColumnAbsolute : “False)  ‘  Writes  the  LHS  of  constraint  (46). 

Cells (lastRow  +  9  +  numNodes  +  i,  17) .Value  =  "*<="  1  Writes  the  <=  for  constraint  (46). 

Cells (lastRow  +  9  +  numNodes  +  i,  18)  =  "="  £  Cells  (13  +  i,  5) .Address (ReferenceStyle :=xlAl,  _ 

RowAbsolute: “False,  ColumnAbsolute : “False)  £  £  Cells(13  +  i,  9) .Address  _ 

(ReferenceStyle :=xlAl,  RowAbsolute :=False,  ColumnAbsolute :=False)  '  Writes  the  LHS  of  constraint  (46). 

Next  i 


End  constraint  (46) 
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Bsheet4(F) 
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I  B  Sheet8  (G) 
j  - B  Sheet9  (flow) 
fi  ThisWorkbook 
B  S  Modules 
:  arc_data 

■#EZ . i 

i  it)  cplexvba 
j  master 

node  .data 


Begin  constrain!  (47) 


For  i  =  1  To  numNodes 

If  Cells (13  +  i,  13)  =  1  Then 

Cells (lastRow  +  numNodes  +  9  +  i  -  k,  3)  =  "(47_"  sis")"  ■  Writes  (47_i)  as  a  constraint  identifier. 

Cells (lastRow  +  numNodes  +  9  +  i  -  k,  4)  “  £  Cells (13  +  i,  11) .Address (ReferenceStyle:=xlAl,  _ 

RowAbsolute: “False,  ColumnAbsolute:“False)  '  Writes  the  LHS  of  constraint  (47) . 

Cells (lastRow  +  numNodes  +  9  +  i  -  k,  5)  “  '  Writes  the  ■  for  the  constraint  (47) 

Cells (lastRow  +  numNodes  +  9  +  i  -  k,  6)  *  "»sub("  £  Range (Cells (lastRow  -  2  *  numNodes  -2+k+c,  2*  lastCol  +  S)  _ 

,  Cells (lastRow  -  2  •  numNodes  -  2  +  k  +  e,  2*  lastCol  +  4  +  numNodes))  _ 

•Address (ReferenceStyle:=xlAlr  RowAhsolute: “False,  ColumnAbsolute:=False)  &  ")"  ’  Writes  the  RH5  of  constraint  (47) 


End  constraint  (47) 


For  i  ■  1  To  numNodes 

If  Cells (13  +  i,  14)  -  1  Then 

Cells (lastRow  +  2  •  numNodes  +  10  +  i  -  k,  15)  “  "(48_"  6  i  &  ")"  'Writes  (48_i)  as  a  constraint  identifier. 

Cells (lastRow  +  2  "  numNodes  +  10  +  i  -  k,  16)  “  "“sum("  &  Range (Cells (lastRow  -  2  *  numNodes  -  2  +  i,  2  *  lastCol  +  5) ,  _ 

Cells (lastRow  -  2  "  numNodes  -  2  +  i,  2  "  lastCol  +  4  +  numNodes))  _ 

•Address (ReferenceStyle:=xlAl,  RowAhsolute :=False,  ColumnAhsolute: “False)  £  ")"  '  Writes  the  LHS  of  constraint  (48) 

Cells (lastRow  +  2  «  numNodes  +  10  +  1  -  k,  17)  =  "<=”  *  Writes  the  <=  for  the  constraint  (48) 

Cells  (lastRow  +  2  *  numNode3  +  10  +  1  -  k,  18)  =  "="  &  Cells(13  +  1,  11)  .Address  (Ref  erenceStyle:=xlAl,  RowAhsolute:  “False,  ColumnAhsolute:  “False) 

£  "+"  £  Cells(13  +  i,  4) .Address (ReferenceStyle:=xlAi,  RowAhsolute: “False,  ColumnAhsolute: “False) 

£  £  Cells (13  +  i,  8)  .Address (ReferenceStyle:=xlAl,  RowAhsolute: “False,  ColumnAhsolute: “False)  ■  Writes  the  RHS  of  constraint  (48) 


(49) 


For  i  =  1  To  numNodes 

If  Cells (13  +  1,  15)  =  1  Then 

Cells (lastRow  +  2  *  numNodes  +  10  +  const48  +  i  -  k,  15)  =  ” (49_”  £  i  £  ")"  'Writes  (49_i)  as  a  constraint  identifier. 

Cells  (lastRow  +  2  *  numNodes  +  10  +  const48  +  i  -  k,  16)  =  "=sumC  6  Range  (Cells  (lastRow  -  2  *  numNodes  -  2  +  i,  2  *  lastCol  +  5),  _ 

Cells  (lastRow  -  2  *  numNodes  -  2  +  i,  2  *  lastCol  +  4  +  numNodes))  _ 

•  Address  (ReferenceStyle:=xlAl,  RowAhsolute :=False,  ColumnAhsolute :=False)  £  ")  ”  •  Writes  the  LHS  of  constraint  (49) 

Cells  (lastRow  +  2  *  numNodes  +  10  +  const48  +  i  -  k,  17)  =  "<=■  ■  Writes  the  <=  for  the  constraint  (49) 

Cells (lastRow  +  2  *  numNodes  +  10  +  const48  +  i  -  k,  18)  =  "="  £  Cells(13  +  i,  11) .Address (ReferenceStyle:=xlAl,  RowAbsolute:=False,  ColumnAhsolute :=False) 
£  "+"  £  Cells (13  +  i,  3)  .Address (ReferenceStyle:=xlAl,  RowAhsolute: “False,  ColumnAhsolute: “False)  _ 

£  £  Cells(13  +  i,  7) .Address (ReferenceStyle:=xlAl,  RowAhsolute :=False,  ColumnAhsolute: “False)  '  Writes  the  RHS  of  constraint  (49) 
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e  Sheetl6  (FI) 

!  ®  Sheet  17  (Gl) 
m  Sheet2 (A) 
m  Sbeet3  (B) 

@Sheet4(F) 

®Sheet5(C) 
m  Sheets  (D) 
fi  Sheet7 (E) 

%  Sheets  (G) 
j  ®  Sheets  (flow) 
ffl  ThisWorkbook 
0  <S  Modules 


. . . . .  . . . .  Begin  constraint  (52)  . . . . . . 

■  Writes  (52)  as  constraint  identifier 

Cells (lastRow  +  10  +  2  *  numNodes  +  const 4 8  +  const49  +  constSO,  15) .Value  =  "'(52)" 

Cells (lastRow  +  12  +  2  *  numNodes  +  const48  +  const49  +  const50,  15) .Value  =  "Constraint  (52)” 

•  Writes  the  LHS  of  constraint  (52) . 

'  Places  the  row  designation  for  x^S'l  +  x_i3“2  matrix 
3=1 

For  i  =  (lastRow  +  13  +  2  *  numNodes  +  const48  +  const49  +  constSO  +1)  To  _ 

(lastRow  +  13  +  2  *  numNodes  +  const48  +  const49  +  constSO  +  numNodes) 

Cells (i,  15)  -  j 

3  =  3+1 

Next  i 

'  Places  the  column  designation  for  the  x_i3~l  +  x_ij~2  matrix 
Ic  =  1 

For  i  =  16  To  numNodes  +  15  'lastCol  +  numNodes 

Cells (lastRow  +  12  +  2  *  numNodes  +  const48  +  const49  +  constSO  +  1,  i)  =  k 


Next  i 


arc_data 


«g|  ~~1 

cplexvba 

master 

node_data 


'  Makes  a  grid  system  for  the  x_i3"l  +  x_i3'"2  matrix 

Range (Cells (lastRow  +  12  +  2  *  numNodes  +  const48  +  const49  +  constSO  +  1,  15),  _ 

Cells (lastRow  +  12  +  2  *  numNodes  +  const48  +  const49  +  constSO  +  1  +  numNodes,  15  +  numNodes) ) .Select 
With  Selection. Borders 

.LineStyle  =  xlContinuous 
■Color Index  =  0 
.TintAndShade  =  0 
.Weight  =  xlThin 
End  With 

■  Shades  the  first  row  of  the  x_i3“l  +  x_i3~2  matrix 

Range (Cells (lastRow  +  12  +  2  *  numNodes  +  const48  +  const49  +  constSO  +  1,  15),  _ 

Cells (lastRow  +  12  +  2  *  numNodes  +  const48  +  con3t49  +  constSO  +  1,  numNodes  +  15)). Select 

With  Selection. Interior 
•Pattern  =  xlSolid 
. PatternColorlndex  =  xlAutomatic 
.ThemeColor  =  xlThemeColorDarkl 
.TintAndShade  =  -0.349986266670736 
. PatternTintAndShade  =  0 
End  With 

’  Shades  the  first  column  of  the  x_i3“l  +  x_i3“2  matrix 

Range (Cells (lastRow  +  12  +  2  *  numNodes  +  const48  +  const49  +  constSO  +  1,  15),  _ 

Cells (lastRow  +  12  +  2  *  numNodes  +  const48  +  const49  +  const50  +  1  +  numNodes,  15)). Select 

With  Selection. Interior 
.Pattern  =  xlSolid 
. PatternColorlndex  =  xlAutomatic 
.ThemeColor  =  xlThemeColorDarkl 
.TintAndShade  =  -0.349986266670736 
.PatternTintAndShade  =  0 
End  With 

1  Places  "0"  in  all  of  the  cells  of  the  x_i3/'l  +  x_i3/'2  matrix 

Range (Cells (lastRow  +  13  +  2  *  numNodes  +  const48  +  const49  +  const50  +  1,  16),  _ 

Cells (lastRow  +  12  +  2  *  numNodes  +  const48  +  const49  +  constSO  +  1  +  numNodes,  15  +  numNodes))  =  0 

1  Finds  the  position  of  the  last  row  and  last  col  (lastRow2  S  lastCol2) 
lastRow2  =  ActiveSheet. Cells (Rows. Count,  "o") .End (xlDp) .Row 
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■  Copies  the  u_ij  matrix  that  was  already  create 
uij  Range .Select 
Selection . Copy 

Range (Cells (lastRow2  -  numNodes  +  1,  lastCol2  h 
ActiveSheet . Paste 


;  datal  module 


Cells  (lastRow2,  lastCo!2  h 


k  numNodes) ) .Select 


reads  the  constraints  i 

•CPLEX  constraints 

Dim  eqConst  As  Range 
Dim  eqConstBound  As  Range 
Dim  geqConst  As  Range 
Dim  geqConstBound  As  Range 
Dim  leqConst  As  Range 
Dim  leqConstBound  As  Range 
Dim  numLeqConst  As  Integer 
Dim  numEqConst  As  Integer 
Dim  const52  As  Range 
Dim  const52bound  As  Range 
Dim  constRow  As  Integer 
Dim  const52ind  As  Integer 

constRow  =  lastRow  +  7 


Set  eqConst  =  Range (Cells (constRow,  4),  Cells (constRow  +  numEqConst  -  1,  4)) 

Set  eqConstBound  =  Range (Cells (constRow,  6),  Cells ( constRow  +  numEqConst  -  1,  6)) 

Set  geqConst  =  Range (Cells (constRow,  10) ,  Cells (constRow  +  1,  10) ) 

Set  geqConstBound  =  Range (Cells (constRow,  12),  Cells (constRow  +  1,  12)) 

Set  leqConst  =  Range (Cells (constRow,  16),  Cells (constRow  +  numLeqConst  -  1,  16)) 

Set  leqConstBound  -  Range (Cells (constRow,  18),  Cells (constRow  +  numLeqConst  -  1,  18)) 

const52ind  =  constRow  +  numLeqConst  -1+5 


CPXaddConstraint  constraint :=eqConst,  Lb : “eqConstBound,  Db : “eqConstBound 
CPXaddConstraint  constraint: “leqConst,  Db: “leqConstBound 
CPXaddConstraint  constraint: “geqConst,  Lb : “geqConstBound 
CPXaddConstraint  constraint: =const52,  Db:=const52bound 

'  '  '  '  . .  I  might  need  to  use  the  code  below  to  copy  the  u_ij  matrix 


.  16  +  numNodes  - 


v  numNodes  -  1) ) 


datal  module 


1  Places  the  u_i j ' s  in  the  appropriate  cells 
for  w  =  p  To  (lastRow  -  numNodes  -  4) 
f  =  Cells  (w,  1) 
t  -  Cells (w,  2) 
uij  =  Cells  (w,  6) 

For  i  =  (firstRow  +  1)  To  (numNodes  +  firstRow)  '  This  checks  all  the  row  values. 

For  j  =  3  To  (numNodes  +  2)  ’  This  checks  all  the  column  values. 

If  (Cells (i,  2)  “  f  And  Cells (firstRow,  j)  “  t)  Then 
Cells (i,  j)  =  uij 

End  If 
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constraints 


1 CPLEX  constraints 

Dim  eqConst  As  Range 
Dim  eqConstBound  As  Range 
Dim  geqConst  A3  Range 
Dim  geqConstBound  As  Range 
Dim  leqConst  As  Range 
Dim  leqConstBound  As  Range 
Dim  numLeqConst  As  Integer 
Dim  numEqConst  As  Integer 
Dim  const52  As  Range 
Dim  const52bound  As  Range 
Dim  constRow  As  Integer 
Dim  const52ind  As  Integer 

constRow  =  lastRow  +  7 

numEqConst  =  2  +  2  *  numNodes  +  Cells (numNodes  +  14,  13)  .Value 

numLeqConst  “  3  +  2  *  numNodes  +  Application. Sum (Range (Cells (numNodes  +  14,  14),  Cells (numNodes  +  14,  16))) 

Set  eqConst  =  Range (Cells (constRow,  4),  Cells (constRow  +  numEqConst  -  1,  4)) 

Set  eqConstBound  =  Range (Cells (constRow,  6),  Cells (constRow  +  numEqConst  -  1,  6)) 

Set  geqConst  =  Range (Cells (constRow,  10),  Cells (constRow  +  1,  10)) 

Set  geqConstBound  =  Range (Cells (constRow,  12),  Cells (constRow  +  1,  12)) 

Set  leqConst  =  Range (Cells (constRow,  16),  Cells (constRow  +  numLeqConst  -  1,  16)) 

Set  leqConstBound  =  Range (Cells (constRow,  18),  Cells (constRow  +  numLeqConst  -  1,  18)) 

const52ind  =  constRow  ^ 


numLeqConst  - 


CPXaddConstraint  constraint 
CPXaddConstraint  constraint 
CPXaddConstraint  constraint 
CPXaddConstraint  constraint : 
i  t  ■  i  ■  ■  i  ■  ■  ■  i  •  i  i  •  ■  •  j  might  need 


=eqConst,  Lb:=eqConstBound,  Ob : “eqConstBound 
“leqConst,  Ob : “leqConstBound 
“geqConst,  Lb: “geqConstBound 
=const52,  Ob:=const52bound 

the  code  below  to  copy  the  u_ij  matrix  from  the  c 


:  datal  module 


'  Places  the  u_i j ' s  in  the  appropriate  cells 
For  w  =  p  To  (lastRow  -  numNodes  -  4) 
f  =  Cells (w,  1) 
t  =  Cells (w,  2) 
uij  =  Cells (w,  6) 

For  i  =  (firstRow  +  1)  To  (numNodes  +  firstRow)  ■  This  checks  all  the  i 
For  j  =  3  To  (numNodes  +  2)  '  This  checks  all  the  column  values. 

If  (Cells (i,  2)  =  f  And  Cells (firstRow,  j)  =  t)  Then 
Cells (i,  j)  =  ui) 

End  If 
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»J jcPXupdateVariable 


End  Function 


Remove  variable (s)  at  1-based  index  Idx  from  the  model’ 3  variable  list. 

Parameters 

Idx  —  The  1-based  index  of  the  variable (s)  to  be  removed  (as  returned  by  CPXaddVariable) . 
If  the  argument  is  omitted  all  variables  are  removed  from  the  model. 

Returns 

True  on  success  and  false  otherwise. 

Function  CPXremoveVariable (Optional  Idx  As  Variant)  As  Boolean 

CPXremoveVariable  =  Application. Run ("CPLEX.REMOVEVARIABLE",  Idx) 

End  Function 


leaving  them 


i  the  model ’ s 


Add  a  constraint  to  the  model  on  the  active  sheet. 

All  arguments  but  Range  are  optional  and  leaving  them  out  is  the  s 
empty  in  the  GUI. 

Upon  success  the  function  returns  the  1-based  index  of  the  new  constraint 
constraint  list.  This  index  can  be  used  as  input  for  CPXupdateConstraint  or 
CPXremoveConstraint . 

Parameters 

Range  —  The  cells  that  are  to  be  constrained.  This  must  not  be  an  empty  range. 

Lb  —  The  lower  bound (s)  for  the  constrained  cells.  Leaving  this  out  means 

"no  lower  bound" . 

Ob  —  The  upper  bound (s)  for  the  constrained  cells.  Leaving  this  out  means 

"no  upper  bound" . 

Returns 

The  1-based  index  of  the  newly  created  constraint  on  success  or  0  if  the 
constraint  could  not  be  added. 

Function  CPXaddConstraint (constraint  A3  Variant,  Optional  Lb  As  Variant,  Optional  Ob  As  Variant)  As  Long 
CPXaddConstraint  =  Application. Run ("CPLEX.ADDCONSTRAINT",  constraint.  Lb,  Ob) 

End  Function 


Opdate  constraint  at  1-based  index  Idx  in  the  model’s  constraint  list. 

All  arguments  but  Range  are  optional  and  leaving  them  out  is  the  same  as  leaving  them 
empty  in  the  GDI. 

Note:  The  function  will  change  ALL  fields  of  the  constraint  definition  —  even  those 
for  which  you  leave  out  parameters. 

Parameters 

Idx  —  The  1-based  index  of  the  constraint  to  be  changed  (as  returned  by  CPXaddConstraint) . 

Range  —  The  cells  that  are  to  be  constrained.  This  must  not  be  an  empty  range. 

Lb  —  The  lower  bound (s)  for  the  constrained  cells.  Leaving  this  out  means 

"no  lower  bound" . 

Ob  —  The  upper  bound (s)  for  the  constrained  cells.  Leaving  this  out  means 

"no  upper  bound" . 

Returns 

True  on  success  and  false  otherwise . 

Function  CPXupdateConstraint (Idx  As  Variant,  constraint  As  Variant,  Optional  Lb  As  Variant,  Optional  Ob  As  Variant)  As  Boolean 
CPXupdateConstraint  =  Application. Run ("CPLEX.OPDATECONSTRAINT",  Idx,  constraint.  Lb,  Ob) 

End  Function 


Remove  constraint  at  1-based  index  Idx  from  the  model ■ s  variable  list . 

Parameters 

Idx  —  The  1-based  index  of  the  constraint  to  be  removed  (as  returned  by  CPXaddConstraint) . 
If  the  argument  is  omitted  all  constraints  are  removed  from  the  model. 

Returns 

True  on  success  and  false  otherwise. 

Function  CPXremoveConstraint (Optional  Idx  As  Variant)  As  Boolean 

CPXremoveConstraint  =  Application. Run ("CPLEX.REMOVECONSTRAINT",  Idx) 

End  Function 
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LinOrQuad  —  Flag  to  control  whether  to  display  a  warning  message  if  CPLEX 

found  an  unknown  function  and  replaced  that  with  its  presumed  linear 
or  quadratic  equivalent. 

ExportModel  —  Flag  to  control  whether  the  optimization  model  is  to  be  exported 


True  on  success  and  false  otherwise. 

Function  CPXsetSpecial (Optional  Stoplnt  As  Variant,  Optional  LinOrQuad  As  Variant,  Optional  ExportModel  As  Variant,  Optional  ModelFile  As  Variant)  As  Boolean 
CPXsetSpecial  =  Application. Run("CPLEX.SETSPECIAL",  Stoplnt,  LinOrQuad,  ExportModel,  ModelFile) 


Note:  The  function  does  not  check  foi 
of  the  model. 

The  index  returned  by  this  function  c 


duplicate  parameters  in  the  parameter  list 
in  be  used  as  argument  to  CPXupdateParameter  c 


Returns 

The  1-based  index  of  the  newly  created  parameter  in  the  model's  parameter  list  c 
success  and  0  on  failure. 

motion  CPXaddParameter (Number  As  Variant,  Value  As  Variant)  As  Boolean 
CPXaddParameter  =  Application. Run ("CPLEX. ADDPARAMETER",  Number,  Value) 


Update  parameter  a 
Parameters 
Idx  --  The  1-ba 
Number  —  The  numb 


1-based  index  Idx. 


rr.ed  by  CPXaddPai 


True  on  success  and  false  otherwise. 

Function  CPXupdatePar ameter ( Idx  As  Variant,  Number  As  Variant,  Value  As  Variant)  As  Boolean 
CPXupdateParameter  =  Application. Run ("CPLEX. UPDATEPARAMETER",  Idx,  Number,  Value) 


!  returned  by  CPXaddParameter) . 


( "CPLEX.  REMOVEPARAMETER",  Idx) 


n  the  active  sheet. 
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CPXaddParajnecer  =  Application. Run ("CPLEX.ADDPARAMETER",  Number,  Value) 
End  Function 


'  Update  parameter  at  l-based  index  Idx. 

'  Parameters 

'  Idx  —  The  l-based  index  of  the  parameter  to  be  changed  (as  returned  by  CPXaddParameter) . 
*  Number  —  The  number  or  name  of  the  parameter  to  be  set. 

'  Value  —  The  value  to  be  set  for  parameter  Number. 

'  Returns 

'  True  on  success  and  false  otherwise. 

Function  CPXupdateParameter (Idx  As  Variant,  Number  As  Variant,  Value  As  Variant)  As  Boolean 
CPXupdateParameter  =  Application. Run ("CPLEX.UPDATEPARAMETER",  Idx,  Number,  Value) 

End  Function 

Remove  parameter  at  l-based  index  Idx  from  parameter  list. 

Idx  —  The  l-based  index  of  the  parameter  to  be  removed  (as  returned  by  CPXaddParameter) . 
If  omitted  all  parameters  are  removed. 

Returns 

True  on  success  and  false  on  failure. 

Function  CPXremoveParameter (Idx  As  Variant)  As  Boolean 

CPXremoveParameter  =  Application. Run ("CPLEX.REMOVEPARAMETER",  Idx) 

End  Function 


Set  the  objective  function  for  the  model  on  the  active  sheet. 

Target  is  only  optional  if  Sense  is  not  3. 

Parameters 

ObjCell  —  The  objective  function  cell.  This  must  be  a  single  cell. 

Sense  —  The  sense  of  the  objective  function.  Allowed  values  are 

1  -  Maximize  ObjCell 

2  -  Minimize  ObjCell 

3  -  Find  a  solution  so  that  ObjCell  equals  Target 
Target  —  The  target  value  for  Sense  =  3. 

Returns 

Function  CPXsetObjective (objCell  As  Variant,  Sense  As  Integer,  Optional  Target  As  Double)  As  Boo: 

CPXsetObjective  =  Application. Run("CPLEX.SETOBJECTIVE",  objCell,  Sense,  Target) 

End  Function 


Invoices  CPLEX  on  the  model  defined  on  the  active  sheet. 

Note:  If  this  function  returns  true  that  does  not  mean  that  a  solution  was  found. 

CPLEX  may  also  have  successfully  proved  infeasibility  of  the  model. 

NoDialog  —  If  a  true  value  is  passed  the  final  dialog  will  not  be  shown.  In  this  case 
there  is  no  way  to  reset  the  variable  cells  to  their  original  values 
or  perform  a  solution  or  sensitivity  analysis. 

If  the  value  passed  cannot  be  converted  to  a  boolean  it  is  ignored. 

Returns 

True  on  success  and  false  on  failure. 

Function  CPXsolve (Optional  NoDialog  As  Variant)  As  Boolean 
CPXsolve  =  Application. Run ("CPLEX. SOLVE",  NoDialog) 

End  Function 


1  Displays  the  CPLEX  Excel  connector's 
Sub  CPXdialog ( ) 

Application. Run  ("CPLEX. DIALOG") 
End  Sub 


solve  dialog. 
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The  Military  Theater  Distribution 
Network  Design  Problem 


o 

00 


Objectives: 

■  Develop  a  math  programming  model 
to  locate  aerial/sea  ports  of 
debarkation  (A/SPOD)  and  distribution 
centers  (DC)  in  order  to  maximize 
multi-commodity  flows  in  support  of 
contingency  operations. 

■  Develop  a  simple-to-use-and- 
understand  user  input-process. 

•  Demonstrate  the  model  by  using  a 
scenario-based  case  study. 

Background: 

•  The  movement  of  supplies  into  an 
Area  of  Operations  is  a  joint 
endeavor. 

•  The  US  Army  is  responsible  for  the 
ground  movement  of  supplies  from 
A/SPODs  to  the  force  in  the  field. 

■  The  Army  distribution  hierarchy 
supports  the  operational  hierarchy 
to  distribute  supplies. 


Operational  Hierarchy 

JTF 


Literature  Review  pertaining  to: 

■  Facility  Location  Models 

•  Network  Flow  Models 

•  Location  and  Routing  Models 


MAJ(P)  Robert  Craig 
Advisor:  LTC  Brian  J.  Lunday,  PhD 
Readers:  Dr.  Ray  R.  Hill,  PhD,  MAJ  Peter  Nesbitt 

Department  of  Operational  Sciences  (ENS) 

Air  Force  Institute  of  Technology 


Sponsor: 

TRAC-LEE 


Methodology: 

•  Move  supplies  through  a  transportation 
network  as  Mixed  Integer  Linear  Program 

•  Locate  A/SPODs  and  DCs  with  multi¬ 
objective  approach  to  meet  FOB  demands 

•  The  multi-objective  approach  includes: 

•  Risk  to  supplies  being  transported. 

•  Total  distance  travelled  by  supplies 

•  Maximum  per  capita  workload 
supported  by  transportation  assets  at  a 
given  echelon 

•  Classify  flow  as  either  A/SPOD  to  DC 
(Echelon  1)  or  DC  to  FOB  (Echelon  2) 

•  Use  echelon  flow  to  identify  routing 

MILP  Formulation 


h~& LS'W-  U  Metrics 


Network: 

•  Nodes:  135  road  junctions  across 
Nigeria 

•  Arcs:  414  arcs  which  are  the  roads 
connecting  the  135  nodes 

•  Arcs:  Classified  by  their  supply 
throughput  capacity 


Optimal  Solution  for  Minimal  Risk 


Bounds  on 
facilities 


Flow  balance 
constraints 


Flow  bounds 

Variable 

restrictions 


Optimal  Solution  for  Maximum 
per  capita  workload 


Optimal  Solution  for  Minimizing 
Total  Supply  Distance 


Conclusions: 

•  The  model  handles  multiple  objectives 
and  user  priorities  to  locate  A/SPODs 
and  DCs 

•  The  simple-to-use-and-understand 
user  input-process  solves  realistic 
sized  instances  within  two  minutes 

•  The  model  uses  a  weighting  scheme 
to  allow  decision  maker  influence 


Appendix  C.  Story  Board 
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